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APPLICATION  OF  RESEARCH  TO  THE  NEEDS 
OF  THE  U.S.  NAVY 


Mechanical  face  seals  are  used  in  numerous  applications  in  Naval 
machinery.  These  applications  range  from  propeller  shaft  seals  to 
boiler  feed  pump  seals.  In  such  equipment  the  mechanical  seal  plays 
a vital  role.  When  such  seals  fail,  repair  is  costly  both  in  terms 
of  lost  time  and  direct  cost,  so  any  improvement  in  seal  life  and 
reliability  would  be  of  significant  benefit. 

As  more  advanced  equipment  is  designed,  it  is  sometimes  difficult 
to  achieve  desired  performance  in  more  severe  service  environments 
with  the  present  state  of  the  art  of  seal  design.  Thus,  an  improveir«nt 
in  seal  technology  would  serve  this  important  application. 

The  immediate  objective  of  the  research  herein  is  to  further  the 
understanding  of  mecheuiical  face  seal  lubrication  phenomena.  The  ulti- 
mate objective  is  to  develop  the  capability  of  designing  contacting 
face  seals  having  a longer  life,  greater  reliability,  and  for  extreme 
environments.  Thus,  the  objectives  of  this  research  are  compatible 
with  mechanical  face  seal  needs  for  N^val  machinery. 
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CHAPTER  1 


LUBRICATION  IN  MECHANICAL  FACE  SEALS 

INTRODUCTION 

In  applications  where  a rotating  shaft  must  pass  from  one  fluid 
region  to  another,  contacting  mechanical  face  seals*  play  the  essential 
role  of  minimizing  the  transfer  of  fluid  between  the  regions.  Appli- 
cations of  face  seals  range  from  water  pump  seals  to  process  pump  seals 
to  propeller  shaft  seals. 

The  performance  and  reliability  of  contacting  mechanical  face  seals 
are  of  great  importance  for  any  type  of  equipment  where  minimal  leakage, 
high  reliability,  and  long  life  cire  necessary.  Even  for  equipment  where 
these  factors  are  not  so  critical,  seal  failures  and  short  seal  life 
lead  to  high  operating  cost  due  to  down  time  and  maintenance  cost. 

Even  though  mechanical  face  seal  technology  has  been  steadily  im- 
proving over  the  past  several  decades,  further  improvement  in  the  state 
of  the  art  of  seal  design  would  be  most  beneficial.  Although  seals 
having  an  acceptable  life  and  reliability  can  be  designed  for  many  ap- 
plications, further  improvement  in  seal  life  and  reliability  would  re- 
sult in  significant  cost  savings  to  the  user.  Also,  there  are  numerous 
mechanical  face  seal  applications  where  seal  loading,  reliability,  life, 
and  leakage  requirements  are  difficult  to  achieve  within  the  present 
state  of  the  art.  Examples  of  such  applications  are  seals  for  pumps 
for  nuclear  power  plants  and  seals  for  large  diameter  submarine  pro- 
peller shafts.  Additionally,  the  friction  losses  in  face  seals  repre- 
sent a significant  fraction  of  energy  consumed  for  pumping  purposes. 
Within  the  present  state  of  the  art  of  seal  design,  it  is  very  difficult 
to  design  a low  leakage  seal  that  also  has  a low  friction  loss. 

The  main  barrier  to  the  advancement  of  the  state  of  the  art  is  that 
the  mechanics  of  seal  operation  are  not  well  enough  understood  to  be 


*The  class  of  low  leakage  face  seals  where  there  is  definite  contact  and 
wear  of  the  faces  as  opposed  to  hydrostatic  or  hydrodynamic  where  a 
definite  clearance  is  maintained. 


1 


able  to  reasonably  anticipate  seal  performance  as  a function  of  design 
parameters.  There  are  no  well  established  fundamental  theoretical 
bases  that  can  be  used  to  indicate  the  type  of  seal  design  that  will  give 
improved  performance.  Improvements  that  have  been  made  have  been  brought 
about  largely  by  trial  eind  error  combined  with  elementary  sealing  theories. 

In  order  to  be  able  to  predict  the  performance  of  contacting  face 
seals  as  a function  of  design  parameters,  it  is  essential  that  the  lubri- 
cation mechanisms  between  the  faces  be  well  understood.  At  present  it  is 
known  that  hydrodynamic  lubrication  plays  some  role  in  providing  load 
support  for  oil  seals  as  well  as  water  seals.  But,  the  precise  nature 
of  this  lubrication  is  not  known.  Several  theories  have  been  put  forth. 
However,  as  discussed  later,  these  theories  have  not  been  verified  for 
contacting  face  seals,  and  it  is  not  possible  as  yet  to  use  these  theories 
for  the  design  of  contacting  face  seals. 

An  argument  can  be  made  for  abandoning  contacting  face  seal  designs 
and  using  hydrodynamic  or  hydrostatic  non-contacting  designs  for  which 
theories  are  readily  available.  However,  there  are  many  applications 
where  non-contacting  face  seals  cannot  be  used  because  of  high  cost,  high 
leakage,  or  inability  to  control  the  geometry.  The  contacting  face  seal 
offers  low  leakage  at  relatively  low  cost  compared  to  non-contacting  seals. 
Therefore,  an  improvement  in  the  performance  of  contacting  face  seals 
while  retaining  their  advantages  is  a worthwhile  objective. 

In  this  report,  a two  dimensional  theory  of  hydrodynamic  lubrication 
with  contact  and  wear  is  examined  and  developed.  The  hydrodynamic  lubri- 
cation considered  is  the  type  caused  by  circumferential  waviness.  This 
particular  lubrication  mechanism  was  chosen  for  study  for  two  reasons. 
First,  as  later  discussed,  there  is  some  experimental  evidence  that  such 
lubrication  occurs  in  seals.  Second,  waviness  is  a controllable  param- 
eter, and  lubrication  might  therefore  be  improved  through  seal  design. 

In  the  model  developed,  the  interaction  of  the  hydrodyncunic  pressure 
with  elastic  distortion  of  the  seal  ring  is  considered.  The  behavior  of 
hydrodynamic  lubrication,  leakage,  and  wear  in  relation  to  controllable 
seal  design  parameters  such  as  ring  stiffness,  material  properties,  and 
roughness  has  been  explored.  These  results  are  used  to  attempt  to  ex- 
plain the  behavior  of  existing  seals  as  well  as  to  suggest  means  of 
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improving  seal  design.  Also  in  this  report,  a comparison  to  some  ex- 
perimental results  has  been  made.  Roughness  models  have  been  evaluated, 
and  test  rig  plans  are  presented. 


In  a previous  report  [1]  and  paper  f2]  a model  similar  to  the  above 
was  developed  on  a one  dimensional  basis.  The  two  dimensional  results 
herein  are  considered  to  be  much  more  accurate  than  the  previous  results 
[1] . The  present  model  accurately  includes  seal  leakages.  This  param- 
eter was  included  in  only  the  most  approximate  manner  in  the  one  dimen- 
sional model.  The  conclusions  drawn  based  on  the  two  dimensional  model 
are  somewhat  different  than  those  based  on  the  one  dimensional  model. 

The  end  result  of  successful  research  into  this  subject  will  be  a 
mathematical  model  that  can  be  used  to  accurately  predict  seal  lubrica- 
tion and  leakage  behavior  as  a function  of  controllable  design  parameters 
and  to  suggest  means  by  which  the  performance  of  contacting  face  seals 
can  be  improved . 

MECHANICAL  FACE  SEAL  LUBRICATION 

The  mechanical  face  seal  consists  basically  of  two  annular  rings 
which  rotate  relative  to  each  other  and  which  are  pressed  together  by 
spring  and  fluid  pressures  (see  figure  1-1) . The  surfaces  that  rub  to- 
gether are  generally  manufactured  as  flat  as  possible  initially  so  as 
to  minimize  leakage.  The  effective  gap  between  the  faces  is  ideally 
quite  small  (order  of  1 pm)  so  that  leakage  flow  across  the  faces  will 
be  quite  small.  The  difficulty  in  designing  a mechanical  seal  is  in 
maintaining  the  gap  at  a very  low  value  while  at  the  same  time  providing 
a definite  lubricant  film  between  the  faces. 

The  load  that  must  be  supported  at  the  faces  of  a mechanical  seal 
is  due  primarily  to  loading  caused  by  the  sealed  pressure.  The  unit 
face  load  can  be  expressed  as  some  fraction  B of  the  sealed  pressure 
where  B can  be  made  greater  or  less  than  unity  by  geometry  selection. 

The  load  support  at  the  faces  is  derived  from  fluid  pressure  and  mechan- 
ical pressure.  If  the  fluid  pressure  at  the  faces  is  large  enough  to 
support  all  of  the  load,  then  there  will  be  no  contact  and  no  adhesive 
* 

Numbers  in  brackets  refer  to  List  of  References  at  the  end  of  report. 


i 


i 

! 


3 


PRESSURE 

SEAL  RING 
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wear.*  If  none  of  the  load  is  supported  by  fluid  pressure/  the  load  must 
be  carried  by  mechanical  contact  and  the  wear  rate  will  be  large. 
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In  practice,  seals  often  operate  at  one  of  two  extremes.  At  one 
extreme,  a large  gap  will  be  created  by  hydrostatic  or  hydrodynamic  ef- 
fects and  the  seal  will  leak  a large  amount.  At  the  opposite  extreme, 
the  gap  will  close  completely  and  leakage  will  be  very  low.  However,  the 
beneficial  lubricant  film  will  be  lost  and  wear  and  heat  generation  will 
increase. 

Based  on  the  above,  it  can  be  concluded  that  an  effective  seal 
should  operate  between  these  two  extremes — having  both  adequate  lubri- 
cation and  low  leakage.  To  do  this  requires  that  any  fluid  pressure 
generation  mechanism  used  to  provide  lift  and  lubrication  to  the  seal 
must  be  very  carefully  controlled.  At  present,  in  commercial  contacting 
face  seals  this  lubrication  is  left  primarily  to  chance,  and  often  such 
seals  operate  at  one  of  the  extremes  mentioned.  Quite  commonly  such 
seals  will  operate  in  the  low  leakage  condition  where  wear  progresses  at 
a definite  rata.  Such  seals  are  quite  satisfactory  for  many  applications. 
However,  when  pressure  or  speed  requirements  are  increased,  then  this 
mode  of  operation  may  lead  to  rapid  failure. 

In  order  to  be  able  to  disign  into  the  seal  the  proper  balance  be- 
tween lubrication  and  leakage,  the  origin  of  lubrication  mechanisms  must 
first  be  fully  understood.  A considerable  amount  of  research  has  been 
accomplished  in  this  area  over  the  years.  Various  theories  have  emerged 
and  are  briefly  described  below. 

PREVIOUS  RESEARCH 

The  first  observation  to  be  made  about  mechanical  face  seals  is  that 
the  hydrostatic  pressure  drop  across  the  faces  provides  some  amount  of 
load  support  in  all  circumstances.  For  parallel  faces,  a unit  load  sup- 
port of  one  half  of  the  sealed  pressure  will  be  provided  by  this  hydro- 
static effect.  If  the  seal  faces  become  radially  tapered  in  a converging 
sense,  this  fraction  can  become  as  large  as  unity.  If  the  faces  become 
radially  divergent,  this  fraction  can  drop  to  near  zero. 

Over  the  years  much  effort  has  been  made  to  attempt  to  control  this 
radial  taper  to  obtain  an  advantageous  compromise  (a  converging  film) 

* 

There  may  still  be  abrasive  wear  even  if  the  surfaces  do  not  touch. 
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between  fluid  film  load  support  and  leakage.  Davies,  Hooke,  and 
O'Donoghue  developed  a radial  seal  model  based  on  elastic  deflection  [3], 
[4] . Cheng  and  Snapp  [5]  have  shown  the  effect  of  radial  taper  in  face 
seals.  More  recently,  Snapp  and  Sasdelli  [6]  have  experimentally  demon- 
strated the  dramatic  difference  in  friction  and  leakage  caused  by  vari- 
ous radial  tapers.  Metcalf  has  published  extensively  [7] , [8] , [9]  on 
theory  and  methods  to  control  the  radial  taper  in  a seal.  Hie  difficulty 
in  using  a controlled  radial  taper  as  a source  for  lubrication  is  that 
it  is  very  difficult  to  adequately  control  the  taper  because  of  wear, 
thermal  distortion,  and  mechanical  distortion.  As  to  whether  or  not 
this  method  of  providing  relieible  seal  lubrication  will  emerge  as  a 
widely  used  approach  remains  to  be  seen. 

The  second  possible  means  of  developing  fluid  pressure  at  the  inter- 
face is  by  various  hydrodynamic  mechanisms.  The  earliest  demonstration 
that  hydrodynamic  effects  occur  in  face  seals  was  performed  by  Denny 
[10].  Denny's  results  show  that  film  thickness  depends  upon  speed — a 
relationship  that  clearly  indicates  the  presence  of  hydrodynamic  pres- 
sure generation. 

Many  theories  have  been  advanced  to  explain  such  hydrodynamic 
effects  in  face  seals.  A review  paper  by  Nau  [11]  provides  a good 
source  of  comparison  of  the  pressure  generation  mechanisms  proposed  up 
to  the  time  the  paper  was  published.  Nau  states  that  it  is  generally 
accepted  that  the  "majority  of  seals  of  all  kinds  depend  upon  the  process 
of  hydrodynamic  lubrication  for  their  satisfactory  operation."  No  con- 
clusion is  reached  on  what  mechanism  is  predominant,  however.  Mayer  [12] 
also  reports  some  experiments  on  the  measurement  of  hydrodynamic  pressure. 

Several  schools  of  thought  for  hydrodynamic  pressure  generation 
mechanisms  are  evident  in  the  more  recent  literature.  The  first  of  these 
is  that  micropads  or  microasperities  (either  intentional  or  by  virtue  of 
the  materials  themselves)  cause  pressure  generation  by  acting  as  small 
step  bearings.  Hamilton,  Walowit,  and  Allen  [13]  provide  a solution  for 
lubrication  based  on  surface  microirregularities  and  associated  film 
cavities.  The  theoretical  results  agree  qualitatively  with  experiment. 

In  a later  paper  by  Anno,  Walowit,  and  Allen  [14],  tests  were  made  on 
surfaces  with  planned  microasperities  and  the  results  were  compared  to 
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theory.  These  tests  suggest  that  the  use  of  planned  microasperities  is 
an  effective  method  for  lubricating  the  surfaces  of  face  seals.  A later 
paper  by  Anno,  Walowit,  and  Allen  [15]  reported  results  on  protruding  and 
negative  microasperities  as  well  as  the  leakage  from  a microasperity  lu- 
bricated face  seal.  In  a paper  on  generally  the  same  sub ject,  Kojabashian 
and  Richardson  [16]  reported  on  a micropad  model  where  the  statistical 
distribution  of  pads  produced  by  wear  on  a carbon  surface  was  determined. 
Using  this  distribution  of  pads  and  a step-bearing  approximation,  per- 
formance predictions  were  made  which  were  generally  in  agreement  with 
e^qjerience. 

During  the  later  part  of  the  60 's  an  extensive  seals  research  pro- 
gram was  undertaken  by  General  Electric  for  NASA.  Some  of  this  work  is 
reported  in  Reference  [17] , and  there  are  several  other  reports  from  the 
same  study.  This  work  is  representative  of  a second  school  of  thought 
which  is  that  face  seal  hydrodynamics  are  related  to  seal  waviness,  mis- 
alignment, and  eccentricity.  Much  of  the  work  has  been  reported  in  the 
literature  by  Findlay  and  Sneck. 

Findlay's  first  paper  [18]  presents  a mathematical  model  (short 
bearing  approximation)  for  the  pressure  distribution  caused  by  a wavy 
seal  surface  (several  waves  circumferentially  around  the  face)  (see 
figure  1-2)  and  cavitation.  In  a second  paper  [19],  Findlay  presents 
the  results  of  some  experimental  work  carried  out  for  a wavy  (actually 
misaligned)  seal.  This  paper  attempts  to  verify  experimentally  the 
leakage  predicted  by  the  above  cavitation  model. 

Cavitation  of  the  type  predicted  by  Findlay's  theory  has  been  ob- 
served by  several  investigators.  Nau  has  carried  out  extensive  research 
on  cavitation  in  thin  films  [20]  in  general,  and  more  recently  he  has  con- 
ducted a series  of  experiments  directed  toward  seals  using  a glass  plate 
as  one  of  the  seal  faces  [21].  Different  cavity  patterns  could  be  ob- 
served depending  on  the  loading  and  distortion.  Cavitation  has  also  been 
observed  by  Orcutt  [22]  under  much  different  conditions  using  water  as  the 
fluid. 

Sneck  has  published  a number  of  papers  dealing  with  various  hydro- 
mechanical effects  in  seals.  In  his  first  paper  [23],  Sneck  studies  the 
effects  of  misalignment  and  surface  waviness  on  leakage  for  laminar  flow. 


In  a second  paper  [24]  treating  the  same  subject  for  turbulent  flow, 
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Sneck  found  that  the  effects  of  waviness  and  misalignment  are  less. 
Sneck  has  considered  several  additional  related  phenomena  in  other 
papers  [25],  [26],  [27],  [28]. 


More  recently  Stanghan-Batch  [29]  reported  results  of  experiments 
on  hydrodynaunic  pressure  in  a carbon  face  seal  operating  in  oil.  He 
found  that  there  is  a periodic  variation  in  pressure  at  the  face  that  is 
consistent  hydrodyneimically  with  periodic  variations  in  film  thickness. 

He  also  found  that  initially  flat  carbon  rings  wear  in  a manner  such  that 
two  high  spots  are  produced  when  the  ring  is  measured  after  testing  and 
that  there  is  some  correspondence  of  these  waves  with  pressure. 

One  of  the  most  comprehensive  studies  on  seal  lubrication  mechanisms 
has  been  contributed  by  Pape  [30].  Based  on  experiments  and  theory, 

Pape  concludes  that  seal  waviness  is  responsible  for  hydrodynamic  lubri- 
cation in  face  seals . 

In  relation  to  the  waviness  theory,  this  author  has  made  an  extensive 
investigation  of  the  waviness  characteristics  of  face  wear  on  carbon  seal 
rings  [31] , [32] , [33] . It  now  appears  that  most  carbon  rings^  do  not 
wear  uniformly  but  rather  in  some  wavy  pattern,  with  two  peaks  being  most 
common  (see  figure  1-2)  . 

Burton  et  al.  have  published  several  papers  [34]  , [35]  , [36]  , [37]  , 
[38]  concerned  with  the  role  of  thermal  expansion  in  contacting  sliding 
surfaces.  The  relationship  of  thermal  expansion  to  instability  has  been 
excimined  both  for  dry  sliding  and  in  lubricant  films.  Lebeck  [39]  has 
also  investigated  thermoelastic  instability  with  wear  for  the  case  of  the 
ring  geometry  used  in  seals.  These  investigations  are  related  to  the  in- 
vestigation herein  in  that  they  provide  an  understanding  of  how  seal  sur- 
faces distort  under  thermal  loading  and  show  under  what  conditions  wavy 
seal  thermal  instability  might  occur. 

Quite  recently,  Ludwig  [40]  has  summarized  many  of  the  published  seal 
models.  He  has  formulated  a new  set  of  models  based  on  interactions  of 
waviness,  misalignment,  inertia,  and  secondary  seal  friction.  Ludwig  and 

seal  ring  made  entirely  of  carbon,  that  is,  having  no  rigid  metal 
support. 
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Allen  [41]  have  worked  out  the  details  of  a lubrication  model  involving 
secondary  seal  friction  and  misalignment. 

To  determine  which  of  the  hydrodynamic  lubrication  theories  have 
some  validity,  the  experimental  evidence  must  be  examined.  First,  con- 
sidering the  microasperity  load  support  mechanism,  cavitation  streamers 
of  the  type  that  result  from  microasperity  lubrication  have  been  ob- 
served by  Orcutt  [22]  and  researchers  at  Battelle  [42]  by  using  a trcins- 
parent  material  for  one  of  the  seal  rings.  Certainly  these  observations 
suggest  that  microasperity  lubrication  plays  some  role  in  the  lubrication 
of  face  seals,  but  the  actual  extent  or  fraction  of  the  load  support  is 
not  known. 

Several  researchers  [22],  [29],  [30],  [42]  have  determined  the 
existence  of  circumferential  pressure  and/or  film  thickness  variations 
which  in  turn  indicates  that  hydrodynamic  lubrication  and  load  support 
do  occur  due  to  waviness  of  the  seal  rings.  Again,  the  fraction  of  load 
support  developed  by  this  mechanism  in  face  seals  in  actual  service  is 
not  known. 
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PROPOSED  MODEL- WA VI NESS  CAUSED  HYDRODYNAMIC  LUBRICATION  WITH  ELASTIC 
DEFORMATION,  ASPERITY  CONTACT,  AND  WEAR 

There  is  significant  experimental  evidence  for  either  of  the  above 
mechanisms.  However,  in  this  research  program,  the  waviness  theory  has 
been  chosen  for  extensive  further  study.  The  reason  for  this  choice,  in 
addition  to  the  experimental  support  already  cited,  is  that  waviness  is  i 

a controllable  parameter  in  seal  design.  A better  understanding  of  how 
waviness  produces  hydrodynamic  load  support  in  face  seals  could  lead  i 

relatively  quickly  to  improved  designs.  However,  as  discussed  in  chap-  1 

ter  7,  the  phenomenon  of  microasperity  lubrication  must  actually  be  con- 
sidered in  conjunction  with  waviness.  So  the  above  choice  is  not  as 
restrictive  as  it  might  at  first  appear. 

In  much  of  the  previous  research  on  waviness  certain  idealizations 
have  been  made  to  obtain  solutions.  What  is  needed  by  way  of  further  j 

research  into  this  area  is  to  include  the  effects  of  wear,  surface  j 

roughness,  and  elastic  deformation.  These  factors  become  particularly  ] 

important  for  the  class  of  heavily  loaded  seals  of  greatest  interest. 
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Much  of  the  previous  experimental  work  on  waviness  has  been  for  lightly 
loaded  seals  where  the  more  simple  models  retain  some  of  their  validity. 


In  the  heavily  loaded  seals  of  interest  here,  wear  of  the  seal  faces 
acts  as  the  primary  limitation  for  better  seal  performance.  The  term 
wear  is  used  broadly  to  include  a gradual  wearing  away  and  deterioration 
as  well  as  heat  checking  and  blistering  of  the  faces.  All  of  these  types 
can  be  attributed  in  part  to  inadequate  lubrication  or  inadequate  fluid 
film  pressure  load  support.  Hydrodynamic  lubrication  due  to  waviness, 
as  the  results  will  show,  may  or  may  not  play  a part  in  the  lubrication 
of  such  seals.  However,  in  order  to  establish  the  role  of  this  mechan- 
ism, the  factors  mentioned  must  certainly  be  included  in  the  model. 

Now,  examining  the  various  elements  of  the  proposed  model,  the  pres- 
ence of  asperity  to  asperity  contact  in  seals  can  be  justified  in  two 
ways.  First,  some  asperity  to  asperity  contact  must  exist  to  account  for 
the  wear  found  in  seals  operating  in  non-abrasive  environments.  Secondly, 
it  can  be  shown,  particularly  for  seals  sealing  low  viscosity  fluids  or 
for  heavily  loaded  seals,  that  the  film  thickness  required  for  equili- 
brium using  hydrodynamic  pressure  only  is  of  the  same  order  of  magnitude 
as  the  surface  roughness.  In  some  cases  the  film  thickness  is  less  than 
the  surface  roughness.  Thus  some  asperity  to  asperity  contact  would  be 
expected  and  is  required  to  support  that  fraction  of  the  load  not  sup- 
ported by  hydrodynamic  pressure. 

Models  of  asperity  load  support  in  conjunction  with  hydrodynamic  sup- 
port have  been  developed  for  bearings  by  Christensen  [43]  , [44]  and 
Thompson  and  Bocchi  [45] . These  models  predict  a friction  and  wear  be- 
havior for  thin  film  thickness  bearings  that  is  consistent  with  experi- 
mental results.  These  same  models  have  been  extended  herein  to  study 
asperity  contact,  hydrodynamic  load  sharing,  and  wear  in  mechanical 
face  seals. 

For  small  film  thicknesses,  the  effects  of  surface  roughness  on 
hydrodynamic  lubrication  must  be  considered.  Christensen  [44]  has  de- 
veloped a statistical  method  that  can  be  used  to  calculate  hydrodynamic 
load  support  where  the  surface  roughness  effects  are  included  in  the  film 
thickness  variation.  Thus,  models  are  available  for  the  purposes  of 
accounting  for  these  important  effects. 
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In  summary,  the  model  developed  is  based  upon  a complex  Interaction 
among  the  following  phenomena: 


1) 

Initial  waviness 

2) 

Elastic  deflection  of  waviness  due  to  asperity 

and 

fluid 

pressures 

3) 

Wear 

4) 

Hydrodynamic  lubrication  due  to  waviness  including 

the  ef- 

fects  of  rough  surfaces  eind  cavitation 

5) 

Asperity  to  asperity  contact,  load  support  and 

friction 

6) 

Hydrostatic  effects. 

The  interactions  of  these  phenomena  are  all  time  dependent  because  of 
wear. 

SOURCE  OF  WAVINESS 

Waviness  of  the  seal  faces  is  essential  to  the  hydrodynamic  lubri- 
cation theory  developed.  In  the  discussion  of  previous  research,  it  was 
pointed  out  that  experimental  evidence  shows  the  existence  of  such  wavi- 
ness, emd  there  are  several  possible  mechemisms  that  can  act  to  produce 
waviness  in  mechanical  face  seals. 

This  author  [31]  has  made  extensive  measurements  of  the  waviness  of 
initially  flat  carbon  seals  after  some  period  of  operation.  The  conclu- 
sion reached  is  that  a large  portion  of  the  waviness  measured  (on  the 
order  of  100 's  of  y-in  in  many  cases)  is  due  to  the  fact  that  the  seal 
ring  is  distorted  into  a wavy  configuration  due  to  driving  forces  [32], 
[33].  After  a period  of  operation,  the  face  wears  relatively  flat — while 
the  ring  is  distorted.  Therefore,  after  the  ring  is  unloaded  and  returns 
to  its  undistorted  shape,  it  appears  to  be  wavy.  In  addition,  Pape  [30] 
has  found  that  a significant  waviness  generally  remains  after  seal  faces 
are  lapped  flat.  Lebeck  [31]  has  also  measured  significant  waviness  in 
as  lapped  seals.  There  are  other  sources  of  waviness  as  well.  Any 
change  in  operating  conditions,  such  as  temperature  or  pressure,  can  lead 
to  changes  in  waviness.  Also  localized  heating  due  to  waviness  itself  can 
change  the  amount  of  waviness.  Finally,  seals  with  a planned  amount  of 
waviness  can  be  designed  so  emy  amount  of  waviness  can  be  introduced 
into  a seal  both  initially  and  over  a longer  term  (see  chapter  6) . 
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From  these  various  sources  there  is  clearly  enough  waviness  in 
face  seals  to  create  temgentially  converging-diverging  regions  in  the 
film  thickness.  What  remains  as  the  question  to  be  answered  herein  is 
the  precise  nature  and  limits  of  such  lubrication  and  how  much  waviness 
is  required  to  cause  significant  hydrodynamic  effects.  Also,  can  wavi- 
ness be  designed  into  a seal  and  provide  improved  performance.  These 
euid  many  other  questions  are  explored  using  the  model  developed  herein. 


CHAPTER  2 


f 


TWO  DIMENSIONAL  EQUILIBRIUM  MODEL 

Based  on  the  previous  discussion,  the  hydrodynamic  problem  to  be 
solved  is  shown  in  figure  2-1.  The  seal  is  subjected  to  outside  and 
inside  pressure.  The  wavy  rough  face  (shown)  is  fixed,  while  the  smooth 
face  rotates.  Given  two  waves,  for  example,  there  will  be  two  points  of 
minimum  film  thickness  between  the  faces  as  shown.  The  seal  may  actually 
contact  in  regions  near  these  points  as  shown.  Angular  rotation  of  the 
mating  face  causes  hydrodynamic  pressure  to  build  up  in  the  converging 
regions  as  the  point  of  minimum  film  thickness  is  approached.  In  the 
diverging  regions,  cavitation  will  occur  as  the  pressure  attempts  to  be- 
come negative.  This  cavitation  is  essential  for  a net  lifting  force,  and 
for  an  outside  pressurized  seal,  the  cavity  will  appear  as  shown  in  fig- 
ure 2-1.  The  fluid  will  flow  in  streamers  across  the  cavity,  and  a full 
film  will  again  develop  at  the  downstream  cavity  boundary. 

The  hydrodynamic  pressure  distribution  for  this  problem  has  been 
solved  by  various  methods  by  Findlay  [18] , Pape  [30] , and  Stanghan-Batch 
and  Iny  [46]  for  perfc jtly  smooth  faces.  These  results  show  that  even  a 
small  waviness  produces  sufficient  load  support  for  complete  liftoff. 
However,  as  discussed  in  chapter  1,  it  is  known  that  many  more  heavily 
loaded  or  low  viscosity  seals  do  not  operate  with  full  hydrodynamic  load 
support  and  complete  separation  because  a definite  wear  results.  Sur- 
faces of  such  seals  contact  during  operation,  and  hydrodynamic  pressure 
is  affected  by  interactions  with  surface  roughness.  Seals  of  this  type 
operate  in  a mixed  friction  regime.  Hydrodynamic  load  support  due  to 
waviness  continues  to  provide  a significant  fraction  of  the  load  support, 
but  asperity  contact  must  provide  the  balance. 

The  problem  to  be  solved  is  to  determine  the  fluid  and  mechanical 
pressure  distributions  for  the  given  configuration  considering  the  effect 
of  roughness.  Given  the  pressure  distribution,  then  load  support  and 
leakage  can  be  calculated  and  the  effect  of  various  parameters  can  be 
studied.  In  the  following  sections  the  various  parts  of  the  model  will 
be  developed  in  detail. 
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Rotation  of  Mating  Face 


h = h^+  hflCos  ne 


Figure  2-1.  Seal  Geometry  and  Interference  Regions 
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MODEL  DEVELOPMENT 
Reynolds  Equation 

Even  though  no  greater  error  is  introduced  for  most  seals  by  using 
rectcingular  coordinates , it  was  decided  to  formulate  the  problem  in 
polar  coordinates  in  case  a wide  face  seal  became  of  interest.  The 
Reynolds  equation  in  polar  coordinates  is  given  by 


Flow  equations  are 


_ -h  3p  ra)h 

12rn  30  2 


= dll  iE 

12n  3r 


where  coordinates  are  defined  in  figure  2-1. 


(2-1) 


(2-2) 


(2-3) 


Surface  Roughness 

Considering  now  the  effects  of  surface  roughness,  several  investi- 
gators [43],  [44],  [47],  [48],  [49],  [50],  [51]  have  developed  models  of 
hydrodynamic  lubrication  which  account  for  the  interaction  of  surface 
roughness  on  hydrodynamic  pressure.  The  work  of  Christensen  and  Tonder 
[43] , [44] , [47]  , [48] , [49]  has  been  selected  as  the  basis  of  the  model 
herein.  The  Christensen  model  has  been  previously  applied  to  the  mixed 
friction  case  where  some  roughness  interference  occurs  [44] . 

The  Christensen  model  for  mixed  friction  [44]  is  based  on  longi- 
tudinal roughness  [See  figure  2-2]  as  opposed  to  transverse  roughness. 

At  first  glance  this  appears  to  be  appropriate  for  mechanical  face  seals 
because  significant  longitudinal  roughness  patterns  are  observed  on  the 
faces  of  contacting  seals.  However,  it  is  known  that  some  transverse 
roughness  does  exist  in  face  seals,  and  the  assumption  of  zero  transverse 
roughness  does  introduce  some  error  into  the  model  because  theoretically 
there  will  be  zero  side  leakage  in  regions  of  roughness  interference. 


As  the  results  will  later  show,  the  use  of  the  longitudinal  rough- 
ness assumption  does  cause  some  unexpected  behavior  in  seal  performauice , 
and  results  using  this  approach  are  later  compared  to  results  obtained 
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Figure  2-2.  Longitudinal  Roughness  and  Geometry 
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using  an  isotropic  roughness  assumption.  Isotropic  roughness  (see 
figure  2-3)  represents  a mathematically  ideal  case  in  which  flow  re- 
sistance and  roughness  statistics  are  equal  in  all  directions.  The 
reality  for  seal  faces  is  that  roughness  is  anisotropic,  possibly  con- 
sisting of  longitudinal  valleys  and  ridges  of  short  length,  thus  giving 
rise  to  side  leakage  in  regions  of  contact.  However,  models  for  aniso- 
tropic roughness  were  not  available  at  the  time  when  this  work  was 
performed.  The  isotropic  model  also  proposed  by  Christensen  [43]  emd 
Tender  [49]  has  been  used  as  a first  approximation  to  actual  seal  rough- 
ness. 

As  to  the  mathematical  correctness  of  the  longitudinal  and  iso- 
tropic models  of  Christensen  and  Tender,  it  appears  that  the  longitudinal 
model  is  reasonably  valid.  However,  there  is  some  question  about  the 
accuracy  of  the  isotropic  model.  Again  because  no  other  models  were 
available,  the  isotropic  results  herein  are  based  on  the  Christensen 
and  Tonder  model  and  must  be  viewed  with  the  possibility  that  some  cor- 
rection will  be  required  as  improved  models  become  available.  These 
questions  are  addressed  in  chapter  7. 

Now  proceeding  with  the  development  of  the  model,  according  to 
Christensen  [44] , the  total  film  thickness  for  rough  surfaces  cem  be  de- 
scribed by  a function  of  the  form 

H = h(r,0)  + hg(  C)  (2-4) 

where  h(r,6)  is  the  nominal  film  thickness  which  varies  due  to  waviness 
and  radial  taper  and  h^  is  the  random  part  of  film  thickness  due  to  sur- 
face roughness.  A polynomial  approximation  of  the  Gaussian  density  dis- 
tribution function  is  used  to  describe  the  surface  roughness 

f(h  ) = — ■ (c^  - hf)  , - c < h < c (2-5) 

c corresponds  to  three  standard  deviations  and  represents  the  peak 
roughness  eunplitude.  It  is  assumed  that  the  roughness  distribution  re- 
mains constant  with  time  and  coordinates  r and  9.  It  is  also  assumed  that 
the  roughness  occurs  on  one  surface  only  as  shown  in  figures  2-2  and  2-3. 
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Figure  2-3.  Isotropic  Roughness  and  Geometry 


Longitudinal  Roughness 


According  to  Christensen  [44] , longitudinal  roughness  effects  (see 
figure  2-2)  can  modify  the  Reynolds  equation.  In  dimensionless  form, 
the  equations  for  the  two  dimensional  flow  region  (ii  > 1)  become: 


h > c 


where 


(2-6) 


(2-7) 

(2-8) 


P = 


2 

r u)n 
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‘’r 

q = 

r r ton 
o 

- ^ ^9 

r liic 
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(2-9) 

(2-10) 

(2-11) 

(2-12) 

(2-13) 


0)  = angular  velocity 

and  where  E is  the  expectancy  operator  which  is  equal  to 


c 

E(  ) = / ( )f(hg)dhg  (2-14) 

a 

where 


a = -c 

where 

h > c 

(2-15) 

a = -h 

where 

h < c 

(2-16) 
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When  h < 1,  the  roughness  interference  region,  radial  flow  becomes 
zero  because  the  longitudinal  roughness  acts  as  a series  of  dams  against 
the  radial  flow.  In  this  case  tangential  flow  is  given  by 


•(- 


E(H  ) 3p  . rE(H)  \ r 


(2-17) 


where  b.  is  the  fraction  of  the  area  available  to  flow.  The  b,  term  re- 
n n 

presents  a correction  to  Christensen's  original  work  [44]  and  reflects 
the  condition  of  reduced  flow  area. 


The  Reynolds  equation  for  h < 1 is  obtained  from  the  continuity 


condition: 


(2-18) 


So,  from  equation  (2-17) 


- ff- 

de  \ 


E(H  ) dp  . rE(H)  \r 


)\  -0 


(2-19) 


Isotropic  Roughness 

The  Reynolds  and  flow  equations  as  modified  by  Christensen  [43]  and 
Tonder  [49]  for  isotropic  roughness  are: 


j w(lf 


(E(H)bj^) 


- ^ /_  EiHJ.|£^ 


- ae 

12r 


\ ‘2  a:)  " 


(2-20) 


(2-21) 


(2-22) 


For  the  given  roughness  distribution  function,  equation  (2-5)  , the 
needed  expectancy  functions  for  the  longitudinal  and  isotropic  cases  as 
determined  from  equatioi.  (2-14)  are 


For  h > c 


E(H)  = h 

E(H^)  = + h/3 


(5h^-l)ln 


26h  + 


30h^ 


= 1.0 
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For  h £ c 

E(H)  = 


E(H^ 
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(2-23) 

(2-24) 

(2-25) 

(2-26) 

(2-27) 


(2-28) 

(2-29) 


Boundary  Conditions 

Witli  reference  to  figure  2-1,  at  tlie  inside  and  outside  boundaries 
the  boundary  conditions  are 


p = p 

at  r = r 

(2-30) 

o 

o 

p = Pi, 

at  r - r. 

1 

(2-31) 

The  type  of  cavitation  expected  in  face  seals  is  similar  to  that  en- 
countered in  journal  bearing  where  flow  separates  into  streamers  at  some 
upstream  condition  where  pressure  would  otherwise  become  negative.  A 
full  film  is  reestablished  at  some  point 'downstream  where  continuity  is 
satisfied  (see  figure  2-1) . The  cavity  boundary  conditions  used  herein 
are  those  given  by  Ja)cobsson  and  Floberg  (52l.  For  a non-asperity  smooth 
surface  contact  condition,  at  the  upstream  boundary,  the  cavity  boundary 
conditions  are 


P 


= P, 


(2-32) 


(2-33) 
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where  n is  the  direction  normal  to  the  cavity  boundary  as  shown  in  fig- 

I 

ure  2-1,  Now,  for  the  case  of  interest  here  where  the  cavity  boundary  , 

may  pass  through  a region  of  roughness  interference,  the  upstream  bound- 
ary conditions  remain  the  same.  The  second  downstream  boundary  condition 
must  be  re-derived  based  on  the  flow  equations  (2-7) , (2-8)  and  (2-17)  for 
longitudinal  roughness  and  equations (2-21)  and  (2-22)  for  isotropic  rough- 


ness and  continuity.  The  result  for  longitudinal  roughness  is 


(2-36) 


where  the  subscripts  and  refer  to  the  upstream  and  downstream 

boundary  conditions  as  shown  in  figure  2-1.  The  functions  bj^  are  as  de- 
fined previously  for  h ^ c and  h > c.  For  the  special  case  where  the 
downstream  cavity  boundary  falls  in  the  roughness  interference  region, 
then  the  above  equation  is  modified  accordingly. 

For  isotropic  roughness  the  result  is 
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^ E(H^(2)  = 6?  cosaJE(fi)^  b 

3n  Q L.  Qy 


• E(hI.  b. 


(2-37) 


In  the  isotropic  roughness  case  the  special  case  where  -►  -ir  when  the 
downstream  cavity  boundary  is  in  the  roughness  interference  region  does 
not  apply. 
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Roughness  Interference  Region 


In  regions  where  there  is  roughness  interference,  additional  load 
support  is  obtained  by  the  flattening  of  the  peaks  of  the  asperities. 
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It  is  assumed  that  the  asperity  pressure  p^  is  equal  to  the  compressive 
strength  of  the  softer  material  [44] , cind  that  this  pressure  is  constant. 
Thus,  the  average  pressure  due  to  asperity  contact  across  the  seal  face 
is  equal  to  the  fraction  of  the  interference  area  times  the  asperity 
pressure  or 


P„  = P„(l  - b,  ) 
mm  n 

a 


(2-38) 


For  this  relationship  to  be  valid,  it  is  assumed  that  surface  roughness 
is  homogeneous  and  that  after  the  faces  are  pressed  together  in  regions 
of  contact,  the  roughness  distribution  is  simply  truncated  giving  rise 
to  asperity  load  support.  It  is  further  assumed  that  this  truncated 
distribution  does  not  change  with  time  (except  by  wear) . The  validity 
of  this  assumption  depends  greatly  on  the  wear  model  assumed.  This 
topic  is  discussed  in  chapter  7. 

Load  Support,  Leakage,  and  Friction 

Total  load  support  is  determined  by  the  hydrodynamic  and  asperity 
pressure  distributions 


where 


2ir 


W=  / / [p  + p^(l-b^)]?d?d0 


0 r, 


(2-39) 


W = 


Wc 
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(2-40) 


The  load  support  required  for  equilibrium  is  determined  by  the 


sealed  pressures  p^  and  p^  and  the  balance  ratio  B.  Thus 
W*  = TT(1  - r^)  [p^B  + p^(l  - B)  ] 


(2-41) 


assuming  spring  pressure  is  negligible. 

Leakage  for  the  longitudinal  roughness  case  is  given  by  integration 
of  the  flow  equation  (2-8) 


(2-42) 
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where 
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where  the  lea]cage  function  is  zero  for  the  edges  adjacent  to  the  cavity 
and  in  the  roughness  interference  region. 

For  the  isotropic  roughness  case,  leakage  is  given  by 


5i  = 


2t!  -3 

, _ E(HJ.  3£  ae 
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In  the  full  film  region,  tangential  fluid  friction  shear  stress  is 
given  by 


(2-47) 


Assuming  a friction  radius  r^  where 
- ■'i’ 


2 2. 

then  the  total  fluid  friction  force  in  the  full  film  region  is 

-2f 


(2-48) 
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In  the  cavitated  film  region,  the  above  expression  becomes 


F = // 


cavitated 


r El  — 


b^drde 
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(2-50) 


‘ f area 


. (i> 

where  the  ratio  — accounts  for  the  fraction  of  full  film. 


The  average  friction  shear  stress  due  to  asperity  contact  is  given 
by  an  expression  similar  to  the  asperity  load  support. 


T = p^(l  - bj^) 


(2-51) 


where  p^  is  the  shear  strength  of  the  asperities.  The  total  friction 
force  due  to  asperity  contact  is 
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a c - , 

r^  roughness 

interference 
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(2-52) 


also, 


(2-53) 


The  coefficient  of  friction  is 


F,  + F + F F,.  + F + F 
f c a f c a c 


(2-54) 


Two  additional  expectancy  functions  are  needed  to  evaluate  equations 
(2-49)  and  (2-50) . These  are 


_/l  V 35  r -2.  h+1  11  ^ 11  - ^ 23  -2 


a -3  27  r4  ^ -5  ^ 147  -6  - 

-3h-^h+h+—  h ,h<l 


(2-55) 


(1-h^)  + II  h - h^  + 2h^  , h > 1 (2-56) 


The  limits  on  the  expectancy  integral,  equation  (2-14) , when  de- 
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riving  equation  (2-55)  are  from  -h  + A to  c where  A is  a small  non  zero 
positive  number  to  limit  the  value  of  E(l/H).  Otherwise  viscous  fric- 
tion goes  to  infinity  as  the  film  thickness  becomes  zero.  A realistic 
value  for  A can  be  obtained  by  setting  the  maximum  viscous  shear  stress 
equal  to  the  shear  strength  p^.  This  gives 


A = 

r 

Ps  ° 


(2-57) 


Film  Thickness  and  Deflection 

In  general,  nominal  film  thickness  is  a function  of  0 and  r.  For 
the  cases  here,  it  is  assumed  that  the  faces  are  radially  parallel.  The 
radial  variation  case  is  treated  later.  It  is  assumed  that  the  tan- 
gential film  thickness  variation  is  of  an  arbitrary  character  within 
each  of  n equal  periods  arovind  the  seal  face.  The  nominal  film  thickness 
is  represented  by  a Fourier  series  as  follows. 


h = h 
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(2-58) 


Since  n can  take  the  value  of  one,  this  film  thickness  description  is 
completely  general. 

The  terms  h^  represent  the  initial  waviness  and  the  modification  of 
initial  waviness  due  to  wear.  The  terms  h^  represent  the  contribution 
to  waviness  caused  by  elastic  deflection. 

The  amplitude  of  the  nth  harmonic  of  the  waviness  is  given  by 

il/2 
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(2-59) 


The  load  support  as  a function  of  6 is  obtained  as  follows; 


p(0) 
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(2-60) 
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(2-61) 


when  p(0)  is  expressed  as  load  per  unit  angle.  This  function  can  be 
represented  as  a Fourier  series 


p(0)  = p + I p cos  nj9  + I P.  sin  nje 


D=1  D 


1=1  3 


<2-62) 


where  p. , p and  p.  are  the  Fourier  coefficients.  The  deflection  for 
o a . b . 

a seal  ring  3an  then^be  obtained  from  reference  [33].  For  a face  pres- 
sure distribution  represented  by  a Fourier  series,  reference  [33]  shows 
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where 
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(2-65) 
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J is  the  moment  of  inertia  about  a radial  axis  through  the  cen- 

X 

troid  of  the  ring,  E is  the  modulus  of  elasticity  for  the  ring,  and  r^ 
is  the  radius  to  the  centroid.  For  the  case  under  study,  it  is  assumed 
that  the  stiffness  of  the  mating  ring  is  infinite  so  that  only  the 
stiffness  of  the  seal  ring  is  required.  However,  in  general 


(2-67) 


primary 

ring 


In  many  of  the  solutions  to  follow,  only  the  combined  or  net 

waviness  h is  used  as  a parameter,  since  this  parameter  alone  is 
n 

sufficient  to  describe  the  total  effect  of  waviness.  Later  in  chapter  5, 
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equations  (2-63)  and  (2-64)  are  used  to  predict  the  initial  waviness 
required  to  obtain  a given  net  waviness. 

METHOD  OF  SOLUTION 

The  problem  to  be  solved  is  to  find  the  complete  pressure  distribu- 
tion for  the  configuration  and  boundary  conditions  shown  in  figure  2-1 
for  an  arbitrary  film  thickness  h.  Then  load  support,  leakage,  euid  fric- 
tion can  be  calculated  using  previously  given  equations,  euid  the  mean 
film  thickness  h^  can  be  adjusted  so  that  the  applied  load  equals  the 
supporting  load. 

The  most  difficult  aspect  of  this  problem  is  that  the  cavity  shape 
(figure  2-1)  is  part  of  the  solution  and  not  known  beforehand.  In  fact, 
the  boundary  conditions,  equations  (2-32) , (2-33)  and  (2-34)  and  (2-36) 
or  (2-37)  require  the  boundary  to  assume  a particular  shape.  In  spite 
of  this  difficulty,  the  solution  has  been  obtained  by  both  Findlay  [18] 
and  Pape  [30] , except  that  surface  roughness  effects  were  not  included. 
However,  both  of  these  investigators  indicated  that  the  numerical  meth- 
ods they  used  were  somewhat  time  consuming  because  of  the  difficulty  in 
finding  the  correct  cavity  boundary  shape.  In  fact  Pape  [30]  used  a 
method  of  trial  and  error  by  hand  calculation  to  obtain  the  correct 
cavity  shape. 

Because  of  a need  in  the  present  research  to  be  able  to  converge 
on  pressure  distribution  and  cavity  shape  in  a short  amount  of  computer 
time  so  that  equilibrium  solutions  could  be  obtained,  a different  approach 
to  the  solution  of  the  outlined  problem  was  developed.  This  approach 
includes  the  effects  of  surface  roughness  as  required  for  the  solution 
of  interest,  but  has  application  to  any  type  of  two  dimensional  unde- 
fined boundary  problem  such  as  solved  by  Jakobsson  and  Floberg  [52]  for 
journal  bearings  using  more  time  consuming  methods. 

The  method  developed  works  as  follows:  The  p > p^  pressure  region 
is  solved  using  normal  successive  overrelaxation  methods.  The  upstream 
boundary  conditions  (figure  2-4)  , equations  (2-32)  and  (2-33)  are  satis- 
fied by  setting  p = p^  when  otherwise  p would  be  less  than  p^.  The 
downstream  boundary  is  established  by  moving  the  right  hand  boundary 


29 


ROUGHNESS  r-UPSTREAM  BOUNDARY  ^ DOWNSTREAM 


I 


3:  o:  o ( 

(S  LlJ  I— ■ 

=3 1-  ta 
O Z LU 

a:  M Of 


• ® ® 
= 00 
o /0i0  0 0 

/o  3 

0 0 0 0 
0_0  0 0 
•\0i@  0 0 

HH 

;•  ^000 


I- 


o z 

o. 

LlJ  O 

QC  0^0. 
=D  LU 

00  LU  Z 
00  Ul 

UJ  =>  O Z 
00  O 

0-00  O 
UJ  Z O 

z q:  S 

^ a.  uj  z 

•-I  Qc:  < 

u-  >-  >—  UJ 
»—  00  a 
—i»-*z»— 
— I > 3 00 
=)  < Q Ci- 
u.  o o 2 


0S^<^ 


30 


CAVITY  BOUNDARY  POINT 


i 

I 

i 


I 

i 

i 

i 


\ 

i 

\ 

I 

! 

f 

t 

\ 

I 

i 

j 

I 

I 

I 

i 

I 


I 


nodes  to  the  left  or  to  the  right  so  that  the  pressure  at  the  boundary  I 

is  reduced  to  zero  and  the  normal  derivative  condition  equation  (2-36)  I 

or  (2-37)  is  also  satisfied.  | 

Now  considering  the  solution  in  greater  detail,  figure  2-4  shows  | 

the  basic  problem  to  be  solved  in  discrete  form.  The  distance  repre-  j 

sented  by  I = 1 to  I = IMAX  represents  one  of  the  n periods  around  the  j 

seal.  It  is  assumed  that  only  one  cavity  occurs  within  this  period.  \ 

The  cavity  is  defined  as  shown  by  connecting  a series  of  comer  points 
lUP(J)  2ind  IDWN(J)  which  are  the  last  points  where  p > p^  before  enter- 
ing the  cavity.  It  is  also  assumed  that  p^  = p^  so  that  the  cavity  is 
open  to  the  inside  of  the  seal. 

For  the  solution,  the  numerical  formulation  of  the  Reynolds  equa- 
tion, equations  (2-6  ) and  (2-19)  for  longitudinal  roughness  and  equa- 
tions (2-20)  for  isotropic  roughness,  was  obtained  using  central  dif- 

2 

ference,  (p^^  - p_j^)/2AX  for  the  first  derivative  and  (pj^  - 2p^  + p_j^)/AX  i 

for  the  second  derivative.  Equations  (2-6)  and  (2-19)  or  (2-20)  can  i 

thus  be  reduced  to  the  following  numerical  form. 

P(I,J)  = AIMl  * P(I-1,J)  + AIPl  * P(I+1,J) 

+ AJMl  * P(I,J-1)  + AJPl  * P(I,J+1)  + B(I,J)  (2-68) 

where  the  coefficients  are  functions  of  the  expectancy  equations,  mesh 
size,  and  mean  film  thic)cness  at  each  point.  The  relationships  for  the 
coefficients  are  given  in  Appendix  A.  For  the  longitudinal  roughness 
case,  there  will  be  a shift  between  equations (2-6)  and  (2-19)  when  moving 
from  point  0 to  point  0 in  figure  2-4.  To  ensure  continuity,  an 
extra  column  of  points  was  added  through  point  O at  the  edge  of  the 
interference  region  and  special  equations*  were  derived  for  points 
(D,  0 , and  (3)  . 

Angle  Definition 

Before  the  iterative  method  can  be  described,  the  angles  shown  in 


The  necessity  for  doing  this  was  discovered  later  when  it  was  found  that 
load  support  as  a function  of  h^  was  discontinuous. 
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figure  2-4  must  be  defined.  Angle  definition  for  the  downstream  corner 
points  is  obtained  as  shown  in  figure  2-5.  Because  of  the  angle  defini- 
tion used,  a movement  of  IDWN(J)  to  the  left  or  right  does  not  cause  the 
emgle  for  the  given  J to  change.  This  helps  to  stabilize  the  solution. 
From  the  geometry  of  figure  2-5,  an  angle  is  assigned  to  each  of  the 
downstreeun  corner  points  IDWN(J) . 

Figures  2-6  and  2-7  show  how  eui  angle  is  assigned  to  the  first  and 
last  corner  points.  The  angle  for  the  top  cavity  point  (ITOP,  figure  2-4) 
is  assigned  the  value  of  a = —n/2  so  that  the  normal  direction  is  ra- 
dially inward. 

Given  an  angle  definition  for  each  of  the  corner  points  and  the  top 
cavity  point,  angle  values  for  all  other  downstream  cavity  points  ^ 

(see  figure  2-4)  are  found  by  interpolation  between  the  adjacent  corner 
points . 

It  was  found  in  the  solution  to  the  problem  that  with  each  iteration 
the  corner  point  angles  would  change  significantly,  and  this  created  some 
instability  in  the  solution.  To  minimize  this  problem,  corner  point 
angles  were  always  taken  as  the  average  of  the  values  at  four  successive 
iterations.  As  the  solution  approaches  a stable  cavity  position,  the 
error  caused  by  this  averaging  goes  to  zero. 

An  angle  definition  is  not  actually  needed  for  the  upstream  boundary. 

Selection  of  JTOP  euid  ITOP 

JTOP  (figure  2-4)  was  selected  as  the  bottommost  row  of  non  cavity 
pressure.  For  cavity  angle  definition,  it  was  assumed  that  the  cavity 
boundary  extends  to  this  row  as  shown  in  figure  2-4. 

In  the  isotropic  case  ITOP  was  selected  as  being  midway  between 
lUP(JTOP-l)  and  IDWN(JTOP-l)  as  shown  in  figure  2-6.  This  procedure 
was  also  followed  for  longitudinal  roughness  where  contact  does  not 
occur.  For  the  longitudinal  roughness  case  where  contact  does  occur, 

ITOP  was  selected  as  point^on  figure  2-4.  JTOP  becomes  equal  to 
JMAX-1.  It  was  found  for  contacting  longitudinal  roughness  cases,  point 
(^in  figure  2-4  would  always  assume  cavity  pressure  unless  restricted 
by  the  above  choice.  Numerically  speaking,  a large  and  erroneous  inflow 
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IDWN(J+1)  IDWN(J-l) 


Figure  2-5.  Angle  Definition  - General  C?se 
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would  then  occur  at  point^)L  Based  on  numerical  experimentation,  it 
appears  that  the  cavity  boundary  for  the  longitudinal  case  approaches 
the  outer  seal  boundary  very  closely.  This  situation  is  difficult  to 
simulate  numerically,  so  the  best  approximation  available  was  used — 
that  of  requiring-non  cavity  pressure  at  node @ in  figure  2-4.  Al- 
though this  approximation  allowed  the  inside  and  outside  flows  to  bal- 
ance, some  results  show  unusual  pressure  values  at  this  point. 

This  particular  problem  occurs  only  because  of  the  zero  side 
leakage  assumption  in  the  longitudinal  roughness  case.  Given  even  a 
small  side  leakage,  this  discontinuity  would  not  occur. 

Iterative  Method 

The  sweeping  pattern  used  for  successive  over  relaxation  starts  at 
the  top  left  corner  of  figure  2-4.  The  sweep  is  performed  row  by  row 
so  that  the  cavity  is  swept  in  the  direction  of  motion  as  opposed  to 
being  swept  radially.  Each  complete  sweep  of  all  rows  is  called  one 
successive  over  relaxation  iteration. 

Before  the  iterative  method  can  be  initiated,  some  initial  guess 
must  be  made  for  the  cavity  boundary.  This  is  accomplished  by  using 
equation  (2-68)  in  a successive  over  relaxation  for  a few  iterations 
completely  disregarding  the  cavity  boundary  conditions.  Then  all  points 
where  pressure  is  less  than  p^  are  assigned  to  the  cavity  and  the  corner 
nodes  are  assigned  accordingly. 

After  an  initial  guess  is  obtained,  the  row  by  row  relaxation  is 

continued,  but  the  following  operations  are  performed  as  each  row  is 
swept . 

1)  Start  at  I = 1 and  solve  for  P(I,J)  using  successive  over 
relaxation  and  the  basic  equation  (2-68) . 

2)  As  I increases,  at  some  point  the  pressure  will  become  less 
than  the  cavity  pressure.  Set  all  such  pressures  equal  to  the  cavity 
pressure.  Call  the  node  just  to  the  left  of  the  first  pressure  set  to 
the  cavity  pressure  lUP(J).  Setting  the  pressure  to  the  cavity  pressure 
satisfies  the  upstrecun  boundary  condition  of  p = p^,8p/8n  = 0,  equations 
(2-32)  and  (2-33). 
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3)  From  (lUP(J)  + 1)  to  (IDWN{J)  - 1)  the  pressures  are  set  to  the 


cavity  pressure  p^. 

4)  When  IDWN(J) (point(4)in  figure  2-4)  is  reached,  special  equa- 
tions are  used  to  solve  for  the  pressure  at  IDWN{J)  and  those  boundary 
points  just  to  the  right  of  IDWN(J)  which  are  identified  with  triangles 
in  figure  2-4  (points and(^ , Based  on  the  downstream  flow  boundary 
condition,  equation  (2-36)  or  (2-37) , depending  on  the  type  of  rough- 
ness, 3p/3n  can  be  determined  for  IDWN(J)  and  the  points  to  the  right 
by  noting  E (H)^^ corresponds  to  E(H)  for  6(IUP,J).  The  angle  is  de- 
termined as  mentioned  in  the  previous  sections  for  the  corner  point  and 
for  the  remaining  adjacent  cavity  points  to  the  right.  Given  3p/3n,  then 

= r ^ cos  a (2-69) 

3n 

sin  a (2-70) 

3r  3n 


These  derivatives  are  approximated  for  the  corner  and  cavity  points  by 
first  differences.  For  the  corner  point, 

lE  = p(I,J)  - p(I-l,J)  ,,  7 

30  A0 


3p  ^ P(I,J) 
3r 


AR 


P(I,J-1) 


(2-72) 


These  equations  are  substituted  into  the  general  form  of  the  numerical 
equation  given  as  equation  (2-68)  in  order  to  eliminate  P(I,J-1)  and 
P(I-1,J)  from  the  general  equation.  This  results  in  the  following 
equation . 

P(I,J)  * (AIM1(I,J)  + AJM1(I,J)  - 1)  + P(I+1,J)  * AIP1(I,J) 

= -B(I,J)  + A0  * If  * AIM1(I,J)  + AR  * |£  * AJM1(I,J) 
c)y  dr 

- P(I,J+1)  * AJP1(I,J)  (2-73) 

where  for  point  I,J  the  right  hand  side  is  completely  known.  For  each  of 
the  cavity  points  (/^  in  figure  2-4),  a similar  equation  can  be  written; 
however,  only  the  radial  derivative  3p/3r  can  be  specified  for  these 
nodes.  For  example,  for  the  first  node  to  the  right  of  the  corner  node 
(I+1,J)  the  appropriate  boundary  equation  becomes: 
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P(I+1,J)  * (AJM1(I+1,J)  - 1)  + P(I+2,J)  * A1P1(I+1,J)  + P(I,J) 

* AIM1(I+1,J)  = -B(I+1,J)  + AR  * * AJM1(I+1,J) 

OT 

- P{I+1,J+1)  * AJP1(I+1,J)  (2-74) 

where  the  right  heind  side  is  completely  known.  On  the  last  of  the  cavity 
nodes  to  the  right,  the  term  P(I+2,J)  * AIP1(1+1,J)  will  move  to  the 
right  hand  side  of  the  equation  as  a known  quantity. 

Thus,  the  corner  node  plus  the  cavity  nodes  to  its  right  become  a 
system  of  n equations  and  n unknowns  which  depend  on  all  of  the  sur- 
rounding nodes  on  adjacent  rows.  This  set  of  equations  can  be  solved 
simultaneously;  but  in  some  cases  the  number  of  cavity  nodes  in  a row 
becomes  large,  so  relaxation  is  used  instead.  Starting  at  the  right 
most  cavity  point/^  and  using  an  equation  such  as  (2-74) , the  pres- 
sure P(I+1,J)  is  solved  for.  Then  the  pressure  at  the  next  point  to 
the  left  is  solved  for  and  the  process  repeated  until  the  corner  point 
pressure  is  calculated.  This  procedure  is  repeated  several  times  until 
the  comer  point  pressure  stabilizes. 

Pre':,sure  values  for  the  corner  point  P(I,J)  and  the  cavity  points 
to  the  right  are  now  known.  Then  equation  (2-71)  is  applied  to  de- 
termine P(I-1,J)  (point  @ in  figure  2-4). 

P(I-1,J)  = P(I,J)  - A0  (2-75) 

Now  if 

P(I-1,J)  < PC  and  P(I,J)  > PC  (2-76) 

then  the  corner  point  position  is  acceptable.  That  is,  the  second  of 
the  downstream  boundary  conditions,  p = p^,  is  satisfied  as  close  as 
possible  within  A6.  Now  if 

P(I,J)  > PC  P(I-1,J)  > PC  (2-77) 

the  corner  point  is  shifted  one  to  the  left  and  the  entire  process  of 
computing  the  value  for  P(I,J)  is  repeated  with  the  new  corner  point 
position. 

If 

P(I,J)  < PC  and  P(1-1,J)  < PC  (2-78) 
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then  the  corner  point  is  moved  one  to  the  right  and  the  entire  process 
is  repeated. 

If  condition  (2-77)  shows  up  and  the  point  is  moved  to  the  left 
one  node  and  then  condition  (2-78)  shows  up,  or  vice  versa,  this  indi- 
cates that  condition  (2-76)  cannot  be  satisfied  within  A0  and  that  i 

either  of  these  positions  can  be  considered  sufficiently  close  to  the 
correct  solution.  (This  explains  why  in  the  results  some  of  the  nodes 
adjacent  to  the  cavity  have  small  negative  pressures.)  In  order  to 
help  stabilize  the  solution,  the  corner  node  is  moved  back  to  its  orig- 
inal position  when  this  situation  occurs.  Otherwise,  the  point  will  I 

shift  back  and  forth  one  node  with  each  iteration. 

5)  Once  the  corner  point  position  is  selected,  relaxation  begins  i 

just  to  the  right  of  the  last  boundary  point  determined  by  the  special 
condition  (point(]^in  figure  2-4)  where  general  equation  (2-68)  applies. 

Then  the  remainder  of  the  pressures  in  the  row  are  relaxed  and  the  iter- 
ation is  completed  row  by  row. 

Boundary  point(^  and  those  to  its  right  in  figure  2-4  are  treated 
in  the  same  way  except  that  there  is  no  comer  point  for  this  row.  An 
equation  of  a type  similar  to  equation  (2-74)  is  set  up  for  each  of 
these  points  and  the  pressures  along  this  row  are  relaxed  accordingly. 

Convergence  j 

Normally  the  convergence  criterion  used  for  an  over  relaxation  solu- 
tion is  the  maximum  amount  of  change  in  pressure  at  some  point  between 

i 

one  iteration  and  the  next.  When  this  quaintity  becomes  sufficiently  | 

small,  convergence  is  considered  to  have  been  achieved.  This  criterion  j 

does  not  work  for  the  solution  method  above  because  with  each  iteration  j 

small  changes  in  the  downstream  boundary  position  relative  to  a mecui  | 

position  continue  to  occur.  Thus,  small  perturbations  are  continually  j 

introduced  into  the  pressure  distribution  and  the  criterion  cannot  be 
easily  satisfied.  ] 

Now  it  has  been  estciblished  that  if  the  boundary  is  fixed  at  its  ] 

i 

mean  position,  the  solution  will  converge  using  the  above  criterion  and  j 

the  boundary  conditions  will  be  satisfied  to  an  adequate  degree  of  approxi- 
mation. However,  to  do  this  for  each  solution  requires  many  iterations 
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because  it  takes  a large  number  of  iterations  to  converge  to  an  approxi- 
mate boundary  shape  which  can  be  used  to  establish  the  average  position 
cind  even  more  iterations  to  converge  to  the  final  solution.  Thus,  it  was 
decided  to  seek  a different  type  of  convergence  criterion  based  on  func- 
tion rather  than  individual  pressures. 


The  parameters  needed  from  the  solution  for  studying  the  effects 

of  waviness  are  lecikage,  load  support,  and  the  first  harmonic  of  the 

variation  of  load  support  (p  and  p needed  for  deflection  calcu- 

®1  “l' 

lations.  It  was  found  that  these  parameters  are  not  greatly  affected  by 
the  small  shifts  in  the  boundary  position,  and  in  fact  display  a very 
small  oscillation  about  a mean  value.  Figure  2-8  shows  the  behavior  of 
the  total  load  support  W as  a function  of  the  iteration  number  for  three 
different  over  relaxation  factors  (based  on  the  Findlay  case  discussed 
later).  The  oscillation  becomes  quite  small  after  some  40  iterations. 
Figures  2-9  and  2-10  show  the  behavior  of  the  parameters  p and  p . 

ai 

These  too  reduce  to  small  variations  about  a mean.  Since  p and  p, 

ai  b^ 

reach  a small  oscillation  at  about  the  same  time  as  the  total  load 
support,  the  behavior  of  only  W was  used  for  a convergence  criterion. 
Letting  i equal  the  number  of  the  iteration , the  convergence  criterion 
used  was 


< e 


where 
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(2-79)  I 
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Thus  a current  average  of  three  points  is  compared  to  an  older  average 
of  5 points.  For  best  results  in  converging  on  load  support,  e = 0.001 
was  used. 


In  some  rare  cases  it  was  found  that  the  above  criterion  could  be 
satisfied  but  that  the  solution  had  not  at  all  converged.  This  generally 
occurred  when,  during  over  relaxation,  a temporary  asymptotic  value  such 
as  shown  in  figure  2-8  for  the  0 = 1.5  curve  was  reached.  In  such  cases 
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Figure 


it  was  found  that  there  would  be  a large  discrepancy  between  inside  and 
outside  leakage.  Thus,  an  additional  convergence  criterion  was  added 

< 0,05  (2-82) 

Average  values  over  5 iterations  were  used  for  eind  Q^.  A large  dif- 
ference was  allowed  because  leakage  agreement  is  not  generally  good  at 
the  small  number  of  nodal  points  used  in  most  of  the  solutions  (see 
later  discussion  on  error) . 

As  to  the  speed  of  convergence,  for  the  isotropic  base  case  pre- 
sented later,  a total  of  20  successive  over  relaxation  iterations  was 
required  to  converge  to  the  limits  indicated.  In  order  to  converge  on 
load  support,  a total  of  4 of  these  solutions  were  required  (the  ini- 
tial guess  was  good) . Each  solution  required  fewer  successive  over 
relaxation  iterations  because  the  final  pressure  distribution  from  the 
previous  solution  was  used  for  the  initial  guess  for  the  next  solution. 
The  total  CPU  computer  time  required  for  this  solution  was  73  seconds. 

Although  the  primary  objective  of  rapid  convergence  was  achieved 
for  this  example  and  many  other  cases,  a greater  number  of  iterations 


are  needed  for  some  cases,  particularly  for  low  pressures  p^. 


SOLUTION  FOR  SMOOTH  FACES 

In  order  to  compare  the  solution  obtained  by  the  methods  herein  to 
that  obtained  in  the  first  report  [1]  and  that  obtained  by  Findlay  [18] , 
a solution  for  smooth  surfaces  with  two  waves  was  first  carried  out. 

The  dimensionless  parameters  chosen  were 

p = 0.012 
o 

p =0.0 
c 

= 0.0 

= 0.8  (2-83) 

h =0.5 
o 

= 0.05 

n =2 
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The  mesh  size  was  60  points  in  the  6 direction  emd  12  points  in  the  r 
direction. 


For  the  solutions  for  smooth  faces,  the  solution  method  used  was 
slightly  different  than  that  described  previously.  First  a constant 
number  of  120  iterations  was  used  as  opposed  to  convergence  on  load  sup- 
port. It  was  found  that  120  iterations  were  more  than  enough  for  the 
solution  to  reach  a small  stable  oscillation  about  a mean  (see  figure 
2-8) . Second,  the  special  boundary  point  relaxation  feature  described 
previously  was  not  used.  The  special  boundary  equations  were  relaxed 
along  with  all  other  node  equations.  These  two  features  were  incor- 
porated into  the  program  after  the  smooth  surface  results  were  obtained, 
so  that  all  the  results  for  rough  surfaces  given  later  do  use  the  corn- 
complete  method  as  described. 


Since  the  solution  was  set  up  for  rough  surfaces  with  c ^ 0, 
to  accomplish  solutions  for  smooth  surfaces  it  was  necessary  to  use  an 
arbitrary  value  of  c to  obtain  dimensionless  constcints.  To  make  the 
surfaces  effectively  smooth. 


E(H)  = h 

3 3 

E(H-^)  = h 


etc. , for  all  h. 

The  resulting  pressure  distribution  for  one  of  the  two  waves  is 
shown  in  figure  2-11.  The  extent  of  the  cavity  is  indicated  by  the  zeroes. 
Minimum  film  thickness  occurs  at  90°  and  maximum  film  thickness  occurs 
at  180°.  The  cavity  begins  immediately  after  the  point  of  minimum  film 
thickness  and  ends  just  after  the  point  of  maximum  film  thickness. 

Pressure  builds  up  in  the  converging  region  as  expected.  Looking  at  the 
outside  two  rows,  a radially  inward  flow  occurs  between  the  7th  and  44th 
point,  and  an  outward  flow  occurs  over  the  remainder  of  the  outside. 

The  net  flow  is  radially  inward.  Leakage  and  load  support  for  this 
solution  are 
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Figure  2-11. 


Pressure  Distribution  for  Smooth  Surfaces 
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W = 0.00757 


= -0.00370 

Q = -0.00407 
o 


(2-85) 


The  difference  between  the  inside  and  outside  leakage  represents  error 
in  the  solution  and  this  is  explored  later. 

Comparison  to  Previous  Results 

To  verify  the  computer  program  and  solution  method,  the  results 
obtained  using  the  solution  herein  can  be  compared  to  those  obtained  by 
Findlay  [18]  and  to  those  obtained  using  a more  elaborate  program  pre- 
sented in  a previous  report  [1].  Findlay's  [18]  non  dimensional  input 
pareune  te  r s are : 


A = 


6 nt*) 
^ref 


2 

/r  \ 

I — I 

\h.  I 
' 1 / 


^ ^amplitude'^^mean 


(2-86) 


(2-87) 


^^o^  p _ 
F ^ref 


F ^ 


Pi 


ref 

= 


(Pc) 


F Pref 


Findlay's  numerical  solution  parameters  are: 


w 

_ Findlay 

" 2 
P 

ref  o 


_ i2nQ 

*r  ^3 

h.p  , 
1 ref 


(2-88) 


(2-89) 


(2-90) 


(2-91) 


(2-92) 


The  Findlay  solution  and  the  solution  herein  can  be  compared  by  as- 
suming certain  arbitrary  constants.  Let 
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h.  = 20  u in 
1 


Findlay  mean  film  thickness 


p , = 14.7  psia 
ref 


(2-93) 


c = 40  11  in 

Also  consider  the  Findlay  case  where 


(P^)  = 2 

F 


(Pi)F  = 1 


(2-94) 


(Pc)  = 1 

^ F 


A = 2000 

Then,  from  the  definition  of  A 


nwr^  = 1.96  • lo"®  lb 
o 


(2-95) 


The  pressures  used  in  the  present  solution  are  gauge  pressures , so 


(p  ) = 2 X 14.7  psia 


(p.)  = 14.7  psia 

^ F 


(2-96) 


p =14.7 
o 


p.  = 0 


and  from  equations  (2-93) , (2-94)  , and  (2-95)  , 


p = 0.012 

o 


p.  = 0.0 
1 


(2-97) 


p =0.0 

c 


As  given  previously,  the  waviness  equivalent  to  the  Findlay  case  is  given  by 


h.  = h 
1 


(2-98) 


e = 


o 

h 

n _ 

h 

n 

h. 

1 

h 

o 

II 

h. 

1 

c 

c 

(2-99) 


0.5 


(2-100) 
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Also, 


r = 1 
o 


r.  = 0.8 
1 


(2-101) 


Findlay's  load  support  is  based  on  absolute  pressure.  Load  support 
W used  herein  is  based  on  gauge  pressure.  The  necessary  conversion  is 


3 - 

r nuW 

C-  - Tid  - r^f)  = 

L 1 2 

c P £ 
ref 


For  lecikage  the  necessary  conversion  is 


2 

12nr  wcQ 

= 3'  ■ 

^i  ^ref 


(2-102) 


(2-103) 


Table  2-1  shows  a comparison  between  the  Findlay  results  and  the 
present  solution  for  several  different  cases.  In  general  the  agreement 
is  good,  especially  considering  the  fact  that  the  Findlay  values  were 
talcen  from  small  curves  presented  in  the  paper. 

A comparison  of  the  cavity  shape  obtained  in  reference  [1]  to  the 
present  solution  for  case  2 of  table  2-1  is  made  in  figure  2-12.  The 
present  solution  is  referred  to  as  the  3rd  approximation  in  the  figure. 
Comparison  to  the  reference  solution  is  good.  The  dots  plotted  repre- 
sent the  zeroes  of  cavity  pressure  which  outlines  the  cavity  as  shown  in 
figure  2-11.  The  dashed  line  was  plotted  just  to  the  outside  of  these 
points.  Since  the  Findlay  solution  agrees  closely  with  the  reference 
solution  (see  reference  [1]),  the  present  solution  then  also  agrees 
closely  with  the  Findlay  solution. 

The  above  results  were  considered  as  being  an  adequate  independent 
checlt  on  the  solution  method  developed  herein. 


Limit  Studies 

To  establish  a grid  size  necessary  for  an  accurate  solution,  studies 
were  made  using  various  mesh  sizes.  Table  2-2  shows  the  results  of 
these  studies.  For  the  double  precision  method,  double  precision  was 
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Figure  2-12.  Cavity  Shape  Comparison 


Table  2-2 


Limit  Studies 


Mesh  Size 
IMAX  X JMAX 

Single  Precision 

Partial 

Double  Precision 

w 

2i 

(1) 

2i 

2o 

50  X 10 

.007488 

-.003675 

-.004087 

.007488 

-.003675 

-.004087 

60  12 

.007672 

-.003672 

-.003791 

.007672 

-.003672 

-.003791 

75  X 15 

.007690 

-.003592 

-.003788 

i 

i 

100  X 20 

.007767 

-.003618 

-.003795 

1 

1 

Po  = 

0.012 

r . = 
1 

00 

o 

Pi  = 

O 

o 

n = 
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Pc  = 

0.0 

h = 
o 

0.5 

ii  = 
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0.05 

120 

iterations 
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applied  to  only  the  finite  difference  operations,  i.e.,  the  entire 
program  was  not  written  in  double  precision.  From  the  tables  it  can 
be  seen  that  load  values  approach  a limit.  The  single  precision 
values  are  equal  to  the  double  precision  values,  so  that  single  pre- 
cision can  be  used  to  acceptable  accuracy.  The  inside  and  outside  leak- 
age values  approach  each  other  as  the  mesh  is  made  finer.  This  is  eui 
indication  that  the  solution  is  becoming  quite  accurate  at  the  smaller 
mesh  size.  One  of  the  factors  that  causes  a significant  difference  be- 
tween inside  and  outside  leakage  at  large  mesh  sizes  is  that  leakage  at 
the  inside  boundary  is  defined  by  only  a few  points  because  of  the 
cavity  (see  figure  2-11) . Also,  the  net  leakage  across  the  seal  re- 
presents a small  fraction  of  the  total  circumferential  flow. 

Because  of  the  better  accuracy,  the  grid  size  chosen  for  further 
studies  was  the  100  x 20.  This  was  reduced  to  100  x lo  for  seals  where 
r^=  0.9  since  this  provides  the  same  Ar. 


proximate  Solutions 


The  solution  method  developed  herein  and  termed  the  3rd  approxi- 
mation is  based  on  the  application  of  correct  boundary  conditions  on  the 
cavity  such  that  flows  balance  and  inside  and  outside  leakages  match  to 
limits  imposed  only  by  the  finite  difference  nature  of  the  solution. 


Even  this  method  requires  large  amounts  of  computer  time  to  obtain, 
for  example,  a time  dependent  or  wear  study,  so  it  is  of  considerable 
interest  to  compare  faster  and  more  approximate  methods  to  this  more 
exacting  method  to  assess  the  error.  Two  other  such  methods  have  been 
so  evaluated.  They  are  both  commonly  used  in  lubrication  problem  solu- 
tions. The  first  of  these,  termed  Approximation  1,  is  to  simply  solve 
for  the  pressure  distribution  completely  disregarding  the  cavity.  Then 
all  pressures  which  are  less  than  the  cavity  pressure  are  set  to  the  cav- 
ity pressure.  The  second  of  these  approaches  is  to  enforce  the  condition 


(2-104) 


along  the  entire  cavity  boundary.  This  is  accomplished  by  setting  all 
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pressures  less  than  the  cavity  pressure  equal  to  the  cavity  pressure  as 
soon  as  they  occur  (same  as  for  the  upstream  boundary  of  Approximation 
3). 

Both  of  these  methods  are  numerically  very  fast.  But  table  2-3 
shows  a large  error  in  the  leakage  since  flow  across  the  cavity  does  not 
match.  There  is  also  some  error  in  total  load  support.  Figure  2-12 
shows  how  the  cavity  shapes  compare  for  the  three  different  approxima- 
tions. The  first  and  second  approximations  compare  favorably  at  the  up- 
stream boundary  but  depart  radically  toward  and  at  the  downstream 
boundary . 

Of  interest  in  later  chapters  where  deflection  of  the  seal  ring 
is  studied  is  the  error  introduced  into  the  load  distribution  function 
p(0) . p(6)  has  been  plotted  for  the  three  approximations  and  the  solu- 

tion in  reference  [1]  in  figure  2-13.  The  first  and  second  approxi- 
mations create  a large  difference  which  will  lead  to  a significant 
error  in  the  first  harmonic  content  of  p(6) . Thus,  for  cuiy  precise 
evaluation  of  deflection,  the  first  and  second  approximations  are  not 
satisfactory.  Even  so,  as  discussed  later,  in  order  to  obtain  some 
preliminary  time  dependent  results , it  has  been  necessary  to  resort  to 
the  second  approximation. 


Convergence  on  Load  Support 

To  be  most  meaningful,  solutions  obtained  must  provide  equilibrium 
between  the  seal  closing  force  W*  and  the  total  opening  force  W.  Now 


W = ^<Po'  Pi'  Pm' 


r. , n,  h , h ) 
1 on 


(2-105) 


and 

W*  = f(p  , p.,  B,  r.)  (2-106) 

O 1 1 

where  B is  the  balance  ratio.  All  of  the  variables  above  but  h are 

o 

established  by  the  conditions  of  operation  and  the  waviness  on  the  seal. 

Thus,  in  order  to  make  W = W*,  h^  must  be  varied  until  this  condition 

is  met.  Thus  h becomes  a dependent  variable  in  the  problem,  and  leakage, 
o 

percent  hydrodynamic  load  support,  and  friction  all  become  dependent 
upon  p^,  p^,  p^^,  r^,  n,  h^. 


i 
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Table  2-3 


Comparison  of  Approximate  Solutions 
Equation  (2-83)  Values 


PROGRAM 

MESH 

W 

2i 

2o 

®AV 

APPROX  1 

o 

1— 1 

X 

CD 

.0104 

-.0100 

-.0055 

-.0077 

APPROX  2 

O 

X 

GO 

,0107 

-.0101 

-.0004 

-.0053 

APPROX  3 

48  X 10 

.0074 

-.0032 

-.0041 

-.0036 
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To  find  the  value  of  h^,  a particular  value  is  guessed  and  the 
secant  method  of  finding  a root  is  used.  This  requires  that  many  val- 
ues of  W(h^)  be  found  in  order  to  establish  an  equilibrium  solution. 

All  solutions  presented  hereafter  are  equilibrium  solutions. 


SAMPLE  SOLUTIONS  FOR  LONGITUDINAL  AND  ISOTROPIC  ROUGHNESS 

Returning  now  to  the  original  problem  of  interest,  that  of  the 
effect  of  surface  roughness,  a base  case  for  study  must  first  be  es- 
tablished. Considering  some  typical  values  for  a hot  water  seal  (a 
heavily  loaded  low  viscosity  application) , the  following  non  dimen- 
sional case  was  established  as  a base  1: 


p = 0.02 

0 

p.  = p = 0.0 

1 c 

p =2.0 
m 

p = 0.1  p 
s m 

B ^ 0.75 

r.  = 0.9 
1 

n 3 

r /c  = 0.5  X 10^ 
o 


(2-107) 


h3  = 0.2 

These  non  dimensional  values  correspond  to  the  following  dimensional 
varicUoles 


S'] 


1 


(6.37»10”^  lb*s/ft^)  (2-108) 


I 

peak  roughness  height,  if  the  sur- 
face roughness  were  distributed  according  to  the  parabolic  distribution, 
equation  ( 2-5  ),  then  the  CLA  roughness  is  given  by  0.27c  so  CLA  = 0.27 
ym  (11  yin) . This  is  a fairly  typical  roughness  for  an  as  run  seal  face, 
and  it  will  be  shown  later  that  smaller  values  of  roughness  are  realistic 

when  this  topic  is  discussed  in  detail  in  chapter  5.  i 

} 

The  value  for  p corresponds  to  the  compressive  strength  of  a typ-  | 

- I 

ical  carbon  face  material,  p = 0.1  p was  established  so  that  the  co- 

s m 

efficent  of  friction  will  approach  0.1  if  the  load  becomes  completely 
supported  by  asperity  contact.  B = 0.75  is  a typical  value  for  the  bal- 
ance ratio  for  a high  pressure  seal. 

h^  = 0.2  ym  (8  yin)  represents  a small  amount  of  waviness  that  might 
occur  in  a seal  due  to  various  distortions  or  manufacturing  imperfections. 

As  equation  ( 2-59)  shows,  this  is  the  net  waviness  at  the  seal  face. 

The  amount  of  initial  waviness  required  to  produce  this  amount  of  wavi- 
ness will  be  larger  and  depend  on  the  stiffness  of  the  seal  rings.  This 
subject  is  taken  up  in  detail  in  chapter  5. 

For  longitudinal  roughness  the  pressure  distribution  resulting  for 
the  base  values  (2-107)  is  shown  in  figure  2-14.  The  *'s  indicate  the 
roughness  interference  region.  The  points  of  minimum  film  thickness  are 
at  the  top  and  bottom  columns.  The  cavity  shape  is  somewhat  typical  of 
many  of  the  solutions  obtained.  Generally  speaking,  the  film  pressure 
reaches  the  cavity  pressure  at  some  point  just  after  minimum  film  thick- 
ness is  reached  and  before  entering  the  non  interference  region.  The  up- 
stream cavity  boundary  is  usually  very  close  to  being  a radial  line  as 
shown.  For  roughness  interference  cases,  the  cavity  will  often  extend 


p =2.87  MPa 
o 

c = 1.0  ym 

p = 250  MPa 
m 

n = 3.05*10"^ 

u)  = 188  r/s 

r = 50  mm 
o 

hj  = 0.2  ym 
Recalling  that 


(416  psi) 

(39  yin) 
(36000  psi) 
Pa»s  (93*0  water) 
(1800  RPM) 
(1.97  in) 

(8  yin) 

c is  one  half  the 
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Figure  2-14.  Typical  Pressure  Distribution  for  Longitudinal  Roughness 
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to  the  top  row  as  shown.  At  higher  outside  pressure  the  cavity  will  move 
inward  and  become  smaller  as  expected.  The  full  film  is  reestablished 
along  the  downstream  boundary  as  shown  in  figure  2-14.  For  solutions 
obtained  so  far,  the  downstream  cavity  boundary  does  not  enter  the 
roughness  interference  region. 

For  the  solution  shown  in  figure  2-14  the  mean  film  thickness  j 

ii  =1.1.  With  ii  =0.2,  this  indicates  a maximum  interference  of  0.1 
o a ^ 

at  the  point  of  minimum  film  thickness  under  equilibrium  conditions.  j 

This  represents  a very  small  amount  of  interference  because  99.7%  of  the  | 

load  is  carried  hydrodynamically.  Thus,  even  though  this  seal  contacts  j 

i 

and  can  be  expected  to  wear,  hydrodynamic  load  support  carries  most  of  i 

the  load.  Friction  coefficient  p = 0.00165  and  Q = 0.0583  which  cor-  ] 

responds  to  a leakage  of  1.6  cc/min  for  the  values  cited  previously.  | 

j 

For  isotropic  roughness,  the  use  of  the  previous  parameter  values  | 

results  in  no  cavitation  amd  consequently  no  hydrodynamic  load  support.  ■ 

The  load  is  supported  by  hydrostatic  pressure  and  asperity  contact.  | 

Thus,  in  order  to  study  a typical  case  for  isotropic  roughness,  the  j 

roughness  value  used  was  halved  and  the  following  parameter  values  were  j 

established  as  base  case  2: 

c = 0.5  ym 

p = 0.005 
o 

p = 0.5  (2-109) 

m 

r /c  = 1 X 10^  : 

o I 

h = 0.2  i 

a ’ 

with  other  values  remaining  the  same.  • 

The  resulting  pressure  distribution  for  isotropic  roughness  is  shown 
in  figure  2-15.  The  cavity  has  a shape  similar  to  that  of  the  smooth 
surface  case  of  figure  2-11,  although  somewhat  smaller.  Again  cavita- 
tion begins  immediately  after  the  point  of  minimum  film  thickness.  i 

The  results  show  that  the  fraction  of  fluid  film  load  support  is  | 

72.5  percent.  However,  only  a small  part  of  this  is  due  to  hydrody-  j 

I 

namic  effects.  The  fraction  of  hydrostatic  load  support  for  parallel  | 

faces  is  given  by  | 
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Figure  2-15.  Typical  Pressure  Distribution 
for  Isotropic  Roughness 
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For  the  present  case  where  r^  = 0.9  and  B = 0.75,  W/W*  = 0.69,  just 
slightly  below  the  total  of  72.5  percent.  This  fraction  changes  over  a 
wide  range  as  shown  in  the  next  chapter. 

The  coefficient  of  friction  for  the  isotropic  case  is  0.031. 
Lealcage  Q = 0.0296  inward.  For  the  given  parameters  Q = 0.42  cc/min. 


RADIAL  TAPER  MODEL 

Seal  performance  is  significantly  affected  by  small  values  of 
radial  taper.  To  evaluate  such  effects,  a radial  variation  in  film 
thiclcness  can  be  added  to  equation  (2-58  ) as  follows: 

00 

h = ii  (r)  + y (ii,  + h,  ) cos  njG 

o la.  da. 

3=1  3 3 

oo 

+ y (h.,  + h-,  ) sin  nj6  (2-112) 

. , ID  . OD  . 

3=1  3 3 

It  is  assumed  that  the  radial  variation  is  independent  of  6,  although  the 
general  case  of  h(r,0)  can  be  solved. 

Equation  (2-20)  applies  to  a general  case.  To  solve  for  a radial 
variation,  it  is  necessary  only  to  include  radial  film  thic)cness  varia- 
tion effects  in  the  finite  difference  representation. 

Of  particular  interest  is  the  effect  of  a uniform  radial  taper 
superimposed  over  an  nth  harmonic  waviness.  In  this  case  film  thickness 
becomes 

h = h + ii  (r)  + h cos  n9  (2-113) 

or  n 


! 


i 
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”1 

j’ 

where  I 

h^(r)  = (r  - Doi  (2-114)  i 

ar 

a = — (2-115)  j 

° 1: 

where  a is  the  angle  of  radial  tilt  of  the  seal.  A positive  a indicates 

that  the  film  thickness  is  becoming  smaller  (converging  for  an  outside 

pressurized  seal)  when  moving  radially  inward.  ; 

Results  based  on  the  radial  taper  model  for  isotropic  roughness  are 
presented  in  chapters  3 and  4. 


63 


CHAPTER  3 


PARAMETER  STUDIES 

Now  that  the  basic  model  and  method  of  solution  have  been  developed, 
it  is  possible  to  make  detailed  studies  of  the  effects  of  various  param- 
eters on  seal  performance. 

The  parameters  which  are  considered  as  the  independent  variables 
over  which  a seal  designer  might  have  some  control  are  surface  rough- 
ness, face  width,  waviness  amplitudes,  number  of  waves,  and  compressive 
strength  of  the  material.  Other  parameter  values  such  as  seal  diameter, 
viscosity,  speed,  and  sealed  pressure  are  usually  fixed  by  the  sealing 
requirements.  Thus,  these  studies  will  concentrate  on  examining  the 
effect  of  the  independent  variables  while  holding  other  parameters 
fixed. 

Of  the  independent  variables,  the  number  of  waves  requires  further 
discussion,  n was  fixed  at  a value  of  three  for  most  studies  herein. 

This  choice  is  somewhat  arbitrary,  but  results  based  on  n = 3 are 
expected  to  be  typical  of  what  can  be  expected  at  a higher  number  of 
waves.  This  number  is  typical  of  the  number  of  waves  found  in  actual 
seals. 

As  for  the  waviness  amplitude,  these  studies  are  based  purely  on 
the  net  amount  of  waviness  at  the  seal  faces.  The  relationship  between 
the  net  waviness  and  the  initial  waviness  required  to  produce  this  net 
waviness  is  exeunined  thoroughly  in  chapter  5. 

PERFORMANCE  PARAMETERS 

There  are  three  performance  parameters  of  importance.  These  are 
the  fraction  of  the  load  supported  by  fluid  pressure,  leakage,  and 
coefficient  of  friction.  The  fraction  of  load  supported  by  fluid 
pressure  is  important  because  from  equation  (5-15)  , 

‘^^w  (1  - H )W* 

__o  ^ ^ 

dt  2ti(1  - r.) 
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where 


= average  wear  rate 


= fraction  of  load  supported  by  fluid  pressure. 

Thus,  as  approaches  unity,  the  wear  rate  goes  to  zero.  This  is  of 
course  very  important  to  seal  performance . As  pointed  out  in  chapter  2 , 
has  a significeint  value  (.50  to  .70)  due  to  hydrostatic  load  support 
alone.  Thus,  any  significant  reduction  in  wear  rate  compared  to 
a seal  operating  with  parallel  faces  requires  that  become  signifi- 
cantly larger  than  the  hydrostatic  fraction.  Leakage  is  important  in 
that  it  must  generally  be  held  to  some  minimum  value.  Friction  is 
also  important.  Not  only  does  high  friction  cause  energy  losses,  but 
perhaps  more  importantly,  the  concentrated  heating  effect  caused  by 
high  friction  can  cause  undesirable  and  uncontrollable  thermal  dis- 
tortion in  the  seal  rings. 


LONGITUDINAL  VERSUS  ISOTROPIC  ROUGHNESS 

Parameter  studies  will  be  presented  for  both  longitudinal  rough- 
ness and  isotropic  roughness.  Now  to  put  the  significance  of  these 
results  in  perspective  in  relation  to  actual  seals,  it  is  concidered 
that  pure  longitudinal  roughness  is  an  ideal  situation  which  is  not 
likely  to  occur  in  any  practical  seal.  Since  the  fluid  film  load 
support  fraction  for  longitudinal  roughness  is  always  very  high 
(98  percent  plus) , the  longitudinal  results  may  be  thought  of  as  an 
upper  bound  on  performance. 

Now,  the  results  in  chapter  7 show  that  for  the  Christensen  isotropic 
roughness  model  the  flow  resistance  is  lower  than  predicted  by  the 
method  of  Patir  and  Cheng  [53].  This  means  that  the  hydrodynamic 
pressure  buildup  using  the  Christensen  isotropic  model  will  be  lower 
than  that  expected  based  on  the  possibly  more  exacting  theory  of 
Patir  and  Cheng  [53].  Thus,  the  results  presented  for  isotropic  rough- 
ness can  be  thought  of  as  being  a lower  bound  as  far  as  hydrodynamic 
load  support  is  concerned. 
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of  the  above  two  extremes,  it  is  thought  at  this  time  that  the 
isotropic  results,  even  though  the  load  support  may  be  underestimated, 
are  closer  to  reality.  The  reasoning  behind  this  is  that  actual  seals 
are  not  accurately  characterized  by  the  ideal  longitudinal  roughness 
in  spite  of  the  fact  that  such  faces  appear  to  have  such  longitudinal 
roughness.  There  is  in  fact  a considerable  roughness  component  in  a 
tangential  direction.  Roughness  in  the  tangential  direction  will  al- 
low considerable  side  leakage,  and  this  will  greatly  reduce  the  hydro- 
dynamic  pressure  buildup.  The  longitudinal  roughness  case  represents 
a practically  unachievable  roughness  condition.  In  view  of  this  con- 
clusion, most  of  the  studies  herein  have  been  based  on  the  isotropic 
model. 

LONGITUDINAL  ROUGHNESS  PARAMETER  STUDIES 


The  parameter  values  used  as  a base  point  for  parameter  studies 

for  longitudinal  roughness  are  the  same  as  given  in  equations  (2-107) 

and  (2-108).  The  results  shown  here  have  been  taken  from  reference  [54]. 

Figure  3-1  shows  the  effect  of  varying  the  quantity  nu)  with  all  other 

parameter  values  held  constant.  This  type  of  study  will  give  similar 

results  to  those  obtained  if  the  pressure  p were  varied  with  p held 

o m 

constant,  but  there  will  be  some  small  differences. 

Figure  3-1  shows  a friction  characteristic  that  is  typical  of  a 
transition  from  mixed  lubrication  to  hydrodynamic  lubrication.  There 
is  a point  of  minimum  friction  and  friction  increases  from  this  point 
as  no)  is  either  increased  or  decreased.  One  interesting  feature  of 
this  curve  is  that  minimum  friction  occurs  to  the  left  of  the  point 
where  some  roughness  interaction  begins.  That  is,  the  lowest  friction 
is  encountered  when  some  interference  exists.  From  this  curve,  it  can 
be  concluded  that  seals  will  continue  to  operate  primarily  hydro- 
dynamically  with  low  friction  even  though  there  is  considerable  rough- 
ness interference  (at  the  lowest  value  of  1/p  plotted  where  1/p  = 

o o 

12,5,  ii  = 1.000  so  the  interference  is  0.2).  However,  if  the  curve 
o 

to  the  left  were  continued,  the  friction  coefficient  would  rapidly 
increase  to  a limit  and  % hydrodynaunic  load  support  would  rapidly 
decrease.  Thus,  for  a given  initial  waviness  and  other  parameters. 


BASE  VALUE 


0.001 
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Figure  3-1.  Effect  of  Viscosity  on  Friction  and  Leakage 
Longitudinal  Roughness 


67 


COMPLETE 

SEPARATION 


r ■ — ^ 

there  will  be  a definite  value  of  nti)  below  which  hydrodynamic  load  sup- 
port will  no  longer  be  significant  and  wear  rate  cind  friction  will  be- 
I come  large. 

I Figures  3-2  and  3-3  show  the  effect  of  initial  waviness  on  various 

I seal  parameters.  Figure  3-2,  for  the  base  values  (2-107) , indicates  that 

' leakage  increases  with  increasing  waviness  as  expected.  The  hydro- 

! dynamic  load  support  is  greater  than  98  percent  for  the  entire  curve, 

and  friction  changes  very  little.  However,  as  the  waviness  decreases 

below  0.03,  the  curve  for  h shows  that  the  maximum  film  gap  becomes 

o 

zero.  This  means  that  the  entire  seal  face  is  in  the  roughness  inter- 
ference region.  There  will  be  no  flow  into  the  seal  (theoretically) 
and  the  hydrodynamic  pressure  field  will  become  indeterminate  because 
the  quantity  of  fluid  present  is  unknown.  In  a practical  sense,  hydro- 
dynamic  lubrication  will  cease,  and  the  load  will  become  completely 
supported  by  the  asperity  tops.  Thus,  there  is  a certain  waviness 
cimount  for  a given  set  of  parameters  below  which  there  will  be  no  hydro- 
dynamic lubrication. 

Figure  3-3  has  been  obtained  for  a case  with  two  times  the  rough- 
ness used  for  the  results  in  figure  3-2.  In  this  figure,  it  is  seen 

that  h does  significantly  affect  friction  and  % hydrodynamic  load 
sl 

support.  Noting  that  the  values  of  h plotted  represent  two  times  the 
waviness  of  those  values  in  figure  3-2 , it  can  be  observed  that  friction 
is  larger  and  % load  support  is  smaller  for  the  same  actual  waviness  h^. 

In  general,  other  factors  being  equal,  a smaller  roughness  height  pro- 
vides a greater  % load  support  and  reduced  friction.  Figure  3-3  also 
shows  that  hydrodynamic  load  support  collapse  will  occur  at  a larger 
value  of  waviness  than  for  the  lower  roughness  seal  of  figure  3-2. 

Conclusions  Based  on  Parameter  Studies  for  Longitudinal  Roughness 

Based  on  the  above  results,  several  conclusions  can  be  drawn  for 
the  case  of  longitudinal  roughness. 

1)  Even  in  heavily  loaded  or  low  viscosity  seals,  most  of  the 
load  may  be  supported  by  fluid  film  pressure.  Friction  and 
wear  are  reduced  accordingly. 
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Effect  of  Waviness  on  Friction  and  Leakage  - Longitudinal  Roughness 


2)  Some  amount  of  touching  or  roughness  interference  can  be 
expected  in  heavily  loaded  seals  even  though  most  of  the  load 
may  be  supported  hydrodynamically.  Coirplete  separation  is  not 
expected  until  the  load  is  considerably  reduced. 

3)  Although  only  a small  amount  of  waviness  is  required  to  cause 
a large  hydrodynamic  load  support,  the  effect  of  roughness  is 
that  when  waviness  reaches  a certain  lower  limit,  hydrodynamic 
support  ceases. 

4)  Increasing  roughness,  other  things  being  equal,  leads  to  a 
decrease  in  the  fraction  of  the  load  supported  hydrodynamically. 

5)  The  waviness  amplitude  has  only  a small  effect  on  the  fraction 
of  hydrodynamic  load  support  over  a wide  range  subject  to  the 
limitation  in  3)  above.  However,  leakage  does  increase  with 
increasing  waviness. 

The  above  conclusions  are  based  on  results  that  do  not  consider 
the  initial  waviness  required  to  obtain  the  net  waviness  values  used 
in  the  studies.  This  matter  is  considered  in  Chapter  5.  However, 
since  the  required  waviness  is  quite  small  for  the  above  results  and 
given  that  such  small  amounts  of  waviness  in  seals  is  likely,  it  would 
follow  that  many  seals  would  be  99  percent  fluid  film  supported.  How- 
ever, comparing  wear  rates  of  actual  seals  to  those  predicted  by  the 
above  results  suggests  that  this  conclusion  is  not  generally  true — as 
pointed  out  in  the  comparison  between  the  longitudinal  and  isotropic  models. 
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ISOTROPIC  ROUGHNESS  PARAMETER  STUDIES 

For  the  isotropic  parameter  studies,  two  base  cases  are  used. 
Case  1: 

p = 0.02 

o 

p.  = p =0 
p = 2.0 

m 

h = Pm 

ha  = 0.2 

n = 3 


(3-2) 


r.  = 0.9 
1 

r^/c  = 0.5  • 10^ 


B = 0.75 


These  dimensionless  parameters  correspond  to  the  dimensional  values 
given  by  equation  (2-108).  For  case  2,  the  roughness  was  changed  to  one- 
half  of  that  used  in  the  cibove  case;  other  parameter  values  remained 
the  Scune. 

Case  2: 


p = 0.005 


Effect  of  Net  Waviness 

Figures  3-4  and  3-5  show  the  effect  of  net  waviness  on  seal  per- 
formance. For  case  1 (figure  3-4)  the  % fluid  film  load  support  in- 
creases only  slightly  above  the  hydrostatic  value  obtained  from 
equation  (2-111)  (69  percent  for  r^  = 0.9).  Friction  is  reduced  as  more 
of  the  load  is  supported  by  fluid  film  pressure.  Leakage  increases 
with  increasing  h^,  as  expected.  The  mean  film  thickness  also  increases. 
At  h^  = 1.0,  the  amount  of  penetration  into  the  rough  surface  is  0.34 

euid  at  ii  =0,  0.20. 

3 

In  figure  3-5  the  % fluid  film  load  support  increases  markedly  as 
waviness  increases.  When  compared  to  case  1,  this  result  shows  the 
importance  of  surface  roughness  on  hydrodynamic  pressure  buildup. 
Comparing  h,  = 0 to  h = 1.0  for  case  2,  the  hydrodynamic  effect  causes 
the  % fluid  film  load  support  to  increase  from  the  hydrostatic  value  of 
69  percent  to  92  percent.  At  hg  = 1.0  the  average  wear  rate  would  be 
reduced  to  25  percent  of  that  at  hg  = 0 (see  equation  (3-1)).  Leakage 
increases  a factor  of  ten.  Friction  is  reduced  to  one  third  of  the 
hj  = 0 value.  Thus,  these  results  show  that  a small  amount  of  waviness 
could  have  a large  effect  on  seal  performance  if  the  as-running  rough- 
ness is  small  enough. 


Effect  of  Speed  and  Viscositv 


In  order  to  study  the  effect  of  speed  and  viscosity  independently, 
the  dimensionless  parameter  p^  was  varied  while  keeping  p^  at  a given 
ratio  to  p^.  Figures  3-6  and  3-8  show  these  results  for  cases  1 and  2; 
(nw)  increases  moving  toward  the  left.  For  case  1,  significant  hydro- 
dynamic  effects  do  not  occur  until  nu)  is  nearly  two  times  the  base 
values.  Leeikage  Q is  shown  as  decreasing  with  increasing  no).  Actually, 
given  the  non-dimensional  character  of  Q,  if  oi  only  were  increased, 
then  the  actual  leakage  remains  constant  with  increasing  u.  If  n only 
were  increased,  then  the  actual  leakage  would  decrease  in  relation  to  Q. 


Figure  3-8  shows  that  even  with  the  small  roughness  value,  if  nw 
were  made  one-half  the  value  in  the  base  case  that  hydrodynamic  effects 
would  vanish.  Also,  only  after  nw  is  increased  a factor  of  four  or  so 
does  the  % load  support  begin  a sharp  upward  turn. 
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Figure  3-6.  Effect  of  (nw)  - Isotropic  Roughness  Case  1 
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Generally,  based  on  test  data  for  seals  presented  in  the  literature 
and  summarized  by  Pape  [30] , it  is  to  be  expected  that  friction  will  i 

increase  with  increasing  noi.  The  beginning  of  this  trend  is  shown  in  ! 

I 

figure  3-8.  A large  increase  would  occur  if  the  curve  had  been  ex- 
tended on  toward  the  left.  Figure  3-1  for  the  longitudinal  roughness  ! 

included  a wider  range  and  shows  the  expected  increase  more  clearly. 

Figures  3-7  and  3-9  show  the  effect  of  a variable  pressure  only.  ! 

All  other  parameters  are  held  fixed.  Now  comparing  figures  3-6  and  3-7, 
all  of  the  curves  are  similar  except  h^.  In  figure  3-7,  the  effect  of 
pressure,  the  actual  load  increased.  This  requires  that  the  depth  of 
penetration  increases  in  order  to  develop  the  needed  asperity  load 

support.  Thus,  h^  required  for  equilibrium  decreases.  In  figure  3-6,  | 

no  further  penetration  was  needed  because  the  load  was  constant,  so  j 

h remained  nearly  constant, 
o 

Comparing  the  leakage  curves  of  figures  3-8  and  3-9  shows  that  the 
leakage  under  increased  pressiire  at  p^  = 0.01  is  somewhat  less  for  the 
p^  study  than  for  the  nw  study.  This  occurs  because  the  flow  resistance 
becomes  greater  as  the  mean  film  thickness  becomes  smaller. 

Face  Width 

Figures  3-10  and  3-11  show  that  an  increase  in  face  width  increases 
the  % fluid  film  load  support  and  reduces  leakage.  Friction  drops  off 
accordingly,  h^  increases  with  increased  face  width  because  less  pene- 
tration is  required  to  develop  the  needed  asperity  load  support. 

Although  these  results  suggest  that  wider  faces  would  be  very  useful, 
other  practical  limitations  such  as  heat  generation  often  work  against 
using  a wider  face  seal. 

Number  of  Waves 

Figure  3-12  shows  the  effect  of  increasing  the  number  of  waves  for 
case  2.  Percent  load  support  is  increased  with  only  a very  small  in- 
crease in  leakage.  Friction  drops  off  accordingly. 
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Surface  Roughness 

While  the  effects  of  surface  roughness  are  clear  when  cases  1 and 
2 are  compared,  a separate  study  is  useful  to  illustrate  some  additional 
characteristics  of  surface  roughness.  Figure  3-13  shows  that  the  percent 
load  support  begins  to  increase  very  dramatically  for  roughness  values 
less  than  c = 0.5  w in.  Friction  drops  off  accordingly.  The  cxirve,  if 
extended,  will  give  100  percent  fluid  film  load  support  at  about 
0.25  pm  (10  y in.).  The  required  CLA  roughness  for  this  case  might 
be  as  low  as  0.07  pm  or  3 p in.  Thus,  the  value  of  the  as-running 
surface  roughness  required  for  complete  liftoff  is  quite  small. 

Due  to  the  dimensionless  character  of  Q,  actual  leakage  will  de- 
crease with  decreasing  roughness  even  faster  than  shown  for  Q.  Also, 

the  dimensional  h actually  decreases  as  surface  roughness  decreases, 
o 

Radial  Taper 

Figure  3-14  shows 'Seal  performance  as  a function  of  radial  taper 
for  Case  2.  The  results  show  that  even  a small  converging  taper  of 
a = 20,  which  is  equivalent  to  a total  taper  of  1 pm  (40  p in.)  across 
the  width  of  the  face,  is  sufficient  to  cause  almost  complete  fluid 
film  load  support.  The  leakage  increases  a factor  of  four  compared  to 
the  base  case.  So,  hydrodyneunic  effects  can  easily  be  overshadowed 
by  radial  taper  effects. 

On  the  other  hand,  for  a negative,  the  hydrostatic  load  support 
rapidly  diminishes  and  total  fluid  film  load  support  decreases.  Hydro- 
dynamic  effects  do  persist  as  shown,  but  only  a small  fraction  of  the 
load  is  carried  by  hydrodynamic  effects  because  the  average  gap  between 
the  faces  becomes  larger  with  increasing  divergence.  The  seal  becomes 
supported  primarily  by  the  wave  peaks  at  the  outer  radius.  Leakage  in- 
creases with  increasingly  negative  a. 

Conclusions  Based  on  Parcuneter  Studies  for  Isotropic  Roughness 

Bearing  in  mind  that  the  results  in  this  chapter  are  based  solely 
on  the  results  from  the  ideal  model  and  do  not  consider  such  practical 
factors  as  required  initial  waviness  and  obtainable  surface  roughness 
(to  be  included  in  chapter  5) , several  conclusions  can  be  reached. 
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1)  Concerning  the  hydrodynamic  load  support  and  concomitant 
reduction  of  wear  rate,  the  previous  results  show  that  a 
high  percentage  of  fluid  film  support  can  be  obtained  even 
for  heavily  loaded,  low  viscosity  seals  if  the  magnitude  of 
the  as-running  surface  roughness  is  sufficiently  small.  One 
hundred  percent  fluid  film  load  support  will  occur  for  smooth 

. 

surfaces,  given  even  a small  amount  of  waviness.  The  percent 
fluid  film  load  support  will  never  be  less  than  that  provided 
by  hydrostatic  pressure  which  is  50  percent  or  greater  for 
parallel  faces. 

2)  An  increase  in  waviness  amplitude  causes  the  percent  fluid 
film  load  support  to  increase.  Friction  decreases  and  leakage 
increases  with  increasing  waviness.  Relatively  small  aunounts 
of  waviness  are  needed  to  develop  significant  hydrodynamic 
effects  if  surface  roughness  is  low. 

3)  Increasing  nw  will  cause  an  increase  in  the  % fluid  film  load 
support.  Friction  is  not  greatly  affected  over  the  range  of 
values  calculated  but  does  begin  to  increase  as  nw  becomes 
larger.  Actual  leakage  decreases  with  an  increase  in  viscosity 
n and  remains  constant  with  increasing  u). 

4)  An  increase  in  pressure  p causes  a decrease  in  the  % fluid 

o 

film  load  support  and  an  increase  in  leakage.  Friction  is 
not  greatly  affected  over  the  range  of  interest. 

5)  An  increase  in  face  width  causes  the  % fluid  film  load  support 
to  increase  and  leakage  and  friction  to  decrease. 

6)  An  increase  in  the  number  of  waves  of  equal  magnitude  causes 
% load  support  to  increase  and  friction  to  decrease  while 
leakage  remains  nearly  constant. 

7)  A reduction  in  as-running  surface  roughness  causes  the  % load 
support  to  increase  and  friction  to  decrease.  Leakage  also 
decreases.  At  low  roughness  values,  100  percent  fluid  film 
load  support  can  be  achieved. 

8)  A radial  taper  can  cause  a large  difference  in  seal  performance 
compared  to  parallel  faces.  Hydrodynamic  effects  can  become 
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of  secondary  importance  for  a radially  convergent  taper.  A 

radially  divergent  taper  reduces  both  hydrodynamic  and  hydro 
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static  pressures. 
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CHAPTER  4 I 

COMPARISON  TO  EXPERIMENTAL  RESULTS  | 

I 

To  test  the  predictive  ability  of  the  theoretical  method  developed, 
predicted  results  must  be  compared  to  theoretical  results.  There  have 
been  three  previous  investigations  that  provide  some  useful  experi- 
mental results  for  seals  of  knovm  or  approximately  known  waviness. 

First,  Stanghan  - Batch,  and  Iny  [46]  present  some  leakage  results  for 
a wavy  seal  tested  using  oil  as  the  sealed  fluid.  These  results  are, 
therefore,  primarily  applicable  to  seals  with  parameter  values  p^  which 
are  very  small  where  full  hydrodynamic  load  support  is  expected.  Since 
this  is  not  the  range  of  interest  here,  these  results  have  not  been 
considered  further. 

Although  Pape's  work  [30]  is  again  for  lightly  loaded  seals,  he 
does  show  some  results  where  mixed  friction  appears  to  be  occurring. 

Pape's  friction  data  for  his  test  series  1 for  B = 1 has  been  repro- 
duced in  figure  4-1.  Friction  coefficient  decreases  with  increasing 

p ; however,  there  appears  to  be  a leveling  out  trend  in  the  data.  In 
o 

fact,  Pape  [30]  later  shows  that  based  on  a hydrodynamic  theory  for 
smooth  wavy  faces,  the  friction  coefficient  should  continue  to  decrease 
and  follow  the  general  downward  trend.  Thus,  the  leveling  trend  in 
the  experimental  data  would  appear  to  be  caused  by  mixed  friction 
(touching)  effects. 

To  test  this  hypothesis,  the  mixed  friction  lubrication  model 
developed  herein  was  used.  The  following  parameter  values  were  selected 
based  on  data  presented  in  Pape's  work  [30]. 


c = 0.30  vim  - based  on  reported  CLA  of  0.05-0.1  ym  before  running. 

r = 30  mm 
o 

r.  = 25  mm 


1 

n = 0.85  • 10”^  Pa's  (50°C  Kerosene) 

CO  = 126/s 

p = 3.59  " 10^  Pa  (invar  surface) 

h^  = 0.03  ym  (estimate  based  on  Pape's  measurements) 


(4-1) 


n = 3 


B = 1.0 


P 


s 


0.1  p (assumed  value  - no  data  available) 
m 


89 


90 


The  dimensionless  parameter  values  corresponding  to  the  above 


are 


p = 0.335 
m 

p.  = p =0 

1 c 

p = 0.1  p 
s m 

h =0.1 
3 

r.  = 0.83 
1 

r /c  = 1 * 10" 
o 


In  addition,  if  gauge  pressures  are  used  and  p^  = 0,  then 


= B 

Pape 


Pape  - ~ 2 B 

^ p r ^ B 

o o 


(4-2) 


(4-3) 

(4-4) 


Some  of  the  above  values  are  based  on  rough  estimates.  For  example, 

the  net  waviness  h^  is  not  )cnown  but  is  estimated  based  on  Pape's 

initial  waviness  measurements  and  displacements  during  operation . The 

effective  value  of  roughness  is  also  an  estimate  and  relies  on  an 

assumed  relationship  between  c and  CLA.  Pape's  data  were  generated  by 

varying  n,  w,  and  p^  in  his  actual  experiments.  In  this  comparison, 

only  p was  varied.  As  pointed  out  in  chapter  3,  some  difference  in 
o 

behavior  is  to  be  expected  for  the  model  depending  on  which  parameter 
is  varied.  The  weakest  assumption  in  the  data  is  that  p^  = 0.1  p^. 

The  effective  value  of  p^  is  simply  not  known  at  this  time.  This 
applies  to  all  results  herein  as  well  as  this  present  comparison.  The 
value  used  for  p^  greatly  affects  the  calculated  friction. 

In  spite  of  the  above  limitations,  the  predicted  results  shown  in 
figure  4-1  display  a behavior  which  is  quite  similar  to  the  Pape  experi- 
mental results.  The  leveling  trend  and  slopes  are  similar.  The  dashed 
curve  shows  the  effect  of  reducing  p^  to  0.3  of  the  value  above. 

Friction  at  the  higher  pressure  values  is  reduced. 

Although  the  comparison  in  the  Pape  data  certainly  suggests  that 
the  wavy  seal  mixed  friction  model  herein  may  be  valid,  more  precise 
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and  exhaustive  comparisons  are  needed  before  any  final  conclusion  can 
be  made  (see  chapter  8) , 

The  third  source  of  data  on  wavy  seals  is  provided  by  experiments 
conducted  by  Snapp  and  Sasdelli  [6] . The  focus  of  the  Snapp  and 
Sasdelli  work  was  on  the  effects  of  radial  taper.  It  is  useful  to  also 
compare  this  data  to  predicted  results  using  the  radial  taper  model 
developed  in  Chapter  2. 

Approximate  non  dimensional  parameter  values  for  the  Snapp  and 
Sasdelli  experiments  are: 

p^  = 0.005  at  100%  test  pressure 

p.  = p =0 

1 c 

p =0.32 
m 

n = 3 

= 11.25  (4-5) 

r /c  = 4.9  x 10^ 
o 

p = 0. 1 p 

'^s  m 

r.  = 0.961 
1 

B = 0.74 


One  of  the  interesting  features  of  this  seal  is  that  the  faces 
rotate  toward  each  other  in  a radially  divergent  fashion  as  sealed 
pressure  increases.  So,  depending  on  the  initial  taper,  various 
angles  of  convergence  or  divergence  occur  at  different  pressures.  The 
relationship  for  the  taper  is  given  by 

P, 


a = a.  . . , - --  x 1430  x 10”^  in/in. 
initial  0.005 


• 

— a.  . . - 1.40  • 10  p 

c initial  o 


(4-6) 


(4-7) 


For  the  theoretical  results,  it  is  assumed  that  no  wear  has  taken 
place;  i.e.,  the  radial  profile  is  a straight  line.  Table  4-1  was  pre- 
pared using  the  seal  model  herein  and  compares  the  effects  of  various 
amounts  of  initial  radial  taper.  Data  supplied  for  seals  C2  and  E2  in 
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the  reference  [6]  were  used.  No  significant  waviness  was  present.  The 
table  shows  that  seal  v 2 was  radially  convergent  throughout  the  pressure 
range.  Friction  is  very  low  and  leakage  is  large.  All  of  the  load  was 
supported  by  fluid  pressure.  In  contrast,  the  taper  for  seal  E2  is 
always  divergent  because  of  its  low  initial  taper.  Friction  is  very 
high  and  leakage  is  low.  Only  at  the  smallest  taper  is  the  % fluid 
film  load  support  significant. 

The  results  in  table  4-1  are  compared  to  the  experimental  results 
[6]  in  figures  4-2  and  4-3.  Figure  4-2  shows  that  predicted  leakage 
has  the  same  trend  as  the  experimental  results  but  is  much  larger.  It 
is  thought  that  the  reason  for  this  large  difference  is  that  the  tapers 
on  the  test  seals  are  slightly  worn.  A small  beveling  off  of  the 
initial  taper  would  reduce  the  leakage.  Both  the  theoretical  and 
experimental  leakage  for  seal  E2  are  very  low  and  were  not  plotted. 

Figure  4-3  shows  the  experimental  and  predicted  friction  torque 
for  seals  C2  and  E2. 

Table  4-1 

Predicted  Results  for  Snapp  and  Sasdelli  Seals  C2  and  E2  [6] , 

Effect  of  Initial  Taper 


Seal  C2  = 2460  x 10"®  • 


T , Torque 
(one  seal) 
in-lb 


Q 

2 

oz/min/in. 

Percent 

304 

1.61 

100 

531 

2.81 

100 

728 

3.85 

100 

791 

4.19 

100 

580 

3.07 

100 

.0118 

6.24-10"’’ 

0.7 

.0470 

2.49-  lO"** 

1.4 

.0833 

4.41-10"'' 

2.9 

.0889 

t 

o 

o 

6.7 

.0393 

2.08-10"'' 

20.3 
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Again  here,  the  theory  predicts  the  general  behavior  observed.  The 

predicted  friction  for  seal  E2  is  larger  than  the  experimental  values. 

This  is  probed>ly  due  to  the  assumption  that  p = 0.1  p . The  fraction 

s nt 

used  is  too  large  for  a carbon  versus  hard  face  combination. 

Table  4-2  compares  the  predicted  perform£mce  for  two  radially 
tapered  seals.  One  of  the  seals  has  a definite  planned  waviness  while 
the  second  does  not.  First  looking  at  the  results  for  seal  D2,  the 
initial  radial  taper  in  this  seal  is  small  enough  so  that  the  seal  goes 
from  a radially  convergent  to  a radially  divergent  operation  as  pressure 
is  increased.  Thus,  there  is  an  abrupt  decrease  in  leakage  amd  an 
increase  in  torque  as  sealed  pressure  is  increased  from  40  to  60  percent. 
The  experimental  leakage  is  compared  to  the  predicted  leakage  for  this 
case  in  figure  4-4.  This  transition  point  occurs  at  about  the  same 
place  for  both. 


Table  4-2 

Predicted  Results  for  Snapp  and  Sasdelli  Seals  A2  and  D2, 
Effect  of  Waviness 


Seal  A2 

^i 

= 720  X 

O 

1 

B 

11.25 

■ 

% Test 

Pressure 

a 

T, 

in.  lb 

Q 

Q 

oz/min/in. 

% 

^3 

100 

-348 

0.0996 

9580 

0.07 

0.0004 

2.5 

0.17 

cavitation 

80 

-208 

0.0959 

7380 

0.0006 

6.3 

0.17 

no  cavitation 

60 

- 68 

0.0791 

4560 

0.0004 

23.1 

0.19 

no  cavitation 

40 

73 

0.0142 

550 

39.68 

0.2099 

85.3 

7.45 

cavitation 

20 

213 

0.0060 

120 

62.62 

0.3313 

94.0 

9.94 

no  cavitation 

Seal  D2 

<^1 

= 740  X 

10~^ 

hi  = 

0.0 

100 

-338 

0.1002 

9620 

0.062 

0.0003 

2.4 

80 

-198 

0.0964 

7400 

0.095 

0.0005 

6.3 

60 

- 58 

0.0792 

4560 

0.064 

0.0003 

25.4 

40 

82 

0.0005 

20 

0.525 

0.0028 

100.0 

20 

222 

0.0004 



10 

5.316 

1 

0.0281 

100.0 
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Again,  predicted  leakage  is  larger  than  experimental  leakage  because 
of  probable  wear  on  the  experimental  seal. 

Seal  A2  in  table  4-2  has  about  the  same  initial  taper  as  seal  D2, 
except  that  some  initial  waviness  is  present.  The  transition  between 
high  and  low  friction  takes  place  at  about  the  same  pressure.  However, 
the  transition  for  seal  A2  is  more  gradual.  Also,  the  friction  for  the 
^ 20  and  40  percent  pressure  cases  is  higher  for  seal  A2.  The  reason  is 

that  in  the  converging  cases , the  initial  waviness  causes  some  contact 
1 to  occur  (compare  the  percent  load  support) . Leakage  is  also  higher  for 

{ the  wavy  seal.  Again,  the  reason  is  that  the  waves  cause  the  equilibrium 

! mean  film  thickness  to  be  larger,  causing  greater  leakage.  Leakage  for 

1 seal  A2  is  compared  in  figure  4-4. 

( 

i Although  hydrodyncunic  effects  do  occur  in  seal  A2  (cases  where 

f 

^ cavitation  is  noted) , these  effects  are  much  smaller  than  the  radial 

‘ taper  effects.  It  appears  that  when  the  seal  is  convergent,  the  load 

is  primarily  supported  hydrostatically  and  the  large  amplitude  waves 
touch,  giving  rise  to  greater  friction.  In  the  divergent  cases,  the  net 
waviness  becomes  very  small;  i.e.,  the  seal  flattens  out.  Net  waviness 
is  so  small  that  very  little  hydrodynamic  load  support  would  be  expected. 

As  for  additional  support  to  verify  the  seal  model,  the  general 
friction  and  leakage  behavior  of  the  test  seals  has  been  predicted  by 
the  model.  These  effects  are  primarily  due  to  varying  radial  taper  and 
the  results  support  the  concept  of  load  sharing  between  the  fluid 
pressure  and  asperity  pressure.  As  to  hydrodynamic  effects,  no  support 
to  the  model  is  provided  because  the  hydrodyneunic  effects  in  the  test 
seals  are  overshadowed  by  radial  taper  effects. 

However,  the  concept  herein  of  how  elastic  deflection  in  conjunction 
with  waviness  operates  in  a seal  is  well  supported  by  the  seal  A2  data 
of  figure  4-4.  As  the  seal  becomes  divergent,  it  flattens  considersibly , 
reducing  the  leakage  (see  also  table  4-2) . 

i 
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CHAPTER  5 


ROLE  OF  WAVINESS  IN  SEAL  PERFORMANCE 

The  information  developed  in  chapters  2 and  3 is  purely  theoretical. 
To  determine  the  role  of  waviness  in  non- laboratory  or  as  meuiufactured 
seals,  the  theoretical  results  must  be  interpreted  in  relation  to  cer- 
tain application  parameters.  From  the  previous  results  it  is  known  that 
expected  performcince  depends  to  a large  extent  on  the  actual  values  of 
net  waviness  and  surface  roughness.  So,  in  order  to  determine  if  actual 
seals  operate  hydrodynamically,  representative  values  for  these  param- 
eters must  first  be  determined. 


INITIAL  WAVINESS 

From  previous  work  [32],  [33],  [1]  it  is  known  that  an  initial 
waviness  placed  on  a seal  will  flatten  out  to  a large  extent  as  the  seal 
faces  are  pressed  together.  Thus,  the  net  waviness,  which  is  the  es- 
sential parameter  for  hydrodynamic  effects,  will  always  be  lower  than 
the  initial  waviness.  The  relationship  between  the  two  depends  largely 
upon  the  stiffness  of  the  seal  rings. 


To  quantitatively  assess  this  relationship,  assume  that  the  tan- 
gential varying  pressure  distribution  for  some  particular  set  of  param- 
eters is  described  by  equation  ( 2-62) . Further,  assume  that  only  the 
nth  harmonic  deflection  is  significant;  all  higher  harmonic  deflections 
are  insignificant  as  discussed  in  reference  [1].  The  amount  of  elastic 
deflection  obtained  for  a given  operating  condition  is  then  given  by 
equations  (2-63  ) and  ( 2-64)  where  j = 1. 


2 

n + A 


da. 


K n 


2n^  + T? 


(5-1) 


2 

n + A 


db. 


-6^4  2 

K n - 2n  + n 


(5-2) 
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Define  net  waviness  by  two  terms 


= h.  + h , 
a da^^ 


^3  = ^ib^  ^ ’’db^ 


(5-3) 

(5-4) 


such  that 


(5-5) 


In  all  parameter  studies  made,  h^^  = 0 so  that 


h = h,  = h.  + h,, 
n a la^^  da^^ 


(5-6) 


Solving  for  h.  and  h.^ 
xa^  xbi 


h.  = h , - h 

la^  da^^  n 


h.,  = -h,, 

xb^  db^ 


(5-7) 

(5-8) 


where  h^  and  h^^  are  given  by  equations  (5-1)  and  (5-2)  . Then 


h . = Jh^  + h^ 

X J xaj^  xb^ 


(5-9) 


where  h^  is  the  initial  waviness  required  to  produce  the  net  waviness 


h given  the  elastic  deflection  h,  and  h 
n dai 


dbi 


With  the  above  results  it  is  possible  to  determine  the  initial 
waviness  required  for  any  set  of  operating  conditions.  The  case  2 
parameter  study  of  the  effect  of  net  waviness  (figure  3-5)  has  been  ex- 
tended over  a wider  range  of  net  waviness  and  is  shown  in  figure  5-1. 

The  curves  show  that  the  trends  established  in  figure  3-5  continue  as 
expected  with  percent  fluid  film  load  support  and  lea)cage  becoming 
greater  with  increasing  net  waviness.  Using  the  method  just  discussed, 
figure  5-2  shows  a plot  of  initial  waviness  required  for  the  net  wavi- 
ness curves  of  figure  5-1  for  various  values  of  stiffness  (defined  in 
equation  (2-66).  As  expected  for  larger  stiffness  values,  the  required 
initial  waviness  is  lower.  The  upper  curves  also  show  that,  for  a given 
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initial  waviness  is  lower.  The  upper  curves  also  show  that,  for  a given 
initial  waviness,  there  may  be  three  possible  operating  conditions  (see 
points  A,  B,  and  C on  figure  5-2) . For  the  lower  curves  there  will  be 
only  one  intersection  point. 

To  understand  the  physical  significance  of  the  multiple  intersec- 
tion points  it  is  necessary  to  consider  the  stability  of  the  seal. 
Consider  the  seal  as  a spring  and  assume  that  some  small  additional  ex- 
ternal load  is  applied  to  the  seal  ring.  If  this  additional  load  causes 
the  mean  film  thickness  to  decrease  in  order  to  establish  a new  equi- 
librium position  (displacement  is  in  the  same  direction  as  the  added 
load) , the  seal  will  be  stable.  On  the  other  hand,  if  the  additional 
load  requires  that  the  mean  film  thickness  increase  in  order  to  estab- 
lish the  new  equilibrium  position  (displacement  opposite  added  load) , 
then  the  seal  will  be  unstable  and  operation  at  this  point  will  be 
physically  impossible. 

To  determine  which  of  the  intersection  points  in  figure  5-2  are 

stable,  it  is  necessary  to  evaluate  the  load  support  versus  h^  for  a 

given  and  constant  initial  waviness  h^.  To  do  this  requires  convergence 

on  the  two  deflection  parameters  h,  and  h,,  . That  is,  for  elastic 

dai  dbj^- 

deflection  equilibrium  (given  a particular  h , h.  with  h.,  = 0), 

o 

h , and  h,.  must  be  adjusted  to  satisfy  equations  ( 5-1  ) and  ( 5-2)  • 

CiO 

This  solution  is  obtained  numerically  using  a secant  method  (similar  to 

that  in  reference  [1]  for  three  variables  h , h,  , h..  ). 

o daj^  dbj^ 

Because  of  the  large  amount  of  time  required  for  this  solution,  the 
second  approximation  discussed  in  chapter  2 was  used.  To  determine  if 
stability  results  based  on  the  second  approximation  would  be  valid,  a 
parameter  study,  such  as  shown  in  figures  5-1  and  5-2  which  is  based  on 
the  third  approximation  (the  most  exacting  method  developed) , was  run 
using  the  second  approximation  and  is  shown  in  figures  5-3  and  5-4. 
Comparing  the  results  in  figure  5-3  to  those  in  figure  5-1,  it  is  clear 
that  the  leakage  is  in  considerable  error,  particularly  at  high  net 
waviness  values.  However,  percent  load  support  and  friction  are  in 
reasonable  agreement.  Comparing  the  required  initial  waviness  of  fig- 
ure 5-4  to  figure  5-2  shows  that,  although  the  curves  are  shifted  ver- 
tically, the  general  shape  of  the  curve  is  the  same.  Based  on  this  ob- 
servation, it  was  concluded  that  a stability  evaluation  based  on  the 
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Figure  5-4.  Approximate  Required  Initial  Waviness  - 
Isotropic  Roughness  Case  2 


105 


second  approximation  would  be  valid  at  least  for  general  conclusions. 

Figure  5-5  shows  the  results  of  the  stability  study  based  on  the 
second  approximation.  The  load  support  developed  by  the  fluid  and  as- 
perity pressures  has  been  found  as  a function  of  h^  given  that  the 
elastic 
stant  h. 

iol 

the  = 0 axis.  The  ordinate  is  the  percent  load  support  above  and 
below  the  value  required  for  equilibrium.  A negative  slope  at  the  in- 
tersection point  indicates  stability  (increasing  load  support  with 

decreasing  h ) . 

o 

Now  considering  figure  5-4  in  conjunction  with  figure  5-5  and  the 
case  for  h^  = 65,  Figure  5-5  shows  that  the  first  intersection  is  stable. 
The  intersection  at  the  higher  h^  value  is  unstable.  Considering  the 
case  for  h^  = 17.5,  the  curve  in  figure  5-5  shows  that  an  intersection 
would  be  expected  at  about  h^  = 5.0  which  corresponds  to  the  intersec- 
tion value  in  figure  5-4.  This  root  is  also  stable.  For  = 15  the 

intersection  point  is  a stable  root.  For  h^  = 6.0  the  intersection 
point  is  also  a stable  root. 

Looking  at  figure  5-4  in  conjunction  with  figure  5-5  shows  that  all 
intersection  points  in  figure  5-4  where  the  slope  is  positive  represent 
stable  solutions.  It  can  be  concluded  that  stable  operation  can  be  ex- 
pected to  the  left  of  the  peaks  of  the  upper  curves  in  figure  5-4.  For 
the  lowest  two  curves,  any  intersection  point  represents  a stable  opera- 
ting condition. 

In  actuality  there  will  always  be  three  intersections  for  the  upper 
curves  in  figure  5-4.  The  third  intersection  point  is  beyond  the  range 
of  the  curves.  If  a seal  corresponding  to  the  upper  curve  in  figure  5-4 
and  having  a variable  waviness  were  operated,  as  initial  waviness  was 
increased  from  h.  = 0 to  say,  h.  =65  where  h.  = 0.5,  friction  would  de- 
crease  and  leakage  would  increase  as  the  percent  fluid  film  load  support 
increased.  That  is,  the  seal  would  become  partially  hydrodynamically 
supported.  If  the  initial  waviness  were  increased  beyond  the  top  of  the 
curve,  then  the  net  waviness  would  suddenly  jump  beyond  the  range  of  the 
curve.  Leakage  would  drastically  increase  and  operation  would  become 
stable  at  this  higher  leakage  rate. 


displacements  h^^  and  h^  are  in  equilibrium  and  given  a con- 
and  K.  An  equilibrium  solution  occurs  when  the  curves  cross 
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Figure  5-5.  Approximate  Load  Function  - Isotropic  Roughness  Case  2 


For  the  lowest  two  curves,  no  sudden  change  in  operation  would  oc- 
cur with  increasing  initial  waviness.  However,  the  slope  of  the  curve 


does  change  from  being  very  steep  for  low  values  of  h^  to  being  very 
shallow  for  higher  values  of  h^.  Thus,  the  point  of  operation  on  the 
curve  becomes  very  sensitive  to  the  value  of  the  initial  waviness  after 
the  low  slope  values  are  reached. 

In  order  to  determine  the  effect  of  surface  roughness  on  the  stable 
range  of  operation,  a net  waviness  study  was  made  for  a roughness  value 
one  half  of  that  of  base  case  2.  The  performance  curves  are  shown  in 
figure  5-6.  Corresponding  required  initial  waviness  curves  are  shown  in 
figure  5-7'.  As  expected,  figure  5-6  shows  that  the  percent  of  fluid 
film  load  support  increases  to  100  percent  at  a very  small  waviness 
value.  Friction  decreases  and  leeUcage  increases  with  increasing  wavi- 
ness. The  required  initial  waviness  curves  are  similar  in  shape  to  those 
in  figure  5-4.  For  the  upper  curves  it  is  possible  to  achieve  nearly 
100  percent  hydrodynamic  load  support  before  the  peak  of  the  curve  and 
the  unstable  point  of  operation  is  reached.  The  slope  of  the  lower 
curves  is  very  low  indicating  that  the  operating  condition  could  change 
radically  with  a small  change  in  initial  waviness. 

In  conclusion,  it  appears  that  only  a limited  range  of  initial  wavi- 
ness values  will  lead  to  an  operating  condition  where  the  fraction  of 
fluid  film  load  support  is  significantly  above  the  hydrostatic  value 
while  at  the  same  time  allowing  only  a small  leakage  to  occur.  The  re- 
quired initial  waviness  to  cause  operation  to  occur  near  one  of  the  peaks 
of  figure  5-4  depends  on  the  stiffness  of  the  ring,  a low  stiffness  re- 
quires a very  large  initial  waviness.  This  question  will  be  examined  in 
relation  to  some  applications  later  in  this  chapter. 

SURFACE  ROUGHNESS 

The  results  of  chapter  3 show  that  surface  roughness  is  one  of  the 
most  important  parameters  in  determining  the  extent  of  hydrodynamic  lubri- 
cation in  wavy  seals.  Keeping  in  mind  that  the  roughness  representation 
used  in  the  model  is  an  idealization  of  actual  roughness,  some  typical 
roughness  values  can  be  compared  to  those  used  for  the  studies  of 
chapter  3. 


Table  5-1  shows  some  typical  roughness  values  for  carbon  seal  faces 
after  a period  of  operation  in  water.  The  extent  of  hydrodynamic  load 
support  due  to  waviness  for  these  seals  is  unknown  but  it  was  probably 
quite  small  because  of  the  high  degree  of  flatness  to  which  the  test 
seals  were  lapped.  However,  these  values  are  assumed  to  be  representa- 
tive of  those  that  might  be  obtained  for  a wavy  seal . Figure  5-8  shows 
the  surface  profile  for  the  first  seal  in  table  5-1.  This  profile  is 
typical  of  carbon  seal  faces.  There  are  a number  of  deep  grooves  or 
craters.  The  top  of  the  profile  appears  to  be  somewhat  flattened.  The 
actual  height  distribution  for  this  profile  would  be  greatly  skewed  such 
as  obtained  by  truncating  a Gaussian  distribution.  Figure  5-9  shows  a 
surface  profile  in  a tangential  direction.  The  profile  is  similar  to 
that  in  the  radial  direction  and  this  supports  the  assumption  of  iso- 
tropic roughness  used  for  most  of  the  parameter  studies.  Figure  5-10 
shows  a radial  surface  profile  for  a seal  face  with  a large  roughness. 
There  appear  to  be  definite  deep  radial  grooves  in  this  face.  A tan- 
gential trace  for  the  same  seal  (not  shown)  appears  similar  to  the  ra- 
dial trace  without  the  grooving  effect. 

Table  5-1  shows  two  additional  results  obtained  from  the  litera- 
ture. The  values  reported  by  Pape  [30]  for  a lapped  seal  surface  are 
taken  as  representative  of  what  can  be  expected  for  initial  surface 
roughness.  Kojabashian  and  Richardson  [16]  did  not  actually  report  the 
CLA  value  but  determined  the  average  height  of  microasperities  in  the 
tangential  direction.  Interpreting  this  average  value  in  terms  of  CLA 
is  approximate. 

Considering  the  values  reported  in  table  5-1  in  relation  to  the 
values  of  c = 1.0  cuid  c = 0.5  ym  for  cases  1 and  2 respectively,  it  ap- 
pears that  test  seal  14  could  clearly  operate  in  the  case  1 range.  The 
Kojabashian  and  Richardson  seal  might  operate  somewhere  between  case  1 
and  case  2.  The  as  lapped  seal  could  operate  at  a roughness  even  less 
than  case  2,  indicating  that  100  percent  fluid  film  could  easily  be 
achieved  at  least  until  the  surface  roughness  became  greater. 

What  is  important  to  observe  is  that  based  on  these  very  approxi- 
mate roughness  comparisons  and  the  approximate  roughness  model  used, 
actual  operating  seal  faces  (test  seal  #14  in  particular)  do  have  rough- 
ness values  that  are  low  enough  such  that  hydrodynamic  effects  might  be 
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Estimate  based  on  assuming  that  roughness  conforms  to  the  polynomial 
distribution  of  equation  ( 2-5  ) • 

’ it 

c is  taken  as  one  half  of  the  mean  value  of  6 obtained  for  method  A — 
as  run  seals. 
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Figure  5-8.  Radial  Surface  Profile  of  Seal  #14 


expected  to  occur,  given  sufficient  waviness.  It  is  also  possible, 
given  the  inaccuracy  of  the  model  and  the  approximate  representation 
of  the  surface,  that  only  in  the  smoothest  seal  faces  will  hydrodynamic 
effects  be  significant. 

In  order  to  predict  the  actual  extent  of  hydrodynamic  load  support 
for  a given  seal  with  a given  roughness,  considerable  refinement  in  the 
roughness  model  will  be  required.  Some  of  these  deficiencies  are  dis- 
cussed at  length  in  chapter  7. 

The  fact  that  hydrodynamic  effects  are  very  sensitive  to  surface 
roughness  values  (chapter  3)  in  conjunction  with  the  fact  that  surface 
roughness  values  can  vary  over  a very  wide  range  (table  5-1)  may  explain 
some  of  the  large  variations  found  in  friction  values  and  wear  rates 
in  tests  of  face  seals.  Although  seal  faces  are  lapped  to  very  low 
roughness  values  initially,  little  attention  has  been  given  to  the 
roughness  that  results  after  a period  of  operation. 

For  the  remainder  of  this  work  it  will  be  assumed  that  the  surface 
roughness  is  low  enough  or  could  be  made  low  enough  by  design  so  that 
significant  hydrodynamic  efforts  might  occur.  Given  this  assumption, 
further  studies  can  be  made. 

EVALUATION  OF  MANUFACTURED  SEALS 

Given  the  assumption  in  the  previous  section  that  an  effective 
roughness  of  c = 0.5  ym  (20  yin)  can  be  achieved  or  actually  occurs  in 
some  operating  seals,  then  one  need  consider  the  typical  magnitude  of 
initial  waviness  values  found  in  operating  seals  to  establish  the  con- 
ditions under  which  hydrodynamic  effects  might  be  occurring. 

To  serve  as  a first  example  the  3 5/8  inch  seals  of  reference  [321 
have  been  selected  for  study.  All  case  2 parameter  values  of  equation 
(3-3)  apply  to  this  particular  seal.  To  obtain  the  stiffness  K,  ref- 
erence [32]  gives 

E = 2.06  X 10^°  P (3  X 10^  psi) 

d 

J = 4.83  X lo"^  m'*  (0.0116  in'^)  (5-10) 

X 

A = 6 
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Using  values  in  equation  (2-108)  with  c = 0.5  pm,  K = 2.77  x lo  ^ which 

-7 

will  be  taken  as  K = 2.5  x lo  for  evaluation  purposes.  Figure  5-2 
shows  that  for  h^  = 0.2  (the  beginning  of  hydrodynamic  effects — see 
figure  5-1)  an  initial  waviness  of  h^  = 31  is  needed.  For  c = 0.5  pm, 
then  h^  = 15.5  pm  (610  pin)  is  the  value  needed.  Reference  [32]  provides 
some  experimental  results  for  the  initial  waviness  developed  in  these 
particular  carbon  rings  due  to  drive  forces.  These  data  are  presented 
in  Table  5-2  where  it  has  been  assumed  the  three  standard  deviations 
establish  the  maximum  possible  value  of  waviness. 


Table  5-2 

Experimental  Waviness  - 3-5/8"  Seal  [32] 


Clearly  the  maximum  initial  waviness  available  is  much  less  than  the 
15.5  Pm  need  to  achieve  the  base  case  2 conditions.  This  particular 
seal  will  thus  have  no  hydrodynamic  effects  at  the  base  case  condition — 
at  least  not  due  to  third  harmonic  waviness. 

Now  to  detemine  under  what  conditions  this  particular  seal  might 
have  significant  hydrodynamic  effects,  the  case  2 pressure  study  results 
were  used  (figure  3-9).  The  required  initial  waviness  for  these  results 
is  shown  plotted  in  figure  5-11.  Now  taking  the  maximum  third  harmonic 
waviness  from  table  5-2  of  h^  = 1.36  (h^  = 2.72),  figure  5-11  shows  that 
the  seal  would  operate  with  a net  waviness  of  h^  = 0.2  at  p^  = 0.00025 
which  corresponds  to  0.14  MPa  (21  psi)  for  the  base  case  values.  Fig- 
ure 3-9  shows  that  the  seal  would  be  expected  to  have  a very  large 
hydrodynamic  load  support.  Of  course,  hydrodynamic  effects  would  per- 
sist in  this  seal  even  at  higher  pressures.  The  net  waviness  would 
simply  be  lower  than  the  0.2  for  which  the  results  were  tabulated.  At 
p^  = 0.005  the  net  waviness  would  be  so  low  that  no  significant  hydro- 
dynamic  effects  would  occur. 
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Figure  5-11.  Required  Initial  Waviness  - Isotropic  Roughness  Case  2 


Now  suppose  this  same  seal  were  ten  times  stiffer  than  for  the 
given  data.  This  might  occur  if  the  carbon  were  mounted  in  a steel 
backup  ring  or  if  the  cross  sectional  dimensions  were  increased.  Then, 
for  the  same  h^  = 2.72,  figure  5-11  shows  that  h^  = 0.2  would  occur  when 
p^  = 0.0036  which  corresponds  to  p^  = 2MPa  (300  psi) . Figure  3-9  shows 
that  about  73  percent  of  the  load  would  be  supported  by  fluid  film  pres- 
sure. At  lower  pressures,  with  the  same  initial  waviness,  the  net  wavi- 
ness would  increase  above  0.2  as.d  hydrodynamic  load  support  would  increase 
significantly. 

In  conclusion,  whether  or  not  a seal  will  benefit  (or  performance 
become  unacceptable  because  of  high  leakage)  from  significant  hydro- 
dynamic  effects  depends  largely  on  the  stiffness  and  initial  waviness. 

The  results  show  that  such  effects  would  occur  possibly  at  only  very  low 
pressures  for  the  plain  carbon  seal  ring  used  in  the  example.  For  stif- 
fer rings,  then  hydrodynamic  effects  may  play  a significant  role  even  for 
small  amounts  of  initial  waviness.  It  is  impossible  to  make  any  gen- 
eral statements  beyond  this.  However,  using  the  methods  developed,  a 
prediction  of  hydrodynamics  can  easily  be  made  for  a given  seal. 

Although  initial  waviness  has  been  measured,  the  fact  of  ring  stiff- 
ness has  been  given  little  or  no  attention  in  experimental  seal  research 
in  the  past.  These  results  show  that  a variation  of  stiffness  only,  all 
other  factors  being  the  same,  could  cause  operation  to  shift  from  the 
purely  hydrostatic  to  complete  fluid  film  load  support.  These  results 
may  also  help  explain  why  such  large  differences  in  seal  wear  rate  often 
occur  in  tests  where  all  factors  are  otherwise  the  same.  Differences  in 
initial  waviness  can  lead  to  large  differences  in  seal  performance  under 
some  operating  conditions. 


EFFECT  OF  WEAR 

The  preceding  results  apply  to  waviness  profiles  that  are  a pure 
nth  harmonic  sine  wave.  Wear  will  modify  this  shape.  To  understand 
what  happens  to  the  performance  of  an  initially  hydrodynamic  seal  after 
a period  of  operation,  the  effects  of  wear  on  performance  must  be 
determined. 

Following  the  method  outlined  in  reference  [1] , it  is  assumed  that 


r 


k 

[ 

I 

I 

E. 


i 


i 

f 

\ 

t 

f 

[ 

i' 

I 


wear  rate  is  given  by  the  following  expression* 


dh 

= C p (1  - b )V 
dt  o m h mecm 


(5-11) 


where  C is  a wear  constant,  p (1  - b,  ) is  the  average  mechanical  pres- 
o m h 

sure  and  V is  the  average  speed  at  the  face  of  the  seal, 
mean 

In  dimensionless  form  the  above  expression  becomes 


dh 
w 

dt 


m 


(1  + r^) 
2 


(5-12) 


where 


t = 


^3  2^ 

tr  nw  c 
o o 

3 

c 


(5-13) 


and  other  symbols  are  as  defined  previously. 
Noting  that 


+ w. 


avg 


(5-14) 


where  p (1  - b,  ) is  the  average  asperity  load  around  the  seal,  and  W 
n*  h avg  ^ 

is  the  total  load  supported  by  fluid  pressure,  then  it  can  be  shown  using 

equations  (5-12)  and  (5-14)  that 


“^w  (1  - H.)W 

o _ y 

dt  2ii(l  - r^) 


(5-15) 


where 

dh- 

w 

o 

dt 

is  the  average  wear  rate  and  is  the  fraction  of  load  supported  by 
fluid  pressure.  Then  the  average  wear  rate  is  reduced  as  the  fraction 
of  fluid  film  pressure  load  support  increases. 

Wear  h is  represented  by  a Fourier  series 
w 
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h = h + E h cos  n j0  + h , sin  n j6 
w w wa . wb . 

o D J 


(5-16) 


The  change  in  wear  over  time  At  is  thus 

Ah  = Ah  + EAh  cos  nj6  + Ah  , sin  nj0 
w w wa . wb . 

o D D 


(1  + r ) _ 

= p (1  - b,  ) r-^  At 

m h 2 


(5-17) 


The  coefficients  h and  h ^ are  found  by  making  a Fourier  analysis  of 
wa  wb 

p (1  - b,  ) . These  increments  are  then  added  to  running  totals  of  h 
m h 

and  h , . 

wb  j 


To  solve  the  time  dependent  problem,  first  the  equilibrium  problem 
for  a given  initial  waviness  is  solved  at  t = 0 (this  requires  the  simul- 
taneous solution  of  three  functions  in  three  unknowns  h , h,  and 

o daj^ 

h„  to  satisfy  load  support  and  deflection  equilibrium  equations 
dbj 

( 5-1  ) and  ( 5-2  ) ) . Then  a Fourier  analysis  is  made  of  the  mechanical 
pressure  distributions,  and  the  wear  increments  are  calculated  using 
equation  (5-17).  These  increments  are  added  to  the  initial  waviness, 
and  the  entire  process  is  repeated  for  each  time  increment. 

The  simultaneous  solution  for  the  three  varieibles  above  requires  a 
considerable  amount  of  computer  time.  Since  this  solution  must  be  ob- 
tained for  each  time  step,  enormous  amounts  of  time  can  be  consumed  even 
for  a short  time  dependent  study.  To  overcome  this  limitation  to  some 
extent,  the  second  approximation  discussed  in  chapter  2 was  used  rather 
than  the  more  time  consuming  third  approximation. 

Base  case  2 parameters  were  used  for  a time  dependent  study.  As  a 
starting  point,  h^  = 0.5  from  figure  5-3  was  chosen.  For  K = 5 x 10  ^ 
figure  5-4  shows  that  h^  = 16.5  is  the  required  initial  waviness.  The 
time  dependent  results  are  shown  in  figure  5-12.  Load  support  increases 
with  time  initially  as  wear  progresses.  The  explanation  for  this  be- 
havior is  that  initially  wear  induces  favorable  changes  in  the  pressure 
distribution  which  in  turn  reduce  the  elastic  deflection.  The  net  wavi- 
ness actually  increases.  However,  after  a long  period  of  time,  fluid 
pressure  load  support  begins  to  decrease.  Leakage  continues  to  increase 
with  time  because  h increases  with  time,  h increases  because  net 
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Figure  5-12.  Time  Dependent  Behavior  - Isotropic  Roughness  Case 


waviness  increases  as  described  above. 

In  terms  of  real  time,  if  we  assume  a typical  seal  material  that 
wears  0.01  mm  per  100  hours  under  base  case  2 conditions,  assuming  zero 
fluid  film  load  support,  then  100  hours  corresponds  to  t = 3350.  So, 
the  entire  curve  in  figure  5-12  corresponds  to  about  100  hours  of 
operation . 

What  is  expected  to  occur  after  a long  period  of  operation  (beyond 
the  range  of  figure  5-12)  is  that  wear  will  continue  and  the  fluid 
pressure  load  support  will  decrease  until  some  point  at  which  hydro- 
dynamic  load  support  no  longer  acts.  More  studies  are  needed  to  actually 
verify  this  behavior. 

One  of  the  limitations  on  the  accuracy  of  this  method  is  pointed 
out  by  the  results  shown  in  figure  5-13.  The  pressure  distributions, 
p(0) , are  smooth  at  t = 0.  However,  after  a period  of  wear,  t = 3370, 
both  the  hydrodynamic  and  the  mechanical  pressure  distributions  pick  up 
a significant  ripple.  The  harmonic  number  of  this  ripple  appears  to 
correspond  to  one  plus  the  highest  harmonic  term  used  in  equation  (5-17) 
to  calculate  the  wear  increment.  It  is  thought  that  this  ripple  may  be 
leading  to  some  error  in  the  results  and  should  be  corrected  for  further 
studies.  The  reason  for  ripple  is  that  the  series  has  been  truncated 
at  j = 10.  Expanding  the  number  of  terms  in  the  series  would  alleviate 
this  problem  but  this  cannot  be  done  because  too  many  points  in  the  1 
direction  are  required  to  properly  define  these  higher  order  harmonics. 
However,  a method  is  now  under  development  which  eliminates  the  use  of 
the  Fourier  series.  Wear  is  calculated  separately  for  each  point  in 
the  mesh. 

In  conclusion,  the  results  show  that  wear  causes  favorable  changes 
in  seal  performance  in  the  short  run.  However,  the  results  suggest  that 
wear  will  eventually  cause  hydrodynamic  effects  to  vanish. 
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CHAPTER  6 


WAVY  SEAL  DESIGN 

The  results  in  the  previous  chapters  suggest  that  a considerable 
reduction  in  wear  rate  could  be  achieved  by  designing  waviness  into  a 
seal.  For  example,  consider  a 75  percent  balanced  seal  where  the  hydro- 
static load  support  for  parallel  faces  is  69  percent.  Now  suppose  the 
fluid  film  load  support  were  increased  to  85  percent  by  hydrodynamic 
means.  According  to  equation  (5-15) , average  wear  rate  would  be  re- 
duced a factor  of  2-1  con^ared  to  the  hydrostatic  case.  If  90  percent 
load  support  were  achieved,  wear  rate  would  be  reduced  a factor  of  3-1. 

Even  though  considerable  gains  ccin  be  achieved,  the  results  of 
chapter  5 show  that  the  effect  of  any  initially  planned  waviness  will 
decrease  with  time.  However,  if  the  wear  can  be  distributed  evenly 
by  moving  the  waviness  relative  to  the  seal  face,  the  effect  of  initial 
waviness  can  be  prolonged  indefinitely.  Reference  [1]  describes  the 
principles  of  doing  this. 

Now  given  the  possibility  of  being  able  to  utilise  initial  waviness 
and  hydrodyneimic  effects  to  reduce  wear  rate  for  the  entire  life  of  the 
seal,  there  are  several  design  factors  that  must  be  considered. 

SELECTION  OF  VARIOUS  PARAMETER  VALUES 

The  results  of  chapter  3 show  clearly  that  load  support  increases 
with  decreasing  roughness  (figure  3-13) . Leakage  also  decreases  with 
decreasing  roughness.  For  best  performance,  it  is  necessary  that  ma- 
terials be  selected  that  provide  the  lowest  possible  as-running  surface 
roughness.  The  results  in  table  5-1  suggest  that  there  may  be  wide 
variations  in  roughness  for  different  material  pairs.  Further  research 
is  needed  to  determine  what  material  factors  or  parameters  are  related 
to  as-running  roughness  and  to  establish  particular  material  pairs  that 
give  low  roughness.  For  the  purposes  of  the  discussion  in  this  chapter, 
it  will  be  assumed  that  a roughness  value  as  low  as  c = 0.5  ym  is 
practically  achievable. 

Figure  3-10  shows  that  the  best  performance  can  be  achieved  by 
increasing  face  width.  Percent  load  support  increases  and  leakage 
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decreases  with  increasing  face  width.  It  has  been  a practice  in  the 
seal  industry  to  use  a narrow  face  on  seals  in  the  order  of  r^  = 0.9. 

The  reason  for  this  is  that  the  heat  generated  at  the  seal  interface 
is  roughly  proportional  to  face  area.  Thermal  coning  or  the  rotation 
of  the  face  due  to  thermal  expansion  (causing  a radially  converging 
taper)  is  in  turn  proportional  to  the  total  heat  generated.  Thus,  a 
larger  radial  taper  will  result  for  a greater  face  width.  This  may 
lead  to  high  leakage  and  uneven  radial  wear.  Also,  the  increased 
operating  temperature  at  the  seal  face  can  lead  to  problems  resulting 
from  vaporization  of  the  sealed  liquid  at  the  interface.  At  present, 
it  appears  that  it  is  best  to  keep  the  face  width  narrow,  in  spite  of 
the  theoretical  hydrodyneunic  advantages  of  a wider  face. 

Figure  3-12  shows  that  improved  performance  can  be  obtained  by 
increasing  the  number  of  waves.  Load  support  increases  while  leakage 
stays  about  constant.  Barring  any  difficulties  in  producing  a large 
number  of  waves  (presently  under  investigation) , it  appears  that  n 
should  be  made  large  (the  limits  of  this  behavior  are  being  studied) . 

Now  limiting  the  discussion  to  a case  where  c=0.5  ym,  n=3,  and 
= 0.9  (until  further  information  is  available  on  the  practical  limits 
of  c and  n)  and  the  base  case  2 operating  conditions,  figure  5-1  shows 
the  full  range  of  performance  that  can  be  expected  using  a hydrodynamic 
seal.  The  remaining  variable  over  which  one  has  control  is  the  net 
waviness  h^.  Table  6-1  shows  some  of  the  various  operating  conditions 
extracted  from  figure  5-1.  The  table  shows  that  significant  reductions 
in  wear  rate  are  possible.  The  amount  depends  on  the  acceptable  level 
of  leakage.  For  many  applications  the  larger  rates  shown  in  the  table 
would  be  acceptable. 

Now  to  see  what  is  practically  achievable  from  the  standpoint  of 

stability,  figure  5-2  shows  that  conditions  1 eind  2 (table  6-1)  can  be 

reached  using  any  of  the  stiffness  values  shown.  To  reach  conditions  3, 

4,  5 and  6 would  require  the  use  of  the  largest  stiffness  K = 2.0  10 

Now  if  initial  waviness  could  be  controlled  to  5 units  ±1,  with 
— —6  — 

K = 2.0  X 10  , then  h^  would  range  from  0.2  to  3.0.  Table  6-1  shows 

that  this  would  cause  performance  to  vary  over  a very  wide  range.  Either 
more  precise  control  of  initial  waviness  would  be  required  or  perhaps 
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Table  6-1 


Various  Possible  Operating  Conditions 
for  a Hydrodynamic  Seal 


Condition 

^3 

% 

Q 

Q cc/min 

— f 1 1 11  1 LI.  M 

Hydrostatic  Wear  Rate 

1 

m 

72.5 

0.031 

0.003 

0.04 

0.89 

2 

0.5 

82.5 

0.020 

0.007 

0.10 

0.56 

3 

0.75 

87.0 

0.015 

0.013 

0.18 

0.42 

4 

1.0 

90.0 

0.012 

0.021 

0.30 

0.32 

5 

1.5 

93.0 

0.009 

0.048 

0.68 

0.23 

6 

2.0 

95.0 

0.007 

0.096 

1.35 

0.16 

even  larger  stiffness  values  than  shown  in  figure  5-2  (giving  a lower 
sensitivity)  can  be  used  (this  is  still  under  investigation) . 

In  summary,  the  results  show  that  hydrodyneunic  water  seals  could 
be  designed  which  would  have  very  low  wear  rates  and  acceptably  low 
leakage.  Performance  would  be  expected  to  be  consistent  and  pre- 
dictable over  the  entire  life  of  the  seal.  Several  questions  still 
remain  to  be  answered  as  pointed  out.  One  additional  question  that 
requires  further  study  is  how  such  a seal  will  behave  in  a real  system 
where  temperature  and  pressure  changes  cause  variations  in  radial  taper 
as  a function  of  operating  conditions. 


WAVE  SHAPE 

All  of  the  previous  studies  herein  were  made  using  a simple  cosine 
wave  shape.  However,  this  may  not  be  the  optimum  wave  form.  A dif- 
ferent wave  shape  might  improve  hydrodynamic  load  support  while  main- 
taining a given  acceptable  leakage.  Such  profiles  would  be  termed 
optimum.  This  was  shown  to  be  true  using  the  one-dimensional  model  of 
reference  [1],  although  leakage  was  not  considered. 

To  answer  this  question  for  the  present  isotropic  model,  a gen- 
eral shape  described  by  a two  term  Fourier  series  was  used. 

h = h cos  n0  + h cos  2n6  + h.  sin  2n6  (6-1) 

1 ®2  “2 
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is  always  set  equal  to  0.0  since  this  term  merely  positions  the 
wave.  The  use  of  three  terms  allows  the  profile  some  amount  of  freedom 
to  assume  different  shapes  within  a period,  while  at  the  same  time  the 
use  of  only  three  terms  permits  a solution  to  be  carried  out  within  a 
reasonable  amount  of  computer  time.  Some  preliminary  work  was  carried 
out  to  see  what  kind  of  leakage  and  load  support  behavior  occurs  using 
the  Fourier  terms  mentioned  above.  Figures  6-1  and  6-2  show  the  be- 
havior of  the  hydrodynamic  load  support  and  leakage  functions  respec- 
tively for  a constant  value  of  h =0.2  and  various  values  of  h 

®1  ®2 

cuid  h,  . A comparison  of  the  two  curves  shows  that  a Q = constant  curve 
°2 

can  be  followed  and  some  small  increases  in  percent  load  support  can  be 
achieved. 

A three  variable  optimization  program  was  used  to  determine  the 

optimum  wave  shape  as  described  by  h , h , and  h,  . The  optimization 

®1  “2 

routine  used  was  a search  method  referred  to  as  the  pattern  search 
method,  combined  with  axis  rotation  and  "successful  direction  remember- 
ing [55]."  Basically,  h and  h,  were  set  equal  to  some  initial 

(starting  point)  value.  Then  the  h required  to  yield  the  given  leak- 

^1 

age  was  found  numerically.  Once  h was  found,  the  pattern  search 

®1 

routine  optimized  on  h and  h.  . Any  time  either  h or  h,  were 

^2  ^2  ^2  ^2 
chcinged,  while  being  optimized,  the  h required  to  achieve  the  desired 

^1 

leakage  was  recalculated. 

Two  different  starting  points  (designated  SPl  and  SP2  in  figure  6-1) 
were  chosen  and  two  different  "answers"  were  found.  The  approximate 
locations  of  these  starting  points  are  shown  in  figure  6-1  as  are  the 
general  directions  the  optimization  routine  proceeds  from  each  starting 
point.  It  can  be  seen  from  figure  6-1  that  the  ultimate  direction  the 
optimization  routine  proceeds  is  a function  of  the  starting  point.  Note 


that  figures  6-1  and  6-2  are  for  a constant  value  of  h =0.2. 


The 


three  variable  optimization  program  varied  h as  well  as  h and  h,  . 

Therefore  the  general  directions  shown  in  figure  6-1  are  actually 

projections  onto  the  h =0.2  plane. 

®1 

Figure  6-3  depicts  a portion  of  the  optimization  procedure  that 
occurred  in  quadrant  IV  (refer  to  figure  6-1) . The  wave  shapes  with  the 
higher  numbers  are  the  successively  improved  shapes,  as  dt‘;ermined  by 
the  optimization  routine.  Each  numbered  wave  shape  in  figure  6-3  is  a 
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Figure  6-3.  Wave  Shape  During  Optimization  (for  Quadrant  IV) 


4 


single  wave  over  one  of  the  n periods  (n  = 3 for  this  solution) . It 
appears  that  the  optimization  routine  is  attempting  to  place  two  waves 
where  there  once  was  one.  That  is,  there  originally  were  three  waves 
around  the  seal  face.  The  optimization  program  is  trying  to  put  six 
waves  around  the  seal  face.  Very  similar  behavior  was  found  when 
optimizing  in  quadrant  III.  Consequently  the  following  analysis 
applies  to  both  quadrants  III  and  IV  even  though  only  the  results  in 
quadrant  IV  are  shown  in  figure  6-3.  Table  6-2  shows  the  results  of 
the  optimization. 

Table  6-2 


Fluid  Film  Load  Support  for  Wave  Shapes  of  Figure  6-3,  Q = const 


Wave  Shape  Number 

Quadrant  III 

(not  shown) 

Quadrcuit  IV 

(Fig.  6-3) 

1 

81.1% 

83.8% 

2 

82.5% 

86.0% 

3 

84.6% 

87.5% 

4 (not  shown) 

85.9% 

For  comparison  to  the  above,  if  the  number  of  waves  around  the  seal 
face  is  changed  from  three  to  six,  and  the  shape  is  changed  back  to  the 
original  cosine  wave  shape,  the  resulting  percent  hydrodynamic  load 
support,  for  the  same  desired  leakage,  is  87.5%.  Therefore,  it  is  ap- 
parent that  the  optimization  routine  is  changing  the  number  of  waves 
from  three  to  six  in  order  to  maximize  load  support  rather  than  changing 
the  wave  shape  itself. 

Based  on  this  limited  study,  it  appears  that  no  improvement  in 
load  support  at  constant  leakage  can  be  achieved  beyond  that  associated 
with  the  simple  cosine  shape. 
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CHAPTER  7 

EVALUATION  OF  ROUGHNESS  MODELS 

In  this  chapter  certain  questions  are  raised  concerning  the  basic 
lubrication  and  load  sharing  models  for  rough  surfaces  used  in  this 
work.  These  questions  are  raised  in  order  to  place  the  user  of  this 
work  on  notice  of  its  limitations,  as  well  as  to  delineate  certain 
areas  of  further  investigation  whose  results  may  significantly  affect 
and  improve  the  type  of  lubrication  model  developed  herein. 

APPLICABILITY  OF  THE  REYNOLDS  EQUATION  TO  ROUGH  SURFACES 

It  is  generally  accepted  that  the  Reynolds  equation  is  applicable 
to  rough  surfaces  if  the  wave  length  of  the  roughness  is  large  compared 
to  the  film  thickness  [50].  For  the  work  herein,  the  film  thickness 
is  of  the  same  order  as  the  roughness  height.  A typical  roughness 
trace  such  as  figure  5-8  shows  that  wave  length  is  many  times  greater 
than  roughness  height  so  that  the  above  condition  is  satisfied.  If 
this  condition  were  not  met,  then  the  Reynolds  equation  would  not  be 
applicable  and  other  methods  of  analysis  such  as  the  use  of  Stokes 
equations  [56]  would  be  required. 

REYNOLDS  EQUATION  FOR  ROUGH  SURFACES 

Godet  [57]  has  derived  the  Reynolds  equation  for  fluids 
between  two  rough  surfaces  from  first  principles.  This  formulation 
differs  from  the  classical  approach  in  that  the  individual  surface 
roughnesses  are  not  lumped  onto  one  surface.  In  general,  the  statisti- 
cal combination  of  two  surface  roughnesses  is  not  easily  accomplished 
nor  is  such  a combination  an  adequate  description  of  the  actual  film 
thickness.  The  elements  of  the  Godet  formulation  will  be  outlined  here 
to  obtain  the  Reynolds  equation  for  rough  surfaces. 

Consider  two  vertically  aligned  surface  segments  Sj  and  at  sur- 
face heights  Hj  and  H2  (figure  7-1) . The  fluid  flow  between  the 
segments  is  governed  by  the  simplified  Navier-Stokes  equations: 


i 

i 

i 

I 


J 


!] 


I 

i 

! 


1 


133 


Figure  7-1.  Geometry  for  Derivation  of  Reynolds 
Equation  for  Rough  Surfaces 


Figure  7-2.  Geometry  of  Rough  Surfaces 


L3' 


I 


u = Uj, 

V = v^, 

w = Wj 

(7-4) 

^2^ 

u = U2, 

CM 

> 

II 

> 

w = w 

2 

(7-5) 

Mclking  the  Reynolds  assiomptions  that  density,  pressure,  and  viscosity 
do  not  vary  in  the  z direction  by  dimensional  considerations,  then 
equations  (7-1)- (7-3)  may  be  integrated  to  solve  for  the  fluid  velocity 
functions  u and  v.  These  functions  may  in  turn  be  integrated  to  yield 


the  leakage  fluxes  q and  q ; 

X y 

“2 

OudZ  - I (H,  * - H,)  - (Hj  - H,) 


= / 


(7-6) 


H 


= /' 


pvdz  = j (Vj  + V^) (H^ 


(7-7) 


H 


The  equation  of  continuity  is  applied  in  the  form: 
H.  H„  H.  H 


r ^ dz  t r 

dz  + f 

J 3t 

J 

3x  J 

3 y J 

H 

H 

H 

H 

1 

1 

1 

1 

After  performing  the 

integrations 

in  equation  (7- 

3 (pw) 
3z 


dz  = 0 (7-8) 


results  (7-6)  and  (7-7) , the  unabbreviated  Reynolds  equation  is  obtained. 
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For  the  sealing  conditions  discussed  herein,  the  density  and  viscosity 
are  constant.  Further,  the  surface  velocities  = V2  = 0,  and  Uj,  U2 
are  constant.  The  Reynolds  equation  becomes; 


_3_ 

3x 


3^ 

3x 


H.)  . 

1 3y 


(Hi  + H^) 

= enCUi-u^)  3— 


+ I2n 


* ^ [ ‘"1  "2'] 


(7-9) 


Because  of  the  appearance  of  the  quantity  (Hj  + H2)  in  (7-9),  Godet  [57] 
has  pointed  out  that  equation  (7-9)  agrees  with  the  usual  Reynolds 
equation  if  Hi  = 0.  This  is  equivalent  to  lumping  the  surface  profile 
Hi  onto  H2  as  was  done  in  this  report,  which  may  be  incorrect  for  seal 
surfaces. 

For  surfaces  having  roughness,  the  total  surface  height  H can  be 
considered  as  being  composed  of  a nominal  height  h and  a roughness 
deviation  6.  Referring  to  figure  7-2: 


H 

1 


H2  = h^  + 


(7-10) 


The  point  velocities  Wi , are  evaluated  with  the  aid  of  figure 
7-3.  It  is  evident  from  this  figure  that: 

3H,  3H. 

^ = “i  17  ^ TT  ' i = 1'  2 (7-11) 


For  the  case  of  interest,  only  the  roughness  61  and  6^  vary  in  time; 
the  nominal  film  thickness  h is  fixed  in  time  and  varies  only  in  space. 
Hence,  3hi^/3t  = 0 and  introduction  of  (7-10),  (7-11)  into  (7-9)  yields 


s If]" 


/iVi 

6n(Ui  +02)-^+  6n(U2  - Ui) 


3{6j  - 62) 


(7-12) 
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Figure  7-3.  Evaluation  of  Point  Velocity  for  Surface  i 


NO  FLOW 


Figure  7-4.  Model  to  determine  ifi 
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where 

h “ h2  “ h ^ 

Equation  (7-12)  is  the  appropriate  Reynolds  equation  for  rough  surfaces 
and  has  been  used  by  Elrod  [58]  and  Cheng  [53] . 

AVERAGE  REYNOLDS  EQUATIONS 

Equation  (7-12)  describes  hydrodynamic  lubrication  effects  for 
rough  surfaces.  In  order  to  utilize  the  Reynolds  equation  in  a practi- 
cal way  for  rough  surfaces,  it  is  necessary  that  the  solution  be  based 
upon  average  roughness  characteristics  rather  than  individual  rough- 
ness configurations.  This  method  predicts  the  expected  values  of 
pressure  and  flow  based  on  a minimum  number  of  surface  characteristics. 

There  are  two  methods  of  obtaining  the  average  solution  to  equation 
(7-12)  currently  in  use.  One  method  is  to  average  the  Reynolds  equation 
(7-12)  itself  prior  to  solution  of  the  pressure  distribution. 

Christensen  and  Tonder  have  developed  theories  for  averaged  equations 
such  as  those  used  herein  [43],  [44],  [47],  [48],  [49]. 

It  now  appears  that  although  the  Christensen  theory  for  longitudi- 
nal and  transverse  roughness  is  considered  to  be  accurate  [48] , the 
Christensen  [43]  and  Tonder  [49]  theory  for  isotropic  roughness  is  not 
accepted  (Elrod,  [58]).  In  a preferred-direction  roughness  arrangement, 
the  lealcage  flow  perpendicular  to  that  direction  does  not  vary  randomly. 
This  observation  allowed  Christensen  to  average  the  Reynolds  equation 
term  by  term  in  a mathematically  valid  fashion.  However,  in  the  case 
of  isotropic  roughness,  Christensen  and  Tonder  did  not  treat  both 
lealcage  flows  q^  and  q^  as  random  variables;  hence,  their  averaged 
Reynolds  equation  is  not  correct. 

Quite  recently  Patir  and  Cheng  [53]  have  developed  the  Average  Flow 
Model  based  on  the  second  concept  of  obtaining  average  results  for 
pressure  and  flow.  This  concept  involves  solution  of  the  Reynolds 
equation  (7-12)  for  a number  of  roughness  configurations  all  having 
identical  statistical  properties  and  then  averaging  the  results.  The 
model  is  practical  and  the  results  are  consistent  with  the  limiting 
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cases  of  longitudinal  and  transverse  roughness  evaluated  by  Christensen. 
This  method  is  considerably  more  powerful  than  the  Christensen  concept 
because  the  model  is  extended  to  include  anisotropic  roughness  con- 
figxirations  quite  easily.  However,  the  Average  Flow  Model  requires 
artificial  surface  generation  (or  knowledge  of  real  surfaces)  which 
is  not  required  in  the  Christensen  concept  of  equation  averaging. 

In  order  to  compare  the  two  concepts  of  obtaining  averaged  results 
for  equation  (7-12)  just  discussed,  the  model  developed  by  Cheng  and 
Patir  is  briefly  reviewed.  A control  volume  with  base  AxAy  and  film 
thickness  H = H2  - is  selected  to  analyze  the  expected  flow 
(figure  7-4).  It  is  assumed  that  AxAy  is  large  enough  to  contain 
numerous  asperities  but  small  relative  to  the  bearing  dimensions.  The 
unit  flows  are  given  by  equations  (7-6),  (7-7).  The  expected  unit  flows 
are  given  by  the  integral  averages  (q^  only) : 

y+Ay  y+Ay 

y y 

- ikt  “2  - "I’l'*''  ” 

Now  Cheng  and  Patir  define  the  pressure  flow  factor  and  the  shear 

flow  factor  A such  that: 
s 


2 / s 12n  X 3x  V 2 / 


(7-14) 


where  the  bars  above  the  symbols  indicate  average  values,  a is  the 

standard  deviation  of  the  combined  roughness.  The  $ term  is  eliminated 

s 

by  considering  the  case  of  pure  rolling  where  = ^2  = U.  Using  the 
boundary  conditions  shown  in  figure  7-4  with  the  nominal  film  thickness 
h constant,  it  is  easily  shown  that: 
y+Ay 

^ (7-15) 

h3  hR 

3x 
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where 


3x  Ax  (7-16) 

Hence,  i|>^  is  obtained  by  generating  random  surfaces  having  the  desired 
statistics  (standard  deviation  and  autocorrelation) , solving  equation 
(7-12)  for  p(x,y)  and  evaluating  (7-15)  given  (7-16). 

Using  this  method,  is  found  for  longitudinal  roughness  and 
transverse  roughness.  The  Christensen-Tonder  results  are  verified  in 
the  limit  by  the  Patir-Cheng  results.  However,  for  isotropic  roughness, 
the  Christensen-Tonder  equation  is  in  error  according  to  the  Average 
Flow  Model.  In  fact,  the  Patir-Cheng  results  show  that  a greater  flow 
resistance  will  be  encountered  for  isotropic  roughness  than  is  pre- 
dicted by  the  Christensen-Tonder  isotropic  model.  Thus,  the  Christensen- 
Tonder  model  underestimates  the  hydrodynamic  load  support,  and  this  is 
a major  inadequacy  of  the  model  used  herein. 

It  is  evident  from  this  portion  of  the  discussion  that  several 
fundamental  concepts  concerning  the  lubrication  of  rough  surfaces  are  in 
dispute.  However,  the  recent  work  by  Patir  and  Cheng  [53]  and  the  review 
by  Elrod  [58]  suggest  that  a better  understanding  and  more  precise 
models  will  be  forthcoming. 

REAL  SURFACES 

It  is  well  known  that  manufactured  surfaces  have  a microroughness 
structure  that  in  many  cases  is  closely  Gaussieui  in  the  distribution 
of  peaks  and  valleys.  The  rough  surface  structures  used  in  lubrication 
models  both  here  and  in  the  literature  are  also  assumed  to  be  Gaussian 
(or  a polynomial  approximation  thereof) . There  are  two  major  problems 
with  this  assumption  that  may  affect  predicted  seal  performance. 

First,  the  wear  effects  on  the  surface  structure  of  a seal  face 
during  run-in  certainly  alters  the  perfect  symmetry  of  the  assumed 
Gaussian  distribution  and  could  change  the  distribution.  It  is  expected 
that  surface  peaks  are  flattened  during  the  collisions  between  opposing 
asperities,  resulting  in  a severely  skewed  distribution.  If  the  peaks 
are  indeed  severed  at  a statistically  significant  height  value,  the 
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initially  Gaussian  pattern  will  not  be  skewed  but  actually  truncated.  j 

Such  a truncated  distribution  might  severely  affect  both  load  support  | 

and  leakage.  This  limitation  of  present  roughness  models  can  be  over-  : | 

come  using  the  approach  of  Patir  and  Cheng  applied  to  actual  worn  j 

j 

surfaces.  | 

A second  mechcmism  of  altered  surface  structure  is  termed  "grooving."  j 

Regardless  of  how  roughness  appears  on  one  seal  face,  it  is  possible  1 

that  due  to  ideal  wear  behavior  the  roughness  on  the  opposing  face  is  ; I 

a mirror  image  of  the  mating  face  (at  least  the  longitudinal  grooving) . : I 

This  implies  that  there  is  a correlation  between  the  two  surface  rough-  i J 

1 4 

nesses,  or  that  6 and  6»  are  not  independent.  For  an  ideal  correlation  | ■ 

^ M 

of  1 (perfect  mirror  image) , there  is  no  net  transverse  roughness  at  | | 

all.  For  any  non-zero  correlation,  the  concept  of  two  independent  sur-  | 

face  roughnesses  6^^  and  62  is  incorrect.  An  understanding  of  the  role 

of  grooving  will  require  an  analysis  of  many  measurements  obtained  from 

both  faces  of  seals.  < 

] 

EFFECTIVE  FILM  THICKNESS 

In  the  model  herein  it  has  been  assumed  that  when  h < 1,  the  rough-  ■. 

ness  distribution  is  simply  truncated  and  that  the  maximum  valley  depth  ! 

is  reduced  accordingly.  That  is  to  say  that  the  effective  film  thick-  j 

ness  E(H^)  is  continually  reduced  with  decreasing  h.  Furthermore,  it  | 

is  assumed  that  this  truncated  distribution  remains  constant  for  a j 

« I 

given  h < 1 with  both  time  and  wear.  1 

A fundamental  question  is  now  raised.  Assume  for  the  moment  that 
the  hard  moving  face  is  perfectly  smooth  and  that  all  wear  and  rough-  : 

ness  occurs  on  the  soft  face.  Now  assume  that  the  hydrodynamic  con- 
ditions are  such  that  asperities  must  carry  a load  such  that  30  percent 
of  the  rough  surface  must  be  in  contact.  Figure  7-5  shows  two  possible  ; 

distributions  that  could  provide  this  amount  of  bearing  area.  Case  a)  ; 

shows  the  present  assumption.  The  distribution  curve  retains  the  orig-  j 

inal  assumed  form  with  a maximum  range  of  2c.  The  peak-to-valley  height 
is  considerably  reduced  below  2c  and  therefore  E(H^)  becomes  smaller. 

Case  b)  shows  a possible  situation  which  would  give  quite  different 
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results.  Here  the  distribution  is  changed,  giving  a peak-to-valley 
height  that  is  much  larger  than  Case  a) . 

The  basic  difference  between  these  two  cases  is  that  in  Case  b) 
it  is  assumed  that  contact  changes  the  distribution  in  such  a way  as  to 
retain  a constant  or  near  constant  peak-to-valley  height.  In  fact,  the 
peak-to-valley  height  could  conceivably  remain  the  same  regardless  of 
bearing  area.  In  Case  a)  it  is  assumed  that  the  peak-to-valley  height 
is  reduced  with  contact. 

Such  different  behavior  could  arise  depending  on  the  precise  nature 
of  the  wear  involved  in  generating  the  surfaces.  Assume  that  the  wear 
is  due  to  abrasion  by  hard  foreign  particles  of  a given  size.  In  this 
case,  one  would  expect  Case  b)  to  be  more  correct.  That  is,  the  larger 
particles  plow  the  same  depth  of  groove  relative  to  the  hard  face 
regardless  of  the  bearing  area.  On  the  other  hand,  if  softer  abrading 
particles  of  the  soft  face  itself  are  responsible,  then  it  is  more 
likely  that  Case  a)  would  be  correct.  The  larger  particles  would  be 
crushed  very  quickly  and  not  have  an  opportunity  to  plow  deep  grooves. 

The  significance  of  this  discussion  is  that  the  pressure  buildup 

3 

in  a wavy  seal  depends  very  much  upon  E(h  ) when  particle  contact 

3 

occurs.  Given  the  possible  difference  of  E(H  ) for  the  same  amount  of 
contact,  the  type  of  wear  distribution  in  effect  could  significantly 
alter  the  hydrodynamic  load  support. 

Therefore,  in  hydrodynamic  lubrication  of  rough  contacting  sur- 
faces, it  appears  that  the  nature  of  the  wear  process  occurring  and  the 
resulting  roughness  distribution  is  of  considerable  importance.  This 
question  has  not  been  previously  addressed  in  the  above  context;  how- 
ever such  studies  are  being  made  as  a part  of  this  research  program. 

MICROCAVITATION 

Several  studies  have  shown  that  significant  hydrodynamic  load 
support  results  from  the  forced  lubricant  flow  over  the  tops  of  asperi- 
ties on  seal  faces  [13],  [14],  [15],  [16].  Cavitation  on  the  down- 
stream side  of  the  asperity  is  essential  to  developing  load  support. 
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Such  cavitation  will  be  termed  microcavitation  as  opposed  to  the  macro- 
cavitation associated  with  waviness  studied  herein. 

In  regions  of  a wavy  seal  where  average  or  expected  pressures  are 
relatively  small,  such  microcavitation  would  be  expected  to  occur. 

Present  rough  surface  models  which  predict  average  pressure  changes  due 
to  large  scale  nominal  film  thickness  changes  (such  as  the  model  herein) 
do  not  account  for  the  additional  load  support  caused  by  microcavitation. 
However,  load  support  is  developed  by  both  forms  of  cavitation.  A 
correct  model  must  account  for  both  types. 

This  shortcoming  can  be  corrected  with  a knowledge  of  the  distri- 
bution of  pressures  about  the  already-predicted  mean  value.  Given  this 
information,  regions  of  predicted  negative  pressure  due  to  microcavita- 
tion can  be  eliminated  from  the  distribution  and  a true  mean  pressure 
determined.  It  is  expected  that  using  the  method  of  Patir  and  Cheng 
[53] , a first  approximation  for  these  distributions  can  be  obtained. 

CONCLUSIONS 

The  questions  raised  in  the  previous  sections  indicate  that  the  de- 
velopment of  mathematical  models  which  describe  hydrodynamic  lubrication 
and  wear  for  rough  contacting  surfaces  is  still  underway.  Much  more  re- 
mains to  be  done  before  such  models  are  fully  developed.  To  date,  very 
little  experimental  verification  of  such  models  has  been  carried  out. 

The  use  of  the  experimental  seal  test  rig  described  in  the  next  chapter 
will  provide  some  data,  but  only  on  a level  of  measured  performamce,  not 
on  a level  of  understanding  the  microscopic  effects  themselves,  which 
are  pertinent  to  such  models.  Thus,  these  investigations  remain  to  be  done. 

In  spite  of  the  theoretical  limitations  of  the  model  as  pointed  out 
herein,  the  model  may  in  fact  be  an  excellent  tool  for  the  prediction  of 
overall  effects  such  as  friction,  leakage,  wear,  and  the  effect  of  sur- 
face finish.  This  question  can  be  answered  by  use  of  the  test  apparatus 
described  in  the  next  chapter.  If  the  model  proves  to  be  a valid  pre- 
dictive tool  for  mechanical  face  seals,  then  one  of  the  primary  objec- 
tives of  this  research,  that  of  developing  a useful  model  for  seal  de- 
sign purposes,  will  have  been  achieved.  However,  work  must  continue  on 
developing  a better  understanding  of  the  fundamental  questions  raised 
herein  in  order  to  completely  understand  face  seal  lubrication. 
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CHAPTER  8 


SEAL  TEST  RIG 

As  pointed  out  in  chapter  4,  only  very  limited  studies  have  been 
made  of  wavy  seals  operating  in  water.  To  provide  additional  experi- 
mental verification  of  the  model  described  herein  and  to  provide  fa- 
cilities for  further  studies  related  to  other  seal  problems,  the  design 
of  a mechanical  seal  test  rig  has  been  complete'd  as  a part  of  this 
project. 


BASIC  REQUIREMENTS 

A test  rig  designed  for  research  must  meet  more  requirements  than 
test  rigs  used  for  study  of  commercial  seals.  First  it  is  necessary 
that  a wide  reuige  of  speed  and  sealed  pressure  be  available  in  order  to 
verify  model  behavior  over  a wide  range  of  conditions.  Certain  other- 
wise uncontrolled  variables  such  as  shaft  runout  and  waviness  must  be 
very  carefully  controlled.  For  testing  the  model  of  interest  it  is 
necessary  that  a controllable  waviness  be  applied  to  the  seal  faces,  and 
in  order  to  reduce  the  effects  of  wear  this  waviness  must  be  moved  cir- 
cumferentially with  respect  to  the  primary  ring.  Leakage  from  the  seal 
must  be  measured.  Torque  must  be  measured  in  a manner  such  that  the 
seal  torque  is  measured  independently  of  bearing  torque  or  torque 
arising  from  a second  seal  used  in  a tandem  arrangement.  Many  of  the 
distortions  commonly  found  in  commercial  seals  must  be  minimized. 

To  meet  all  of  these  requirements,  the  design  of  such  a test  rig 
must  depart  somewhat  from  that  found  in  commercial  practice.  On  the 
other  hand,  it  is  important  that  the  rig  be  capable  of  simulating  com- 
monly found  seal  installations.  In  fact,  it  is  useful  to  be  able  to  in- 
stall a commercially  available  standard  seal  in  the  test  rig. 
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SPECIFICATIONS 


The  test  rig  specifications  which  satisfy  the  above  basic  require- 
ment are  set  forth  below. 


Configuration : 


Horizontal  shaft,  single  seal  with  torque  measured 
on  the  seal  housing.  Rotating  secondary  ring. 


Size: 


Accepts  up  to  3-5/8  inch  nominal  shaft  size  seal. 
4 inch  nominal  seal  face  diameter. 


Speed: 
Pressure: 
Bearings : 


500-4000  RPM.  9-70  ft/sec  for  a 4 inch  seal. 

0-1000  psi  continuous.  Vessel  is  rated  at  2000  psi. 


Five  ABEC-7  bearings  on  precision  spindle.  Bear- 
ings are  rated  at  a maximum  of  12000  lb  thrust  at 
4000  RPM  (corresponds  to  a 4 inch  seal  at  1000  psi) 


Motor: 
Fluid: 
Balance : 
Seal  Rings: 


5 Hp  variable  speed  drive  (belt  type) . 

Mate<^?  als  are  selected  for  seawater  service. 
Balance  ratio  of  1.0  (lower  values  can  be  used) 


Zero  moment  design.  Designed  to  minimize  the 
effects  of  pressure  caused  radial  taper. 


Waviness : 


Temperature : 


Seal  waviness  is  controllable  by  variable  ring 
distortion.  It  can  be  moved  relative  to  the 
primary  ring. 

Can  be  controlled  at  any  given  temperature  above 
100°F. 


Quantities 

Measured: 


Speed,  pressure,  temperature,  waviness,  torque, 
and  lealcage. 


DESIGN  FEATURES 


The  general  layout  of  the  test  rig  is  shown  in  figure  8-1.  The 
primary  and  secondary  seal  rings  are  shown  as  parts  ^7'  and  ^^31  . 

Shaft  1 turns  on  five  precision  angular  contact  bearings  designed  to 
operate  under  the  large  thrust  load.  Cooling  water  is  circulated  through 
the  test  chamber  41  through  openings  not  shown.  The  entire  assembly 
is  supported  cantilever  style  by  supports  attached  to  the  bearing 
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housing.  The  torque  transducer  acts  as  the  support  for  the  vessel. 

The  pressure  vessel  is  designed  to  be  readily  taken  apart  by  removal  of 

, A 

bolts  (43' . The  secondary  ring  support  rotor  is  removed  by  removing 
capscrew  . To  minimize  axial  runout  at  the  seal  face,  the  end  of 
shaft  is  ground  true  with  the  bearings.  Rotor  (^)  is  ground  both 
where  it  mates  with  the  shaft  and  secondary  seal  ring.  The  two  secondary 
seal  ring  faces  are  ground  parallel. 


i 


f 

i: 


The  primary  ring  is  shown  in  figure  8-2.  The  material  selected  for 
initial  tests  is  Pure  Carbon  P658RC.  The  unusual  proportions  of  the 
ring  result  from  adjusting  the  geometry  so  that  rotation  about  a circum- 
ferential axis  is  not  caused  by  pressure  variation.  The  primary  ring  is 
driven  by  lugs  which  engage  two  notches.  Two  as  opposed  to  a greater 
number  of  notches  were  used  because  this  arrangement  gives  the  primary 
ring  maximum  freedom  to  float  radially.  Second  harmonic  waviness  is 
minimized  by  driving  the  ring  through  the  centroid  of  its  cross  section. 
Any  remaining  second  harmonic  waviness  produced  by  drive  forces  will  be 
more  easily  flattened  compared  to  the  third  harmonic  waviness  built  into 
the  ring,  so  is  not  considered  important.  The  secondary  seal  is  lo- 
cated at  the  left  end  on  the  inside  diameter.  The  balance  ratio  for  the 
design  shown  is  unity.  Springs  to  the  left  of  the  seal  provide  preload. 


The  secondary  seal  ring  is  shown  in  figure  8-3.  Rings  are  fabri- 
cated from  tungsten  carbide  and  ceramic.  This  ring  is  also  of  a zero 
pressure  moment  design.  An  unusual  feature  of  the  particular  design  is 
that  the  mechanical  force  at  the  right  hand  side  of  the  seal  has  been 
reduced  to  nearly  zero  by  placing  the  Oring  seal  on  the  right  hand  face 
as  opposed  to  the  outside  diameter  (see  figure  8-1) . It  is  useful  to 
minimize  the  axial  force  to  reduce  coning  effects  caused  by  friction  at 
the  load  bearing  point  in  conjunction  with  radial  deformation  due  to 
pressure  and  temperature.  The  secondary  seal  ring  is  held  in  place  by 
spring  clips  (not  shown)  which  engage  the  notches  shown  in  figure  8-3. 


The  waviness  drive  mechanism  is  shown  in  figure  8-1.  A worm 
drives  gear  which  is  attached  to  cylinder  . Fixed  to  the  op- 
posite end  of  the  cylinder  is  the  waviness  ring  shown  in  figure  8-4. 

An  0-ring  is  placed  into  each  of  the  six  holes  on  the  ring.  These  O- 
rings  seal  between  the  primary  ring  and  the  cylinder  . A controlled 
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Figure  8-2.  Primary  Ring 
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Figure  8-4.  Waviness  Ring 


I 

L 


li 

J 


151 


pressure  is  supplied  to  each  of  the  six  cavities  by  channels  in  cylinder 
which  in  turn  are  supplied  by  a rotating  seal  in  the  left  hand  end 
plate  . The  controlled  pressure  in  the  six  cavities  creates  bending 
moments  in  the  seal  ring  which  cause  three  waves  on  the  face.  Since  the 
primary  ring  is  not  allowed  to  rotate  while  cylinder  rotates,  the 

wave  pattern  is  slowly  rotated  relative  to  the  primary  ring.  The 
0-rings  are  located  symmetrical  cd^out  the  centroid  so  as  to  produce  a 
zero  average  rotation  about  the  centroidal  axis. 


The  torque  transducer  is  shown  in  figure  8-5.  The  difficulty  in 
designing  such  a torque  transducer  is  that  torque  to  be  measured  is 
small  whereas  the  etxial  load  due  to  the  thrust  of  the  seal  is  large. 
Several  design  concepts  were  developed.  Figure  8-5  shows  what  appears 
to  be  the  best  compromise.  The  0.120  x 1 inch  members  carry  the  axial 
load  in  tension.  Due  to  their  small  moment  of  inertia  about  a radial 
axis,  these  members  have  in  combination  a low  torsional  stiffness.  The 
torsion  load  is  sensed  by  the  two  reduced  section  be2uns.  Strain  gages 
are  attached  to  these  members.  The  axial  load  effect  is  isolated  from 
these  members  by  the  cutout  shown.  The  torque  transducer  does  have 
some  small  aunount  of  sensitivity  to  axial  loads  because  the  0.120  x 1 
members  do  attempt  to  straighten  (when  twisted)  due  to  the  axial  load. 
This  effect  is  largely  offset  by  the  fact  that  the  torque  sensing  beams 
have  a very  high  torsional  stiffness.  Thus  the  actual  torsional  de- 
formation of  the  unit  is  quite  small. 


The  spindle  is  designed  using  ABEC-7  angular  contact  bearings.  Four 
bearings  in  tandem  are  required  to  support  the  thrust  load.  Recirculated 
oil  is  supplied  to  the  housing  from  an  oil  cooler.  The  right  hand  fifth 
bearing  serves  to  preload  the  bearings. 


The  vessel  cylinder  euid  end  plates  are  made  of  monel  400  for  good 
resistance  to  seawater.  All  other  parts  exposed  to  the  working  fluid 
are  made  of  316  stainless. 


Pressure  to  the  vessel  is  supplied  by  pressurizing  an  accumulator 
with  nitrogen  gas  at  the  specified  pressure.  Makeup  water  to  offset 
leakage  is  supplied  by  a slow  acting  air  driven  piston  pump.  Flow  cir- 
culates through  a heat  exchanger  which  is  cooled  by  tap  water.  Tem- 
perature is  controlled  by  regulating  the  cooling  water  flow.  Waviness 
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Figure  8-5.  Torque  Transducer 


L 


control  pressure  is  supplied  by  nitrogen  gas 


EXPECTED  PERFORMANCE  ! 

Using  the  isotropic  waviness  model  of  chapter  2,  expected  per- 
form2uice  curves  have  been  obtained  for  the  test  seal  for  one  set  of 
operating  conditions.  These  conditions  are: 

= 1070  psi  ' 

= Pc  = 0 

= 200/s 
= 20  yin 
= 2.088  in 
= 1.900 
= 38000  psi 

= 9.9  • 10  ^ lb  s/in^  (100®  F water) 

= 0.0441  in'* 

= 3 X 10^  psi 
= r 

O 

= 8.14  I 

(' 

1 

= 3 ! 

i 

= 1.0 

The  corresponding  dimensionless  parameter  values  are 

p = 0.005  I 

o I 

p = 0.176  I 

m I 

r.  = 0.910  i 

^ I 

r /c  = 1.04  10^ 

O I 

K = 3.1  X 10  ^ 

The  performance  curves  are  shown  in  figure  8-6.  At  this  high 
operating  pressure  the  percent  fluid  film  load  support  increases  from 
52  percent  to  only  about  80  percent  as  h^  increases.  At  lower  operating 
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Figure  8-6.  Predicted  Performance 
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Curves  for  Test  Seal 


pressures  the  fluid  film  load  support  would  be  higher.  Figure  8-7 
shows  the  required  initial  waviness  for  these  same  conditions.  For 

- »7 

the  test  seal,  K = 3.1  x 10  . This  curve  shows  that  a condition  of 

instability  will  be  reached  approximately  when  h^  » 1.0.  A wide  range 
of  leakage  and  friction  behavior  can  be  obtained  over  the  stable 
operating  range. 

STATUS  OF  THE  RIG 

As  of  this  report  date  the  entire  test  rig  assembly  shown  in  fig- 
ure 8-1  has  been  designed  in  detail  and  parts  are  being  made.  The 
auxiliary  systems,  support  stand,  and  controls  are  nearing  design  com- 
pletion. It  is  expected  that  the  test  rig  will  be  in  operation  on  a 
limited  basis  in  May  of  1978.  Full  operational  capabilities  are  ex- 
pected in  June  or  July,  1978. 


CHAPTER  9 

SUMMARY  AND  CONCLUSIONS 


Many  detailed  conclusions  have  been  made  throughout  this  report. 

The  purpose  of  this  chapter  is  to  pull  together  some  of  the  various 
ideas  and  make  more  generalized  conclusions  and  to  provide  a short 
summary  of  the  most  important  features  of  this  work. 

1)  A hydrodynamic  lubrication  seal  model  has  been  developed.  The 
effects  of  waviness,  roughness,  asperity  contact,  wear  and  elastic 
deflection  have  been  included.  A numerical  technique  has  been 
developed  to  accurately  and  rapidly  solve  for  the  seal  equilibrium 
solution.  Given  initial  waviness,  period,  viscosity,  speed,  size, 
pressure,  material,  euid  roughness,  the  model  can  be  used  to  find 
the  fraction  of  fluid  film  load  support,  relative  wear  rate, 
pressure  distribution,  friction,  and  leakage. 

2)  Results  based  on  assuming  longitudinal  roughness  show  that  98% 
plus  of  the  load  is  carried  by  fluid  pressure  for  heavily  loaded, 
low  viscosity  seals.  This  result  is  unrealistic  and  shows  that 
pure  longitudinal  roughness  is  not  a correct  assumption  for  seals 
in  spite  of  their  appearance. 

Results  obtained  using  the  isotropic  roughness  assumption  appear 
to  be  realistic.  Using  the  isotropic  model  developed  herein,  it 
is  expected  that  predicted  fluid  pressure  load  support  may  be  slightly 
lower  than  in  reality.  Certain  refinements  to  the  roughness  lubri- 


cation model  are  needed  to  provide  more  accurate  results.  Theo- 
retical and  numerical  developments  which  will  correct  these  de- 
ficiencies are  being  made  by  these  authors  and  others. 

3)  Results  obtained  using  the  isotropic  model  have  been  compared 
to  experimental  results  available  in  the  literature.  The  comparison 
is  generally  good.  However,  more  complete  experimental  data  are 
needed  for  a thorough  comparison. 
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4)  The  results  from  the  model  show  that  surface  roughness  is  the 
most  important  parameter  in  determining  the  amount  of  fluid  film 
load  support  caused  by  waviness  in  conjunction  with  hydrodynamic 
effects.  If  the  roughness  is  large,  mechanical  contact  will  sup- 
port most  of  the  load.  If  roughness  is  low,  much  of  the  load  will 
be  supported  hydrodynamically  even  with  low  waviness.  The  range 
of  roughness  values  examined  herein  is  representative  of  the  range 
of  values  found  in  seals.  This  suggests  that  a considerable  varia- 
tion in  hydrodynamic  performance  would  be  expected  in  seals  due  to 
differing  surface  roughness  effects  alone.  Surface  roughness  and 
the  material  parameters  that  affect  it  will  require  much  more  ex- 
perimental study  before  the  effects  of  materials  on  seal  hydrodynamics 
are  completely  understood. 

5)  The  second  most  important  parameter  affecting  seal  hydrodynamics 
is  the  net  waviness  at  the  seal  faces.  Generally  as  net  waviness 
increases  in  amplitude,  friction  decreases  and  the  percent  fluid 
fluid  film  load  support  and  leakage  increase.  As  waviness  is  in- 
creased, there  is  a distinct  transition  between  pure  hydrostatic 
operation  and  partial  hydrodynamic  operation.  If  surface  rough- 
ness is  sufficiently  small,  the  load  will  be  supported  completely 
by  fluid  pressure  even  at  a small  net  waviness.  However,  even  at 
large  fractions  of  fluid  film  load  support,  the  faces  still  touch 

at  the  tops  of  the  waves. 

6)  Hydrodyneimic  effects  can  easily  be  overshadowed  by  hydrostatic 
effects  caused  by  radial  taper.  For  example,  with  a slight  con- 
verging radial  taper,  the  hydrostatic  load  support  increases 
markedly  over  its  parallel  face  value  and  may  easily  cause  as  much 
additional  load  support  as  the  hydrodynamic  effects.  On  the  other 
hand,  a divergent  taper  will  reduce  hydrodynamic  effects.  The  load 
is  carried  mechanically  at  the  outer  radius  and  the  average  film 
thickness  increases.  This  reduces  the  hydrodynamic  load  support. 
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7)  Whether  or  not  significant  hydrodynamic  effects  occur  in  com- 
mercial seals  depends  on  several  factors.  If  as  running  roughness 
is  low,  the  seal  rings  are  stiff,  and  a moderate  amount  of  waviness 
is  present,  it  is  likely  that  hydrodynamic  effects  are  occurring. 
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even  in  low  viscosity  heavily  loaded  applications.  On  the  other 
hand,  if  the  rings  are  not  stiff  (e.g.,  made  of  solid  carbon)  it 
is  likely  that  hydrodynamic  effects  are  not  occurring,  even  for 
large  values  of  initial  waviness.  The  waviness  simply  flattens 
out  under  load. 


For  higher  viscosity  applications,  then  hydrodynamic  effects 
might  occur  even  in  low  stiffness  seals.  The  question  can  be 
answered  for  a given  seal  system  by  using  the  model  developed 
herein.  For  the  particular  carbon  seals  examined,  it  was  con- 
cluded that  significant  hydrodynamic  effects  would  not  be  expected 
at  least  at  higher  pressure  values. 
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It  was  found  that  low  stiffness  high  waviness  seals  could 
operate  under  two  equilibrium  conditions,  one  at  low  leakage  and 
one  at  very  high  leakage.  High  stiffness  seals  do  not  have  such 
instabilities  but  the  operation  point  does  become  very  sensitive 
to  the  value  of  initial  waviness. 


8)  Hydrodynamic  load  support  slowly  decreases  with  time,  after  an 
initial  increase,  due  to  the  effects  of  wear.  Thus,  the  long  term 
effect  of  waviness  in  a steady  state  operation  is  zero.  The  time 
period  over  which  this  occurs  depends  on  many  parameters.  This 
behavior  can  be  overcome  by  design. 

9)  As  to  the  potential  of  utilizing  hydrodynamic  effects  to  ad- 
vantage by  design,  the  results  show  that  wear  rate  and  friction 
can  be  reduced  significantly  while  maintaining  leakage  at  accep- 
table levels.  This  can  be  accomplished  by  moving  the  waviness 
relative  to  the  ring  so  that  the  effect  of  wear  is  minimized. 

10)  The  effects  of  roughness  on  hydrodynamic  lubrication  are  not 
completely  understood.  Better  models  of  this  interactive  behavior 
are  now  under  development.  Some  questions  raised  herein  are  still 
to  be  addressed. 

11)  A test  rig  has  been  designed.  Using  this  apparatus,  it  will  be 
possible  to  obtain  data  needed  to  verify  the  seal  model  herein  as 
well  as  to  explore  many  other  seal  related  questions. 
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APPENDIX  A 


REYNOLDS  EQUATIONS  AND  ASPERITY 


BOUNDARY  RELATIONSHIPS 
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1)  General  form  of  Reynolds  equation: 
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Coefficients  are  functions  of  6,  r.  In  finite  difference  form: 

P(I,J)  = AIM1(I,J)*P(I-1,J)  + AIP1(I,J)*P(I+1,J) 
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2)  Reynolds  equation  for  longitudinal  roughness.  Film  thickness  H 
a function  of  6 only.  (See  Chapter  2 for  actual  equation) . 
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For  finite  difference  purposes; 
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3.  Reynolds  equation  for  isotropic  roughness,  H 
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= H(9)  . 


The  finite  difference  coefficients  are  otherwise  the  same  as  above 
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4)  Asperity  boundary  transition  equations  from  full  film  to  contactinq 
region.  Longitudinal  roughness. 
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ASPERITY  REGION 
H < 1 


Coefficients  A-E  already  defined  for  longitudinal  roughness. 
Finite  difference  coefficients  for  the  boundary  points: 


Point  a)  AIM1(I,J)  = 


(2  - 

2 


AIP1(I,J)  ^2  - X)  (1  - X)(t>[^g2  2A 


r_^  (1  - x)B'i 

- 2A6  J 


AJM1(I,J)  = i - -2- 


B(I,J)  = - T 


Point  s)  AIM1(I,J)  = T- 


1 -2 

1 Ar 

2Ar 

D 

-2 

Ar 

2Ar 

E 

2A 

2 ' 

- x)Ae 

A 

(1  - 0)Ae^ 

B 

AIP1(I,J) 


Point  a) 


Point  s) 


Point  c) 


AIP1(I,J) 

AJM1(I,J) 

AJP1(I,J) 


(2  - X)  (1 

2 


(2  - X)4, 

c 


+ ■ 


= 

AIM1(I,J)  = 

AIP1(I,J)  = 

AJM1(I,J)  = 

AJP1(I,J)  = 

B(I,J)  = 


2A 


(1  - X)A0' 


X(t) 


Ae 

2 


(1  - X)* 


Ar 

C 


-2 

Ar 


+ ■ 


AIM1(I,J) 

AIP1(I,J) 

AJM1(I,J) 


2A 

X(1  - X)Al 

2 

(1  + X)  4i 
2 

X(1  + X)i{) 

AJPl  (I, J) 
E 


B(I,J) 


ooorir»oonnrir>ririor»rioo 


APPENDIX  B 


ISOTROPIC  ROUGHNESS  SOLUTION  - COMPUTER  PROGRAM 
BASE  CASE  2 SOLUTION 

PN0GRAH  NAME  - APPRaxS 

PReGRAM  T2I  calculate  CQJILIBRIJ'1  aPERATINC  C0^3iri0MS  PaR  an 
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FACE RajGHNESS  AMD  RADIALLY  PARALLEL  FACES 
ALL  Parameters  are  dimemsionlcss 
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UEL»1 ./PSBAR/RaaC 

Ma«0,R15 

HME AN*m0 

waPIE*(l,-RI»RI) •(P0*BR^PI*(l,-aR) ) 
nRITE(6,27)m 

f FaRHAT(SX,*REaUlREO  L0AD  Ha*,E12,b) 
CALCULATE  HESH  P0INTS 

0T«2,aPIE/NHAVE/(IMAX-l) 

DRa(R0-RI) / ( JMAX-l) 

tst«pie/nhavb 

TaTST-OT 

00  100  Ial,lHAX 

TaTaOT 


BESl  AVA1U8LE  COef 


B-l 


100  T0(n«T 

HmdlmOU 

[)3  110  J>1,JMAX 

110  RS(J)aR 

D0  10  iBl.MTERiS 
HjA(I)a0. 
in  HlBtniO, 

HIA(1)«.2 

JMM«JMAX-1 

C SET  initial  GJESSCS 
CUESS*(P0+»I)A,5 
od  111  i«i,nAx 
P(I»1)*PI 
P(I»  J*<AX)*P0 
00  111 

P(I» JIiGJESS 

111  C0NTINJE 
NRITECbfDNNAVEfPO 

1 F0HmaT(5X. 'NWAVE  I2»/,5X, 'Pa  '»1X,E10.4) 

wRITE(b,2J^lIA(n#RI 

2 F0«MAT(5X, 'HIAdJ  MX,El0.a,/,5X, '^1  *,1*.E1D,4) 

NRITEfh,3)R00C,PH 

3 PaRMAT(5X,»R/C  lX,E10,a,/,5X, 'PM  »,lX,E10.a) 

CALL  <4L0AD(M,Hai 

CALL  FRICT(PSBAR,RaaC,RF, AMJJ 
i*RITE(b,a)  AMJ 

a F0RMAT{/,3X, 'C0EFFICIENT  0F  F^I C T I aNa ' , E l 2, 4) 

CAUL  PRSHA4 
900  ST0P 
ENi) 


BESI  AVAIWBLE  COPY 
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subRMjrnc  EdNS/wC^r^ESAPfEPs# ^Axir^irEiys, jteo) 

C GIVE^I  H0  THIS  SJBPfT2Si9AH  SOLVES  P0R  THE  PRESSJRE  0I S T ) I B J T I 
C INCLJOlNS  cavitation 

COHMON  0T,3ft,T8(192}»RG(13).  IHAX,  JHAXfHASPtnZ) 

C0MHON  T,R,RI,^0,AI,R0,PC,P{132,1 J) ,PIE,HAAVE»PFAC,PHLS 
COMMON  HHEAN,HIA(13) ,HXB(1Q) ,3EL»NTERHS,IAJP,IA0N 
C0MM0N  H,OH,CH,OCH,EHS,OEHS,EiaHS,9I,OBX,EHL0AO«PSAR(192} 

C0MM0N  AIHKXOS,  IS)  ,AIP1(102|13)  , AJHl(l02«t9)  ,AJP1  (102,  13) 

COMMON  B2HS(i02, IS) ,HL0(i02) ,HLeAO>aiN,oaur,TarL3iEieM 
nlMENSlBN  lUP(lS) ,IOHN(13),HL0O(lO2),aOXN(lO2},930jT(132} 
DIMENSION  ANSAi/(1S)  ,PMAT(SO,S)  ,BMAT(S3)  , AN&2(1S,  «) 
iSrOPalHAX-l 
J5T0P«JHAX-1 

C initialize  all  EQJATION  COEFFICIENTS 
DO  S12  Ja2fJST0P 
00  990  I«l,ISTaP 
TaTSd) 

CALL  NAVSHP 
AAsBI*EHS/9G(J) 

HRb(BI*0EH3^EMS*3SI}/9&(J) 

CCaBI*R6(J)*EHS 

00«BI*EM3 

EEa6,*9S(J)*(ai*3EH»EM*0BI) 

AlJa2,a(AA/DT/DTtCC/D9/Drt) 

A I Ml (I,J)a(AA/3T-93*,5)/DT/Al  I 
AIPI (I, J)«(AA/3Tt83*,5)/0T/AI  J 
AJHl (I, J)b(CC/3R-03*,5)/DR/AIJ 
AJPI (I, J)a(CC/DR*D3*.S) /UR/AI J 
BRHSCI, J)a-EE/AlJ 
990  CONTINUE 

bl2  continue 

IMMlalMAX-l 

JM2*JMAX-2 

IF( JTER.NE.O)  30  To  «l 

C this  sequence  FR0H  here  T0  1329  IS  JSEO  THE  FIRST  TlMc  through  This 
C SUBPROGRAM  TO  OSTaIN  AN  INITIAL  3JESS  FOR  THE  CAVITT  SHAPE 
U0  999  J«1,JMAX 
IUP(J)»0 
999  I0wN(J)«0 
jTOPaO 
OMEGAiI.S 
flM0«i  ,»aHE3A 

C this  do  L<)0P  is  a successive  0VERRELAXATI0N  LOOP  TO  OiTAlN  AN 
C initial  GUESS  FOR  The  PRESSURE  DISTRIiJTION 
00  1300  1I1b1,3 
DO  50  JJai.JM? 

JaJMAX-J J 

nil 

I2*IMM1 

DO  «0  Iall,I2 

IMl«I-l 

iJ("5i1e3,o)  iHi.iMMi  BE5T  AVAILABLE  COPY 


B-3 


r*  r»  o oooocin 


T 


IF(IPl,Ea,IMA*)  IPl»l 
TEMPiPCI, J) 

Pdf  J)*P(IMl,  J)*AIPl(I,  J)«PCIP1,  J) 

PCIf  J)»PCIf  (I#J)«»Pd,  J-U 

Pdf  J)»PCI»  J)*AJP1(I,  J)*P(I,  J*l)*8^HSd,  J) 

Pdf  J)>Pdf  J)*a^CGA«0»10*TEMP 
ttO  C0NT1NJE 
50  C0NTINUE 

1000  cbntxnue 

C IN  THE  SEQUENCE  HERE  T0  600,  THE  INITIAL  5JES8  FSR  THE 

C CAVITY  SBUNOARV  IS  DEFINED  BY  ALL'  PRESSURE  POINTS  HhERE  P<PC. 
JMt«JMAX»l 
00  600  JaZfJHl 
NUP«1 

00  700  HlflHMi 

G0  TB  (7Dlf702,6D0}fNUP 

701  IF(P(If J),LT,PC)  S0  T0  703 
60  TB  700 

703  IJP(J)«I-1 
NUP«2 
JT0PbJ»1 
G0  T0  700 

702  IFlPd,  J)  ,GT,PC)  G0  T0  704 
G0  10  700  • 

704  inwM(J)»I 
NJP«3 

700  C0NT1NJE 
600  C0NTI>f'JE  . 

00  1329  J«1,JHAX 
ANGAV(JJ»-.5*PIE 
00  1329  VN»1,4 
1329  ANG2(Jf NV)iO,0 
41  J0UT«1 

BHtGAiBHEGAP 

0M0al,«0HESA 

C 00  710  IS  The  PRIHARY  successive  aVERRELAXATlaV  L00» 

42  00  710  III«lfHAXIT 

JT0P  « THE  B0TTaHH03T  R0H  aF  NON-CAVITY  PRESSURE 

I1C2  - THE  T0PM8ST  CAVITY  R0H  CENTRAL  PaiNT 

TUP(J)  • THE  RIGHTHAND  MOST  V0N  CAVITY  PRESSURE  PalviT  13 

THE  LEFT  0»f  THE  CAVITY 

IDtYN(J)-  THE  LEFTHANO  M0ST  N0N  CAVITY  PRESSURE  P0INT  T0  THE 
RIGHT  0F  THE  CAVITY 
IF(JTaP,LE,l) JT0P«2 
IUP0«1UP( JT0P-1) 

IIC2«(IOHN(J.Ta»-n-IJPO)/2*IJ»(JT0B.l) 

IICl»IOHM(JTaP-l) 

MT0P«1 

THE  SWEEP  IS  R0W  SY  Raw  STARTING  FR0M  THE  BjTStOE  RADIALLY  IVward 
MT0P«1  • DURING  THE  FIRST  CAVITY  R0W  REACHED 
■2  • AFTER  the  FIRST  CAVITY  R0W  REACHED 
00  720  JjaifJMB 
I«0 


I 
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BEST  AVAIUBLE  COPY 


o o n orionoo 


NUPat 

J8JMAX»JJ 

NUPal  . at  the  start  0F  A RSIH  AND  UBBKINB  FBR  P<>C  • NaRHAL 


REwAXATiaN 

ai  • aXTHXN  THE  CAVITY  LeBRiNS  FflR  IOhN(J) 
a)  • re  THE  RI6HT  OF  THE  CAVITY  ••  N0RHAt  REtAXATIBN 
■4  • SEaRCHINB  for  CBNDITlaN  TB  DEFINE  lOHN 

b5  • OPING  special  ROM  JT0a 
DB  752  IlaUlSO 
I8l41 

1F(NUP»EB.5)  SB  Ta  1105 
IFCJ.EO.JTBP)  BB  TB  1104 
GB  TB  1105 

1104  1F(I,6T,IIC2)  NUPaS 

1105  continue 
IHlal.l 

IPialal 

IFClHl.Ca.D)  iHlalHMl 
ZFCIPl.ES.IMAX}  IPial 
GB  TB  (721, 72S|721,7?1, 1103), NJP 
NUPa5  • FRBH  HERE  TB  1150  IS  F0R  A SPECIAL  R04  JaJTBP,  S»EClAt. 

CAVITY  aajNOARY  EOUATIBNS  ARE  USED  IN  THE  RCLAXATiaN 
FRBH  IIC2  TB  ZZCl 
1103  lF(J«CQ,JTaP)  SB  TB  110b 
ZF(J,Ea,2)  GB  TB  1150 
IIC2alDNN(J) 

IIClaI0NN(J«l) 

110b  JFd.EQ.IlCl)  GB  TB  1150 
NINCaXZCl«IIC2 
ifcninc.le.i)  SB  ra  ii5o 

ANGAaANGAV(J) 

IF{J,E9«JTBP)  ANSAao.SaPlE 
c angle  OEFINEO  in  corner  point  B^ERATIBN  is  jseo 
c angle  is  interpolated 

ANGBaANGAV(J-l) 

0AN6aAN68«AN5A 

ZFRACaZ-IlC2 

ZFRACaZFRAC/NiNC 

AN6aAN6AAZFRAC*DANS 

IIiaIUP(J*n 

TarSdll) 

CALL  RAV3H3 
LHlaEH 
BHlaBI 
TaTGd) 

CALL  HAVSna 
BHZaBZ 

PPPNab,aRG(J)aCBS(ANG) *(EH*flHl/BH2aEHl) /EHS 
PPPRaPPPNaSlN(AN3) 

ANUHaP(I-l,  J)aAIHl  (I,  J)»PdAl,  JlaAlPl  d,J).DRaPPaR#AJHi  f I,  J) 
ANUHaANUMtPd,  JAI)  aAJPl  (I,  J)4SRHS(1,  J) 

Pd,  J)aANUH/(l,*AJHl  d,J)) 

Pd,J-l)aPd,  J)-0RaPPPR 

BEST  AVAILABLE  COPY 
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p(i,  j-n»pc 

Ge  TB  750 
llSO  NUPaS 

C 721  TB  723  IS  THE  HELAXATIBN  STEP  APPLIED  Ta  THE  &EHE?At 
C CAVITY  PBXNT 
721  TEHPaP(X,J) 

P(X*J)aAIPl(If  J)*P(IM1,J)«AXP1(XW)*P(1P1» 

Ptif J)aP(X,J)^AJHl(X»J)aPtI, J-l) 

P(Xf  J)aP(XiJ)A4JPl(I»J)*P(Xi  jAl)«'f}^HS(Xi  J) 

P(It J)aP(If J)*3^E6A«aM0*TeHP 
6B  TB  (723, 72b, 1100, 769), SUP 
C SUPal 

723  CINTXSUE 

ZF(P(I,J),lT,PC)  Ga  TB  725 
GO  TB  1100 

C SETS  UP  SEARCH  FBR  I0i4N(J)  IF  SBT  PREYIBUSLY  OCFXSEO 

725  lUP(J)aI«l 
IF(HTiP,B0,l)  JTaPaJai 
NUPb2 

XF(IDHS(J).C9.0)  SJPaa 
IF(I0HS(J).LE.Z)  SJPaa 
P(l,J)aPC 
GB  TB  730 
C NUPa2 

726  IFd.Ea.ZOHNCJ))  63  TB  727 
P(Z,J)aPC 
GB  TB  730 

C SUPaR 

769  IF(P(Z»J).3T,PC)  63  TB  768 
P(I,J)aPC 
GB  TB  730 
766  10KS(J}al 

C SETS  UP  LB»»ER  rO-^S  VALUES  IF  S3T  OEFISEO 
1F(IDHN(J-1),E3,0)  GB  TB  5326 
GB  TB  5327 

5326  jSaJ.l 
UB  5323  J3a2,JS 
jaaj,>j3«l 

1F(ZDHS(J4) .EOtO)  ZDHS(J4)aX0HS(JA»l)tl 
5323  CBSTZSJE 

5327  lal-i 
NUPa2 
60  ra  730 

c HEGIN  CbRHER  NBOE  AOJJSTHtNr  HERE  • ESO  AT  1133 

727  IFIMTBP.EU.n  63  73  789 

IF(I04H(J)  ,UE,I04S(J4>1)1  GB  TB  788 
GB  TB  789 

788  X0PH(J)«ID4N(J«l)tl 
PCI, J)aPC 
GB  TB  730 

789  SUPaS 

C SSTPbO  • HAS  SBT  HBVEO  CBRNER  SBOE  YET 
C al  - HBVXSG  CaRSER  SBOE  TB  LEFT 


B-f) 


BEST  AVAIUBLE  COPY 


non  V non  non 


■ e - CaO^Ei)  N0OE  T0  <tmHT 

■ 3 . indicates  a reversal  " CaRNER  N0OE  IS  <iavE?  SACK  and 

N0yE'<ENT  IS  ST0PPED 

NSTPaO 

740  lClaI0wN(J«n 

IF(J,Ea.2)  lCt>I0i^M(J) 

tC2«I0NN(J«n 

NSP«i 

IF(HT0P,nE,1)  50  Ta  742 
02?  IUP0«IlIF(J) 

IC2«(tDWN(J)-IJP3)/2tIUP(J) 

742  NlNCaXCl-IC? 

DELFa«2.*0R 
1F(J.E3.2)  OEi.)(*«OR 
DELYa»NlNC*DT*R6(J) 

ANGLE  IS  DEFINED 

ANCaATANECOELYrOELK) 

angles  are  set  phen  Nar  assigned  during  first  FER  iterations 

00  1302  NN«1,4 

lF(ANC2(Jf  NN},E0,0,0}  ANG?(J, NN}aANG 
1302  CONTINUE 

angle  is  defined  as  average  0F  angles  FOR  FaUR  iterations 
AN62(Jf 4)aANG2(Jf 3) 

ANC2(Jf 3}BANG2(Jf2) 

ANG2(J,2}aANG2(J,I) 

ANG2(J»n>ANS 

ANG>0.2S4(AN62(J»  n«ANG2(J,2)  ^ANG2(J,  3)  4-aNG2(  J,  4) ) 

24  ANGAV(J)b^ANG 
XUalUP(J) 

T«T5(III) 

CALL  ^AVSHP 
EHl.Ert 
BHjaai 
IClalDNN(J) 

TaTG(ICl) 

CALL  NAVSH9 
BH2aai 

needed  PARTIAlS  are  found  based  on  CAVITT  BaUNDARV  CavDITIi^NS 
PPPNaS,*R6(J)  *C0S(ANG)*(EH-Bm/BH2*EH1)  / EH3 
744  PPPTaRGCJ) APPPNaCOSCANG) 

PPpRaPPPNaSINtANG) 

A MATRIN  aF  CaEFFICIENTS  FDR  THE  CORNER  NaOE  AND  POINTS  JJST 
TO  ITS  RIGHT  ARE  SET  JP 

PNATCl  ,2jaAIHl (I, J)»AJMI (I, J)-l, 

PHAT(l,3)aAlPUI,  J) 

BHaT (l)a-BRHS(I, J)tOTAPPPT*AlNJ (I , J) tORa^PPRaA JNj ( i, ji 
BHAT(naSHAT(l)«P(I,  J^DaAJPI  (I,  J) 

PPPToaPPPT 

NEOaXOHN(J.n«ID4N(J) 

IF(J.E0,2)  NEOal 
IZTaNEQ*! 

IF(IZT.lE.3)  60  Ta  7300 
00  6001  iZaUJZT 


j 

i 
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BEST  AVAIUBli  COPY 


iRewsiz^i 


C IN  between  ANSLES  are  INTERP0LATCO 
• ANGAsANGAVCJ) 

AN6BbAN6AV(J*1) 

0AN6aAN68*AN6A 

ZOUMiXZ 

NiNCaNCQ 

ANGBZOUNtOANQ/NlNCaANCA 

TaTGtXN) 

CALI  4AVSHR 
BH2bBX 

PPPNa6,*R6(J)*CeS(AN6I«(eH«BHl/SH2*EHn/EHS 

PPPRaPPPN*SIN(AN6) 

PHAT(XRa><,2)aAJHt  (IN,  J)»|, 

PHAT(XReN,l)aAXNX(XN, J) 

PMAT(XRa«l,3)aAXPl(XN,J) 

BHAT(IR0N)a-BRHS(lN,  JJtDRaPPPRaAjNl  (IN,J)-PlIN,J*n*AJ»l(l\,J) 
IFdZ.EQ.XZT)  PMAT(XRBN,3)aO. 

IFdZ.EO.lZT)  3NATdRaKOaBMArdRB<4)-P(IN#l,  J)aAiai  (IN,  j) 

(»oni  caNTiNjE 

C CBRNER  AND  ADJACENT  PRESSURES  ARE  RELAXED  5 TI^eS 
DB  1493  lASaifS 
INaltNEQ-l 

PCIN,  J)a(«aMAT(NEa,  n*P(IN-l,  J)^BHAT(NEO)  J/PNATlNEJ.aj 

IFdZT.LT.E)  Ca  Ta  1492 

DB  6002  XZa2,XZT 

XNbI»NEO*IZ 

lRawaIN«Z4l 

PdN,  J)aPMATdRa4,l)»PdN-l,  J)»PMATdR0a,S)*PdNfl,  J) 

PtlN,  J)a(-PdN,  J)f3MATdR0Nl)/PMAT(IR04,2) 

6002  CaNTlNUE 

1492  CBNTINJE 

: THE  CvJRNER  PaiNT  PRESSURE  IS  CALCULATED 

Pd,  J)a(-PHAT(l,5)»P(Itl,  J)*S4aT(1)  J/PMAT(1»2) 

1493  CPNTINJE 
Ga  Ta  7001 

7000  C0NTINUE 


Pd,  J)a(8MAT(n-PHAT(l,3)*P(Itl,  JD/PHATd.Z) 

7001  PCI-l,  J)aPd,  J)-OT*PPPT0 
if(nstp,eo.3)  sa  ra  749 

: FIRST  caNoirxaNs  aF  corner  point  pressure  eojal  APPRaxiHATEL^ 

C CAVITY  PRESSURE  IS  CHECKED 

IF(P(I, J) ,5T,PC)  Sa  T0  731 
Pd-1,  JJaPC 

C CORNER  PBINT  IS  HOVED  T0  RIGHT 
lalti 


1115 


IDwN(J)al 

iF(j,E0.2)  G0  ra 

IF(I,GT,(I0WV(J.1 
IF(NSTP.EQ.l)  6a 
IF(J,NE,(JTeP-i)) 
IF(NSTP,E0,0)  6a 


Ills 

}-l)}  60  T0 

ra  733 

60  T0  150S 
T8  ISOB 


SOI 


Ta 


best  avaiuble  cm 
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ce  T0  doi 
1506  NST^aa 

Ge  T0  749 
601  1>I«1 

IOWN(J)«I 

MTld6«2 

6i  re  7330 

.733  NST6«3 

IOMN(J)aI 

Ge  re  74a 

C SECOND  CONOXTleN  eF  CONNER  POINT  PRESSURE  ESUAL  APPROtliATELT  TO 
C CAVITY  PRESSURE  IS  CHECKED 
731  1FCP(X*1,J),LT,PC)  BO  TO  749 
P(I-1,J)«PC 

C CORNER  POINT  IS  HOVEO  TO  LEFT 

m-i 

IOWN(J)aI 

lF(I,LT.(I9WN(jai)fl))  GO  TO  833 
IF(NSTP,EQ,a)  30  TJ  741 
lF(J,NE.(JrOP-l))  SO  TO  1509 
IF(NSTP,EQ*0)  go  to  1509 
GO  TO  802 
1509  MSTPat 

GO  TO  749 
802  I*IP1 

IOWN(J)aI 

NTOPa? 

GO  TO  7300 
741  KS7Pa3 

lOWNCJJal 

GO  TO  740 
749  MToPaE 

P(I-l, JJaPC 
GO  TO  7300 
7300  lalOt<NtJ«l)-l 

IF(J.EU.2)  lBl9aN(J) 

N'JPaJ 
GO  TO  739 

C END  OF  corner  node  ADJUSTMENT 
1100  CONTINUE 

730  TF(I,E0,1MM1)  so  TO  1101 

752  CONTINUE 

noi  IF(NuP.Ea.l)  I04N(J)a0 

P(IMAX, J)aP(l, J) 

72U  CONTINUE 

CALL  LKLOAO 

C leakage  and  load  for  this  ITERATION  ARE  TABULATED 
HLOOniDaHLOAD 
00IN(lII)a3lN 
OOeUTdlDaQOUT 
wRirE(6,lSMZX»HuOAO»OJN,OOjr 
pRITEfSf 17} (lOKNCJK), JKal, JMAK) 
lb  F0RNAT(14|3E14,S) 

BESl  AVAIIABIE  COPY 
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iFdii.uT.s)  sa  ra  7io 

tali  ■HLa3(  IT  1*7)  ♦hl0D(1II«6)»hli3o  (111*5)  ♦Ht.aa  ( 1 1 1- 9 ) 

TAtlBfTALI  ♦HL0O(ni*3))/5. 

rALa^CMLnodii-aj^HLaodn-DtHLBodii))/!, 

uyif^AVaDtO 

UQMTAV«Q,0 

Da  13  KFil^S 

AKABKK 

K<K«KA«1 

UUIMAVbOSI  9AV«3QHdn*KKK) 

M OaaTAV«0  jaTAV43QajT(ni-KKK) 

(jaiNAViQJiNAV/AKK 

uai<)TAy«09aTAV/AK< 

UOAVa  (3IS<f3ajT) /2, 

C LkiAD  Ca'yVERSENCE  CHECKED 

lF(AtJS((rAul-TAL2)/TALl),GT,E<»S)  63  70  7(0 

C leakage  iJAUaNCE  checked 

iF(At)s((3QiMAy-a30rAV)/auiNAy) ,le,3.i)  60  ra  2001 

rji)  r.piNTi'JUE 
ITfcRS*!  n 
wRI TE(b,  1) 

I fawMATC'  wa  cayyERCENCE') 

1 aTL()«T  AL2»EHL3AD 

MRITE(bd5}031VAy,  3(307 AV,TarLD,HL0AO«aOAy,0HE&AP 

lb  F?HMAT  ((*fcl  4,b)  • 

wETfJRM 

JJOl  17ERS»II1 

Ta7Lf)«7AL2tEHL0A3 

HHLSaHL0AD/737uD*13O, 

^RI7Efb»17) (inKM( J) » J«l» JMAK) 

17  FaMMA7(ni3) 

MRl7E(b,  la)  l7ERS,33nAV,QO07Ay,Ta7i.O,HLeAD,3aAy,  aiEGAP 
la  FaRMA7i»  iters:', T5»'  SOINAV:  ' ,Eia,(»,  • 330TAy*'»tl  9,b,/, 

1'  TaTLD:',Eia,b,'  HL0AO:',E14,b, ' iSAV: ',Ei4,b, ' aKEGAi'.Eia.b) 

return 

ENT) 


BEST  AVAIUBLE  OPY 
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1 


r i 

i 


SltbRlAJTlMF.  I^<AVSHP 

C SunPRiaGPAM  CAtCJlATES  EXPECTED  </ALUES  AND  0EPIVATIVE8, 

C theta  (T)  MJ8T  BE  SJPPLIEO 

C(!ihh«)N  DT,3R,TB(l02),R6(13)f  liAX,  JNAX,HASP(102) 

C0HH0N  T,R,Rl,R0,PI,P0,PCfP(lO2»  lS),PIE*N««AVE|PFAC,PHtS 
C0MH0N  HNEAN,HtA(tO)  f HIB(IO) » ^ELi  '«>TERHS|  lA  jP,  I AON 
Ci^KMkiN  H,Di,EH,DEH,EH3,DEH3,Et0H3»ai,DBl,ENL0AO#PBAR(132) 

CBMM^n  AIHi (102,13)»AIP1(102(13),AJN1(102|1S),AJP1(102,13) 
BRHS(102,  13)»HU0(102)f  HL0AOf  O2N,Q0Urrr0ruO«E10H 

snoo  M«HMEAN 
r>M«o,o 
!■! 

CCONvi^AVE*! 

CS»C0StCCC*T) 

SN«SIN(CCC*T) 

A (1) tCS 

DH«OH-CCC*HlA(n*SN 

if(h,le,i,j)  ea  ra  as 

I HaH 
DEhsOH 

f h3«h**3^h/3, 
nEH3B(N.aH*H^l ,) aOH/S, 

E|0Hb(1,-H«H)*«3*AL0GC(H^1.)/(H’*2,)) 
t ll(>HB(EiaH*2,BHB(33.»H*M*(15.BH*M-aO,)  ) /15,  J b35,/32# 

E10HiB3,B(l,-HBH)«(5,*HBH-l,) *AL0S(CM^l,)/(H.l,) ) 
E10H3Bl£l0H3-HB(2b,-3O,*MBH) J BiS./Sa, 

BIbI. 

DBl bD, 

Ga  Ta  so 

as  EHB-ij,25«-HBH*t3,l-N*H/S6,) 

EHB(0.12S»HB(lb,i^35,fM»tO,5*HBHBEH))  )*35,/32, 
UEH«ClB»/35,»H*n,tH*H*(»l,*HBHB(,b«H*M/7»)J))*55./32,B0H 
E h3  Bib,/ 35,  25  ( 3» /2B0  • oHAH/BRO.  ) ) ) 

EH3B(,02S*-iBllb,/l35,»HB(,3T5»H*(EH3))))B35./32. 

DEHSbI  ,4>H*H*(-,3fHBHB(3,/3S,-H*H/84,)  ) 
0EH3B(lb./l0S,*HB(,7S«H*(4B,/3S,»HBDEH3)))*3S./32,*!)H 
El0HBf-55,fHB(l32,*HB(345,4M*(*lbO,»Ho(-aO5t»H*(bO,*l47,*H) )) ))) 
E lli!HB(EiaH/bO,*(U»H*H)«*3BAl.aS((U^H)/OCL)  )*3S,732, 

E 10H3b1,E^1O 

HtB(lb.tH*(3S,tHBH*(*3S,^H*HB(21,-S,*HBH)))) /32, 
l)»I»(l.-H*'i)**3*35,/32,*0H 
SO  return 
E'ln 
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SU6l»0UriNE  LXiaAO 

C SUBPR06RAH  CALCULATES  LEAKAGE  AMO  MECHANICAL  A^O  PLJD  RRESSJRE 

c lsao  support 

CflHMiaN  0r,0R,TG(102)|RG(13)  ,IMAX,  JMAX,HASP(132) 

C0MM0N  T,R,Rl,R0,PX,P0,PCrP(132,IS)rPICf NMAVCfPFACfPHLS 
C0HM0N  HNEAN^HlAdO)  »HIB(iQ)  ,OCL»NrERMS,IAjP,lAON 
C0MM0N  H,OH|EH,DEH,EHSf OEHS,E10HS,9I,DBIf EMLaA0(PSAR(l32} 
C0MM0N  AIMl(l02,13)»AIPl(102,l3)f AJM1(i02»1S),AJP1(102»13} 
C0MM2N  BRM9(lO2,13}»HLD(lO2},HL0AOf OlNfQBUTf TBTLOfElBH 
niMENBlBN  EMLD(1Q2) 

HL0AO«O.3 

EHL0AOB3.O 

JST«P*JMAX»1 

D0  100  lalflMAX 

HLU(I) aO.O 

eML0(n"3,3 

PBAR(l)a3.3 

T«.TG(I) 

call  HAV8H» 

D0  SO  Jal«JST0P 

HL0(naHLnCI)^3R»(RG(J)*P(I,  J)tRB(J»n*P(I,  JAl)  W2,/R0 
SO  C0NTINJE 

fcMLD(I)a(l.-an*PfAC 
PHAR(I)aEMtOCI)^HLOCI) 
iF(i,E3.n  r,a  tm  loo 

HLtiA0aMLaA3tMWAVEa(TG(I)-TGtI-l)  )*(HL0(nAHL3(I-l) ) /2,*R0 
fc1LaA0aE'1uaA0*'»WAVE*(EMLf)(n*CML0(I-lJ  J*,5a(rG(n-rS(I-n  ) 
lUO  C0NTIMJE 

00  tOHO  lLCAKil,2 
G0  T0  (tin,120)«lLCAK 
110  R«R1 
Ja2 

OOHb.OR 
UINaO.O 
Ga  T0  133 
IPn  RaR0 

JaJMAX-l 

OORaDR 

QPUTaO.O 

130  00  200  JalflNAX 

220  TaTGd) 

CALL  «tAVSHa 

OPRa(P(I»JAl)-P(I»J-lJ) /2,/OR 

OPORaOPR^03Ra(P(I, J»l)-2,aPfl, J)tP(I, J^1))/0R/0R 

UTNEMa»DP0RaRaCH3aBl 

IFd.EU.l)  60  T0  180 

r,0  TR  (16S|170),ILEAK 

IbS  UINaOIN»SMAVEa(aT8LO«OTNEW)a(rGd)*T6(I«l))/2,/12, 

60  T0  160 
170  C0NTINJE 

U0UTa30iiTtNNAVCa(3TOLD^QTNE4)a(T6(X)«TG(I*l))/2,/12. 

160  UTRLOaSTNEA 
200  C0NTIMJE 
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SUBROUTIi^E  PRSHAR 

SUBPROGRAM  calculates  THE  N-TH  PRESSURE  HARMBRICS  PAl  AM9  PHi 
CBHMON  DT,aR,T8(102) ,RG( 13} , IRAK, JMAK, hasp (132) 

COMMON  T,R,RI,RB,PI,P0,PC«P(lO2f  13)  »PIE*NRAt/E»PFAC|PHLS 
COMMON  HmEAN,HIA(13} ,HIB(10) rOCLf MTERM8*1AjP,XA0R 
COMMON  H,DH»EH,OEH,EH3rOEH3«ElBH3,SZ,OBIf EMLaA0«*SAR(132) 
COMMON  A1H1(102,13)  ,AIP1 (102,  13)f AjMl(l02f 13) I AJ»1 (102,13) 
COMMON  BRHS(102,13) ,HL0(102) ,HL0A0,31N,QBJr,rBTL9,E10H 
DIMENSION  PBA(10),PBB(1Q) 

PBbO.O 

LL«1 

PBA(LL) ■3.9 
P8B(LL)*D,0 
DO  A40  Til,IHAX 
8!U  T2«TG(I) 

ARGbLL*HMA^E*T2 

P8A2BPBAR(n*C0S(ARC) 

PRBPbPBAR(I)*SIN(AR6) 

IF(I.CO.l)  GO  TO  830 

PBA(LL)"»BA(LL)^3.5*(PBA2^P0Al)*(T2-Tl)BNAAyE/PIE 
PBB(LL)bPB3(LL)*3,5b(PBB2*PBB1)*(T2-T1)*NAAVE/PIE 
IFILL.^'E,!)  go  to  830 

P8BPB^0.S*(P8AR(I)fP8AR(I-l))*(T2-Tl)*NHAVE/2,/PIE 
830  TIbT? 

PBAIbPBAE 
PBRl bP8B2 
840  tONTINiJE 

^'RlTE(6,l)aBA(LL),P8B(LL) 

I format C/,5K, 'PRESSURE  H AHMONI CS* ' , 2E 1 4 , b) 

RETURN 

END 
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SUBRUJTIME  ML(tAO(itf,Ha) 

C 8UBPR0GRA*'.  JSES  SECA'JT  METHlflO  Ta  FIND  Ha  THAT  HATCHES  APPl1E3 

c lbao  tb  lbao  sgppaRT 

CaHHaN  0T«3R»TS(102)  »RC(13)  flHAX,  JHAXfHASPdOa) 

CaHMaN  T,R,R1,R8,PI,P8,PC,P(132,  J3)  »PIE,  MHA</E»PFAC,PHLS 

CUMMdN  HMEANr HIA(13) ,HI8(10) , 3ELr HTERHSf  ZA jP; I AOH 

C8HH0N  H,OH,EH,OEH,EH3,OEH3,EiaH3«B1.0BI,EHLaAD,PaAR(132) 

C8HM8H  AlHl (I02f 13)i AIPl (102, 13) ,AJH1 (I0a»l3) ,AJP1 (102, 13) 

C0MH8N  BRHS(102, 13)  ,HL0(102)  ,HLaA0,3X>i,QaUT,  T8TL3,E18H 

DH8i,f)S 

T«L*,005 

IKaO 

H81«Ha 

H0g*M0»OHB 

H83b0,0 

F3«D,0 

jter»o 

na  20  i«i,io 
TFdK.EO.l)  50  za  33 
ZF(H83.E3.H02)  68  T8  32 
MHEAH»h8J 
^RITE (b,23)Mai 
23  FPRMATISX.FID.S) 

CALI.  E3NS8Ld.Z,  ,331, 130,  ITERS,  J TER) 

JTE»»JTER^1 

IF(TTERS,EJ,103)  58  T8  27 
JHHl*JMAX-l 
2t*  Pl«TaTLO-w 

IF(ABS(F1/W) .UT.Tat)  60  T0  13 
IF  (H23,Ea,H0n  68  T0  33 
32  HMtAN«M02 

wRITE(b,23)Ha2 

call  EOHSPtd.7,,031,  !Of),ITERS,  JTER) 

JTtR»JTER*l 

IF(ITERS,E3,103)  sa  Ta  27 
F2«T0TLD-w 

IF(ABS(F2/H),LT,Tau)  60  T0  13 
33  H03»Hai-Fl*fHa2-Hai)/(F2-Fl) 

IF(Fl*F2)2,2,3 
g HMEAHaH03 
IF«1 

wRlTE(6,23)Ha3 

CALL  EQNS0  yd.  7,  ,331,130,  ITERS,  J TER) 

JTERi JTER^I 

lFdTERS,E3,100)  52  Ta  27 
F3*T0TLD»W 

IF(A8S(F3/X) .LT.TBl)  60  T0  13 
IF(FS*F2)31,31,S3 

3 IF(ABS(H03-H01),LE.ABS(H03*Ha2))  68  78  3l 

I«»0 

30  H02aHM3 
F?aF3 
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i\ 

f l»f3 

?0  CJINTIMUE 

wRITE(t»,ll) 

11  F3HMAT(3>«, 'N0  CB'^VERGENCe  IN  4JL0AO*) 

la  wNlTE(b,l3) JTER 

13  F0kMATl5ii, 'HNEAN  C0MVERGENCE-I TE»A TI aNS» • # 1 5) 

^-RITE(6»  12)H*1EAN,  TZTLOfPHtS 

la  F0RMAT(3K, *HNEAN|*,Eia,b,3X, 'raTL0l*#Eia,6,3A, 'PHLJl'fElA.b) 

27  WRlTEIb.n 
1 F0RM/IT(1HI,/) 

100  00  21  Ilal,INAX 
HASP(rn«0,0 
T«T5(II) 

CALL  WAV5H» 

IFlH,LT,l,)HASP(II)«l,D 

21  ^RlTE(6,22)n,  (P(II,JJ)»JJ»l»  J'<AX),HASP(II) 

22  ► 0RMaT(I4,  lOF7,a,  ax,  ID 
2H  RETURN 

E.NO 
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SUHR.dJTI>4C  Fp(ICT(RSBAR,R(40C,RF,  AMJ) 

C SUBPRtrtGRAM  CAtCJLATES  FRICTIBN 

C0MM(9N  r)T,3R,TS(102)  ,RG(13)  , 1 1 AX  , J AX , HASP  ( 1 32) 

C0MM0W  T,R,RI,  R0,PI,R0,PC,P(132,13)»PIE»N<<AVE»RFaC#<»hlS 
C0HM0N  HHEAN,HIA(13) »3EUf NTERMSf lAjP.lAOM 
C0MM0N  H,U4,EH,DEH,EH3,DEH3,ElEH3,3I,0Bl,EHL3A0«»aAR(t  32) 
CBMSBN  AIMl (I02,)3)f AIP1(102,1B), AJMl(t02f 13) ,AJP1 (102.13) 
C0MM0N  BRHsrio2, 13)  >hu3(102)  »>4uflAOf  oiN,Qeur»rerL3f  El  an 
F2«0,n 

00  10  jBltJHAX 
rsuMso. 

RbRG(J) 

NNal 

00  20  iBlflHAX 
G0  T0(3O|43)|NS 
30  IF(P(I,J).EO,PC)  62  TB  60 
G0  T0  40 
60  tiTG(I) 

CAUL  4AW5HP 
I Hl«EH 
NMb2 

40  T«TC(I) 

CALC  -^AVSmP 
lil«l-l 


70 

«0 


VO 

20 


1 on 

If) 


IP1»I^I 

lF(IPl,GT,IMAX)IPlt2 
1F(IMI,LT,1)IM1«I4AX-1 
IF(PtI,J),£0,PC)  SB  70  70 
|)P.(P(IP1,  Jl-PCm,  J))*,5/DT 
F «R*R» (R*E10H*,5/R*OP*EH) 401/RF 
G0  T0  80 

Fa8T*R*R4R*E10H*El41/EH/HF 
F«F*R0aC/PF* (1 .-SI) *P5BAR*R*R 
FNEMmF 

IFCI.EQ.l)  G0  70  93 
TSuMbTSU'1*(F4E-*^F0lD)  *,S*DT 
F0LO«Fi4E4 
CPNTI'^JE 
F nbTSJM 

IFtJ.EO.l)  G0  70  130 

FZ«F2*(F4PF0)4,5*OR 

F0«FM 

C0N7IVJF 

A»1l,l»FZ/TV)7tD/Ra0C*4PAVE 

RE7ijRN 

end 
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REOJIRED  LiSAO  M*0,2a3tt3HE-02 
NXAVE  3 

R2  n.SO0OE*O2 
HIaCU  0,2030E  00 
RX  O.ROOQE  no 
R/C  O.tOaOE  06 
RM  O.SOOOC  00 


0,91500 


1 

0,179433E«02  *0 •4467S5E»01 

-0,2402596-01 

0 

46  36  32  29  0 

0 0 0 

0 

2 

0,1786546-02  -0,38960! E-Ol 

-0. 2085966-01 

0 

52  45  35  3)  23 

0 0 0 

0 

s 

0,1761466-02  -0.349013E-31 

-3.184523E-01 

0 

57  51  44  34  24 

0 0 0 

0 

« 

0,1720726-02  -0,36a421E-Dl 

-0,8931 12E-02 

0 

57  53  46  39  25 

0 0 0 

n 

5 

0,1667166-02  -0.305526fc-01 

-0,2735686-01 

0 

61  56  SI  42  26 

n 0 0 

0 

6 

3,165BOOE-02  -0 . 327 1 666- J1 

-0,2594996-01 

0 

61  57  52  43  27 

3 0 0 

0 

? 

0,164S52E-02  -0,3121506-01 

-0,2814396-01 

3 

62  59  53  43  26 

3 3 0 

0 

6 

0,1638416-32  -0,2665366-31 

-0,2471666-01 

0 

64  60  54  43  0 

0 3 3 

0 

9 

0,162489E-02  -D,2827B3E-31 

-0.2730946-01 

0 

65  61  54  42  0 

0 0 0 

0 

10 

0.1618546-02  -0.2B3H746-31 

-0. 2775396-01 

0 

65  61  S3  41  0 

3 0 0 

0 

11 

0,1615016-02  -0,2829a4e-31 

-3,2738246-01 

0 

65  61  51  40  0 

0 0 0 

n 

12 

0,1613006-02  -0,2638706-31 

-0,2880546-01 

0 

65  61  49  39  0 

3 3 0 

0 

13 

0,1614946-02  -0,2832766-31 

-0.2763326-01 

0 

65  60  48  36  0 

0 3 3 

0 

19 

0,1616376-02  -0,289Sb36-31 

-0,2811926-31 

3 

65  58  46  37  0 

3 0 3 

0 

15 

0,1620956-02  -0,2972136-31 

-3. 2837986-01 

0 

63  57  uH  36  0 

0 3 0 

0 

16 

0,1627186-32  -0,3073396-01 

-0,2843906-01 

0 

63  57  49  39  0 

0 3 3 

T 

17 

0.1632326-32  -0.3060706-01 

-0.2854376-01 

0 

63  57  50  40  0 

0 0 0 

0 

Is 

0,1637016-02  -0,3156306-31 

-0,2906276-01 

0 

62  57  50  41  0 

0 3 0 

0 

19 

0,1640706-02  -n,312M76E-31 

-0,2774436-01 

0 

62  50  51  41  0 

0 3 3 

3 

20 

0,1640316-02  -0,3!144n3fc-01 

-0,2817906-01 

0 

63  59  51  41  0 

3 0 0 

3 

21 

0,1639636-02  -0,3112776-01 

-0.2733586-01 

0 

63  59  51  40  0 

J 0 3 

0 

22 

0,1637176-02  -0,3035j 46-01 

-0,2737426-01 

0 

63  60  51  39  0 

0 3 0 

0 

23 

0,1634176-02  -0,2951426-01 

-0,2723996-01 

0 

64  60  50  39  0 

0 0 0 

0 

0 

64  60  50  39  0 

0 0 0 

0 
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